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EVALUATION  OF  ALTERNATIVE  PROCEDURES  FOR 
ATMOSPHERIC  ABSORPTION  ADJUSTMENTS 
OURING  NOISE  CERTIFICATION 

VOLUME  II:  COMPUTER  PROGRAM 
1 .  INTRODUCTION 

This  Volume  contains  the  source-code  statements  for  the  computer 
program  that  was  prepared  to  perform  the  calculations  for  the  analyses 
described  in  Volume  I.  The  program  was  written  in  the  extended  FORTRAN  IV 
programming  language  for  compilation  and  execution  on  a  large-scale  digital 
computer.  Input  data  can  be  entered  from  an  interactive  terminal,  digital 
magnetic  tape,  or  punched  cards  to  a  disk  file  for  access  during  execution. 
Results  can  be  output  to  a  line  printer  or  to  an  interactive  terminal. 

The  objective  of  the  program  is  to  calculate  adjustment  factors  to 
be  applied  to  test-day  1/ 3-octave-band  sound  pressure  levels  from  aircraft 
flyover  noise  measurements  to  account  for  differences  in  atmospheric 
absorption  along  the  sound  propagation  path  under  test  and  acoustical- 
reference  conditions.  The  band  adjustment  factors  are  calculated  by 
four  different  procedures. 

To  achieve  the  objective,  the  program  reads  in  measured  sound  pressure 
levels  at  0.5-s  intervals.  From  the  measured  sound  pressure  levels, 
the  program  calculates  corresponding  values  of  perceived  noise  levels 
(PNL),  tone-corrected  perceived  noise  levels  (PNLT),  and  A-weighted 
sound  levels  (AL) .  The  sets  of  PNLs ,  PNLTs ,  and  sound  levels  are  searched 
to  find  the  maximum  values:  PNLM,  PNLTM,  and  ALM.  Effective  perceived 
noise  level  (EPNL)  is  calculated  from  the  set  of  PNLT  values;  sound  exposure 
level  (SEL)  is  calculated  from  the  AL  values.  The  range  over  which  the 
PNLT  or  AL  values  are  summed  to  find  EPNL  or  SEL  is  determined  by  a 
specification  of  th*e  10-dB-down  values  from  PNLTM  and  ALM,  respectively. 

The  time  of  occurrence  of  PNLTM  and  of  ALM  during  the  flyover  defines  the 


test-day  spectra  and  propagation  paths  to  be  used  in  determining  band  adjust¬ 
ment  factors. 

To  perform  the  atmospheric-absorption  adjustments,  the  program  reads 
an  input  data  file  for  air  temperatures  and  relative  humidities  at  various 
heights  above  ground  level.  Barometric  pressure  is  also  read  in.  The 
air  temperatures,  relative  humidities,  and  barometric  pressure  represent 
values  that  have  been  previously  interpolated  to  the  time  of  the  flyover 
noise  measurement. 

After  the  test-day  spectra  have  been  adjusted  to  reference-day  meteorolo¬ 
gical  conditions  by  the  four  procedures,  PNL,  PNLT,  and  AL  values  are  calculated 
from  the  reference-day  sound  pressure  levels.  Reference-day  EPNLs  are  found 
from  the  reference-day  PNLTMs  and  reference-day  SELs  are  found  from  the 
reference-day  ALMs  by  adding  the  corresponding  values  of  test-day  duration- 
correction  factors. 

As  output,  the  program  lists  the  measured  test-day  sound  pressure 
levels  and  corresponding  PNLs,  PNLTs,  and  ALs,  EPNLs ,  and  SELs.  Meteorological 
data  are  also  listed.  The  adjusted  reference-day  sound  pressure  levels  and 
corresponding  PNL,  PNLT ,  EPNL,  AL,  and  SEL  values  are  tabulated  on 
separata  pages.  The  program  interrogates  the  input  data  file  for  another 
set  of  data  until  the  data  are  exhausted. 

The  program  was  written  as  one  main  routine  called  TESTREF  with  nine 
subroutines.  There  is  one  labeled  common  block  of  literal  data  used  for 
the  legend  on  the  final  two  output  pages.  The  main  routine  reads  in  the 
input  data,  calls  the  subroutines  to  perform  the  calculations,  and  controls 
and  specifies  the  output  data  listings. 

It  is  noted  that  while  the  program  is  relatively  efficient  (requiring 
on  the  order  of  only  30  seconds  of  central-processor-unit  time  for  all  the 
calculations  for  one  set  of  flyover  noise  measurements) ,  greater  emphasis 
was  placed  on  producing  an  operational  program  that  met  the  objective  of 
the  program  than  on  maximizing  computing  efficiency.  The  program  statements 
include  a  generous  number  of  descriptive  comments  to  assist  the  user  in 
understanding  the  various  steps. 
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The  balance  of  this  Volume  contains  four  additional  Sections.  Section 


2  describes  and  defines  most  of  the  parameters  and  variables  used  in  the 
program  along  with  their  accompanying  symbols  and  preferred  units.  Those 
parameters  and  variables  not  included  in  Section  2  are  unique  to  a 
subroutine  and  are  defined  in  the  program  statements  for  the  subroutine. 
Section  2  also  omits  definitions  for  standard  intrinsic  FORTRAN  sub¬ 
programs  such  as  ABS  for  the  absolute-value  function,  COS  for  the  cosine 
function,  INT  for  the  integer-part  function,  and  LOGIO  for  the  function 
to  calculate  base-10  logarithms.  Definitions  for  ordinary  DO- loop  index 
counters  are  not  included,  nor  are  standard  FORTRAN  logical  or  arithmetic 
functions  such  as  .LE.  for  less  than  or  equal  to,  .GT.  for  greater  than, 
and  so  on. 

Section  3  contains  the  listing  of  all  the  program  statements  for  the 
main  routine  and  all  nine  subroutines.  There  are  a  total  of  2024  lines  of 
source-code  for  the  main  routine  and  the  subroutines. 

Sections  4  and  5  present  the  data  for  a  sample  or  test  case  from  one 
of  the  cases  used  for  the  analyses  in  Volume  I.  Section  4  contains  the 
input  data  in  card-image  format.  Section  5  contains  the  corresponding 
listings  of  output  data. 


2.  PARAMETERS  AND  VARIABLES  USED  IN  MAIN  PROGRAM  AND  SUBROUTINES 

Symbols  Description  and  Units 

1  height,  in  meters,  of  the  airplane  reference  point  above  the  ground 

plane  when  over  the  microphone.  AH  is  read  from  the  input  data 


AHSAM, 

AHSPTM 


ALMB, 

ALMC, 

ALMD, 
ALME 


ANSAB 


array  of  the  average  of  the  relative  humidities,  in  percent, 
measured  at  two  consecutive  heights  above  ground  level.  The 
average  relative  humidites  in  AHR  correspond  with  the  average 
heights  of  the  layers  of  the  atmosphere  in  array  HMI. 

difference  between  the  height  of  the  airplane  and  the  height  of 
the  microphone  at  the  time  when  (a)  for  AHSAM,  the  maximum 
A-weighted  sound  level,  and  (b)  for  AHSPTM,  the  maximum  tone- 
corrected  perceived  noise  level  was  noted  at  the  microphone,  m 

array  of  the  A-weighted  sound  levels  at  0.5-s  intervals  corres¬ 
ponding  to  the  measured  test-day  1/ 3-octave-band  sound  pressure 
levels  at  0.5-s  intervals  during  a  flyover,  dB 

maximum  value  of  A-weighted  test-day  sound  levels  in  array  AL 
(note  that  if  there  is  more  than  one  sound  level  with  the  same 
maximum  value,  then  ALM  is  the  one  that  occurs  last),  dB 

A-weighted  sound  levels  corresponding  to  the  reference-day 
spectra  SPLAB,  SPLAC,  SPLAD,  and  SPLAE  at  the  time  of  occurrence, 
at  the  microphone,  of  ALM,  dB 

airplane  flight  Mach  number  when  over  the  microphone.  AM  is 
read  from  the  input  data  file. 

distance  between  the  airplane  reference  point  and  the  microphone 
at  the  time  the  airplane  was  overhead  (equal  to  the  difference 
between  airplane  height  above  the  ground,  AH,  and  microphone 
height  above  the  ground,  ZM) ,  m 

name  of  a  subroutine  for  calculating  pure-tone  atmospheric 
absorption  coefficients,  in  dB/m,  by  the  method  of  American 
National  Standard  ANS  SI. 26-1978,  given  the  air  temperature,  in 
degrees  kelvin,  relative  humidity,  in  percent,  air  pressure. 


in  standard  atmospheres,  and  the  frequency  of  the  sound  wave  in 
hertz 


AP 


APMB , 

APMC, 

APMD, 
APME 

APNL 


AR 


ARP866A 


AS 


ratio  of  the  A-weighted  squared  sound  pressure  to  the  square  of  the 
reference  pressure  in  a  calculation  of  A-weighte  sound  level,  in 
array  AL,  from  a  set  of  test-day  1/ 3-octave-band  sound  pressure 
levels 

ratios,  as  for  AP  above,  except  used  in  the  calculations  of  ALMB, 
ALMC,  ALMD,  and  ALME  from  reference-day  spectra  SPLAB,  SPLAC, 

SPLAD,  and  SPLAE 

name  for  a  dummy  variable  used  in  the  call  for  subroutine  CPNL 
when  calculating  the  array  of  test-day  perceived  noise  levels 
corresponding  to  the  1/ 3-octave-band  sound  pressure  levels,  SPL(I,J), 
and,  specifically,  to  the  SPL  spectrum  ASPL(J)  for  a  particular 
time  index  I  at  some  0.5-s  interval  after  the  first  set  of  test- 
day  data 

in  method  E  for  adjusting  the  test-day  sound  pressure  levels  to 
reference-day  conditions,  AR  is  the  pure- tone  atmospheric  sound 
absorption  coefficient,  in  dB/m,  returned  by  the  call  for  subroutine 
ANSAB  for  the  defined,  uniform,  reference  meteorological  condi¬ 
tions  with  a  reference  air  temperature  of  TKR  *  298.15  K,  a 
reference  relative  humidity  of  HRR  -  70.0  percent,  a  reference  air 
pressure  of  1.0  standard  atmosphere,  and  at  the  exact  center 
frequency,  FC  in  hertz,  of  a  1/3-octave  band 

name  of  a  subroutine  for  calculating  pure-tone  atmospheric- 
absorption  coefficients,  in  dB/m,  applicable  to  1/3-octave  band 
sound  pressure  levels,  according  to  the  method  of  SAE  Aerospace 
Recommended  Practice  ARP866A,  15  March  1975,  given  the  1/3-octave- 
band  center-frequency  index  counter  J  (where  J  *  1  to  24  for 
nominal  center  frequencies  between  50  and  10,000  Hz),  the  air 
temperature  in  degrees  celsius,  and  the  relative  humidity  in 
percent 

airplane  airspeed,  in  m/s,  when  over  the  microphone.  AS  is  read 
from  the  input  data  file. 


•  .#v  v  •v- 


A3PL 


ATE 

ATK 

ATTAC 

ATTPTC 

AVAL FAC 


array  of  1/ 3-octave-band  sound  pressure  levels  used  with  APNL  in 
the  call  for  subroutine  CPNL  to  calculate  the  test-day  perceived 
noise  levels  at  0.5-s  intervals  corresponding  to  the  test-day 
sound  pressure  levels,  SPL(I,J),  dB 

array  of  average  air  temperatures  over  horizontal  layers  of  the 
atmosphere,  with  a  layer  defined  by  two  consecutive  heights  above 
ground  level  where  test-day  meteorological  data,  appropriate 
for  the  time  of  the  aircraft  noise  measurement,  were  measured,  C 

in  adjustment  method  E,  pure-cone  atmospheric  sound  absorption 
coefficient,  in  dB/m,  as  in  AR  above  but  for  test-day  meteorological 
conditions  along  segments  of  the  sound  propagation  path  defined 
by  the  heights  where  meteorological  data  were  measured 

an  array  like  ATC  above  except  for  air  temperatures  in  degrees 
kelvin 

in  adjustment  method  C,  attenuation,  in  dB,  for  test-day  meteorolo¬ 
gical  conditions  over  a  segment  of  the  sound  propagation  path  as 
well  as  over  the  entire  pathlength  PDAM  corresponding  to  the 
relative  data-sample-time  TRAM  associated  with  the  spectrum 
SPLAM  at  the  time  of  occurrence  of  the  maximum  test-day  A-weighted 
sound  level  ALM  and  with  atmospheric  absorption  calculated  by 
Che  method  of  SAE  ARP866A 

attenuation,  in  dB,  like  ATTAC  above,  except  for  the  propagation 
pathlength  PDPTM  corresponding  to  the  relative  time  TRPTM 
associated  with  SPLPTM  at  the  time  of  the  maximum  test-day  tone- 
corrected  perceived  noise  level  PNLTM 

in  adjustment  method  C,  average  pure-tone  atmospheric  sound  absorp¬ 
tion  coefficient,  in  dB/m,  applicable  to  a  1/3-octave  band  and 
defined  by  the  ratio  of  the  attenuation  to  the  propagation  path- 
length  (i.e.,  the  ratio  ATTAC/PDAM)  and  used  with  the  atmospheric 
absorption  coefficient  RALPHA  (calculated  also  by  the  method  of 
ARP866A  but  for  the  defined,  uniform,  reference  meteorological 
conditions)  to  determine  a  band-adjustment  factor  to  be  added  to 
the  test-day  sound  pressure  levels  SPLAM  at  the  time  of  ALM  to 
yield  the  reference-day  spectrum  SPLAC 


AVALFPTC 


like  AVALFAL  above,  except  calculated  from  the  average  attenuation 
ATTPTC  and  propagation  distance  PDPTM  to  determine  the  band-adjust¬ 
ment  factor  to  be  added  to  the  test-day  sound  pressure  levels 
SPLPTM  at  the  time  of  PNLTM  to  yield  the  reference-day  spectrum 
SPLPTC 

AW  at  the  nominal  center  frequencies  of  the  1/3-octave  bands  between 

50  and  10,000  Hz,  a  data  array  of  24  frequency-response  weighting 
factors  relative  to  the  response  at  1000  Hz  for  the  A-weighting 
scale  from  the  International  Standard  on  sound  level  meters,  dB 

BAAD  array  of  24  band-adjustment  factors  to  be  added  to  the  1/3-octave- 

band  test-day  sound  pressure  levels,  SPLAM,  at  the  time  of  ALM  to 
determine  the  reference-day  spectrum,  SPLAD,  according  to  adjustment 
method  D,  dB 

BAAE  similar  to  BAAD  above,  except  that  BAAE  is  used  to  determine  the 

reference-day  spectrum,  SPLAE,  according  to  adjustment  method  E, 
and  except  that  BAAE  is  not  an  array  but  a  variable  determined  anew 
for  each  of  the  24  bands,  dB 

BAPTD  an  array  like  BAAD  above,  except  that  the  band-adjustment  factors 

are  to  be  added  to  the  test-day  spectrum,  SPLPTM,  at  the  time  of 
PNLTM  to  determine  the  reference-day  spectrum,  SPLPTD,  according  to 
adjustment  method  D,  dB 

BAPTE  similar  to  BAPTD  above,  except  that  BAPTE  is  used  to  determine  the 

reference-day  spectrum,  SPLPTE,  according  to  adjustment  method  E, 
and,  like  BAAE,  except  that  BAPTE  is  not  an  array  but  a  variable 
determined  anew  for  each  of  the  24  bands,  dB 

BP  barometric  pressure  of  the  air  (usually  at  a  height  of  10  m  above 

ground  level)  and  read  from  the  file  of  input  data,  in.  Hg 

C  array,  one  for  each  set  of  test-day  1/3-octave-band  sound  pressure 

levels  at  0.5-s  intervals,  of  the  maximum  values  of  the  tone-correc¬ 
tion  factors,  in  dB,  calculated,  according  to  the  rules  of  the  3 
April  1978  amendment  to  FAR  Part  36,  by  procedures  in  subroutine 
P36TC  when  P36TC  is  passed  a  set  of  test-day  sound  pressure  levels 
in  dummy  array  ASPL(J) .  Subroutine  P36TC  returns  the  maximum  tone- 


CMX 


CMXB, 
C.MXC , 
C>!XD, 
CMXE 

C^NL 


DA-'! 


DATA10 


correction  factor  C(I)  for  the  spectrum  ASPL(J)  at  each  relative¬ 
time  index  counter,  I. 

CMX  is  the  name  used  in  subroutine  P36TC  for  the  maximum  tone- 
correction  factor,  in  dB,  in  a  calculation  of  tone-corrected 
perceived  noise  level  for  any  set  of  sound  pressure  levels 

maximum  tone-correction  factor  returned  by  subroutine  P36TC 
when  passed  the  spectra  SPLPTB,  SPLPTC,  SPLPTD,  and  SPLPTE, 
respectively,  dB 

name  of  a  subroutine  that  accepts  an  array  of  2U  1/ 3-octave-band 
sound  pressure  levels  at  nominal  band  center  frequencies  from 
50  to  10,000  Hz  and  calculates  the  corresponding  perceived  noise 
level,  PNL,  using  data  constants  from  SAE  ARP865A,  15  August  1969 

array  of  the  lengths,  in  meters,  of  the  segments  along  the  sound 
propagation  path  from  the  aircraft  reference  point  to  the  micro¬ 
phone  at  the  time  of  the  maximum  A-weighted  sound  level,  ALM, 
for  test-day  atmospheric  conditions.  Segment  lengths  are  deter¬ 
mined  from  the  thickness  of  the  horizontal  layers  of  the  atmosphere 
defined  by  two  consecutive  heights  above  ground  level  where 
test-day  meteorological  data  were  measured,  by  the  sound  propaga¬ 
tion  angle  relative  to  the  flight  path  direction  ip.  at  the 

time  of  ALM,  and  by  the  flight  path  angle,  y,  relative  to  the 
horizontal.  These  angles  are  called  PSIAM  and  GAMAR  in  the  program 
statements  and  in  this  symbol  list. 

subroutine  INTEG  numerically  integrates,  over  a  specified  duration, 
a  set  of  specified  input  data,  in  dB,  that  vary  with  time.  The 
input  data  can  be  any  desired  physical  variable,  but,  in  program 
TESTREF,  subroutine  INTEG  is  used  to  numerically  integrate,  by 
a  summation  process,  (1)  the  tone-corrected  perceived  noise  level, 
PNLT,  to  determine  the  effective  perceived  noise  level,  EPNL,  and 
(2)  the  A-weighted  sound  level,  AL,  to  determine  the  sound  exposure 
level,  SEL.  The  limits  of  the  numerical  integration  are  specified 
by  the  value  of  DATA10  which  is  10.0  dB  less  than  the  corresponding 
maximum  value  of  the  input  data,  i.e.,  DATA10  »  PNLTM  -  10.0  or 
ALM  -  10.0,  in  dB. 


9 


DATE  a  number,  defined  to  be  an  integer  (not  a  real)  constant  that  is  read  j 

from  the  file  of  input  data  giving  the  day  of  the  month  when  the 
flyover  noise  measurements  were  recorded.  DATE  is  used  for  label 
information  on  the  output  pages  for  the  test-day  sound  pressure 

levels.  ] 


DB  DB  is  the  distance  measured  along  the  aircraft's  flight  path  between 

the  overhead  point  and  the  position  of  the  aircraft  reference  point  \ 
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DCF 

DM 

DO  AM 

DOPTM 

DPTM 

DSEL 

EPNL 

EPNLB, 

EPNLC, 

EPNLD, 

EPNLE 


at  the  point  of  closest  approach  to  the  microphone,  m 

duration  correction  factor  in  the  formulation  of  EPNL  and  determined 
for  the  test-day  data  from  DCF  =  EPNL  -  PNLTM,  dB 

minimum  distance  between  the  microphone  and  the  aircraft  flight  path 
(at  the  closest  point  of  approach) ,  m 

distance,  in  meters,  along  the  flight  path  from  the  overhead  point 
to  the  position  of  the  aircraft  reference  point  at  the  relative 
time  TRAM  corresponding  to  the  maximum  test-day  A-weighted  sound 
level,  ALM 

distance,  in  meters,  like  DOAM  above,  except  at  the  relative  time 
TRPTM  of  the  maximum  test-day  tone-corrected  perceived  noise 
level,  PNLTM 

like  DAM  above,  an  array  of  segment  lengths,  in  meters,  along  the 

sound  propagation  path,  except  at  the  time  of  PNLTM  and  at  propagation 

direction  angle  „  or  PSIPTM 
PNLT, max 

duration  of  the  sound  exposure  level,  SEL,  and,  like  DCF  above  for 
EPNL,  the  difference  between  SEL  and  ALM,  dB 

effective  perceived  noise  level,  in  dB,  as  calculated  by  subroutine 
INTEG  from  test-day  PNLT  data  with  the  corresponding  value  of  PNLTM 
as  DATA10  and  with  specified  values  for  IFP,  ILP,  and  NS 

effective  perceived  noise  level,  in  dB,  calculated  from  PNLTMB  +  DCF, 
PNLTMC  +  DCF,  PNLTMD  +  DCF,  and  PNLTME  +  DCF,  respectively,  where 
DCF  is  the  duration-correction  factor  for  the  test-day  PNLT  values 


FI 


FMX 


FMXB, 

FMXC, 

fmxd, 

FMXE 

FREQ 


GAMAD 


GAMAR 


exact  center  frequency,  In  Hz,  of  a  1/ 3-octave-band  filter.  FC  is 
calculated  from  using  the  international  standard  band 

numbers  (ISBN)  that  range  from  17  for  the  50-Hz  band  to  40  for  the 
10,000-Hz  band. 

array  of  nominal  center  frequencies,  in  kHz,  of  the  1/3-octave 
bands  that  yielded  the  largest  tone-correction  factors,  C(I), 
for  the  test-day  sound  pressure  levels,  SPL(I.J).  For  each  SPL 
spectrum,  the  band  level  that  produced  C(I)  is  identified  by 
frequency-band  index  counter  JMX  in  order  to  locate  the  nominal 
band  center  frequency  FMX  in  the  data  array  FREQ. 

nominal  center  frequencies  for  the  bands  producing  the  reference- 
day  tone-correction  factors  CMXB  to  CMXE,  kH2 


a  data  array  consisting  of  the  24  nominal  center  frequencies,  in 
Hz,  of  the  1 /3-octave-band  filters  and  ranging  from  50  to  10,000 
Hz.  The  frequencies  in  array  FREQ  are  declared  to  be  integer 
constants. 

angle  y,  or  gamma,  in  degrees.  Angle  y  is  the  acute  angle  the 
aircraft  flight  path  makes  with  the  local  horizontal.  For  calcu¬ 
lation  purposes,  the  engine  thrust  axis  is  assumed  to  be  parallel 
to  the  flight  path.  Angle  y  is  positive  for  climbing  flight,  nega¬ 
tive  for  descending  flight,  and  is  measured  relative  to  the  direction 
of  flight. 

angle  y  in  radians 


GROUND  array  of  literal  characters  read  from  the  input  data  file  that  gives 

a  description  of  the  type  of  ground  surface  that  was  under  the  mic¬ 
rophone  at  the  time  of  the  test  (e.g.,  asphalt  or  concrete).  The 
description  is  used  as  label  information  for  the  output  pages 
containing  the  test-day  sound  pressure  levels. 

HM  array  of  heights,  in  meters,  above  ground  level  where  the  input 

meteorological  data  were  measured.  Array  HM  is  read  from  the 
input  data  file. 
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HtlEA 


HMI 


HR 


HREF 


HRR 

HR10 


I  AM 


IE,  IS 


array  of  the  measured  or  test-day  values  of  the  atmospheric  humidity, 
expressed  as  the  molar  concentration  of  water  vapor  in  percent,  and 
determined  at  the  average  heights  in  array  HMI  for  the  average 
test-day  temperatures  in  array  ATK,  the  average  relative  humidities 
in  array  AHR,  and  an  air  pressure  of  1.0  standard  atmosphere. 

At  each  increment  of  height,  the  humidity  values  in  the  array  HMEA 
are  calculated  by  calling  subroutine  MOLAR. 

array  of  heights  that  are  midway  between  two  consecutive  heights 
where  meteorologicsl  data  were  measured,  (i.e.,  the  average  height 
above  the  ground  plane  for  each  layer  of  the  atmosphere),  m 

array  of  test-day  relative  humidities,  in  percent,  measured  at  the 
heights  in  array  HM  and  read  from  the  input  data  file 

array  of  atmospheric  humidities,  expressed  as  the  molar  concentration 
of  water  vapor,  in  percent,  like  the  humidities  in  array  HMEA  and  at 
the  same  heights,  but  calculated  for  defined,  reference  lapse  rates 
of  air  temperature,  relative  humidity,  and  pressure  rather  than  the 
measured  test-day  temperatures  and  relative  humidities  and  a 
constant  air  pressure  of  1.0  standard  atmosphere.  Subroutine  MOLAR 
is  called  to  calculate  the  values  of  HREF  at  each  increment  of  height. 

reference-day  relative  humidity  for  a  uniform  atmosphere  and  set 
equal  to  70.0  percent 

relative  humidity  for  the  test-day  conditions  interpolated,  or 
read  directly,  from  the  measured  vertical-profile  meteorological 
data  for  a  height  of  10.0  m,  percent. 

index  from  the  relative  test-time  indexes,  I,  in  the  array  of  relative 
test  times  TT(I)  for  that  value  of  the  A-weighted  sound  level 
selected  from  the  set  of  sound  levels,  AL(I),  and  identified  as  the 
maximum  A-weighted  sound  level,  ALM,  i.e.,  AL(I)  ■  AL(IAM)  *  ALM. 

If,  in  the  search  from  the  first  to  the  last  relative  time  (over 
the  index  range  from  I  ■  1  to  I  *  NS) ,  more  than  one  AL  value  has 
the  same  maximum  value,  then  IAM  is  set  equal  to  the  last  value  or 
the  last  time  when  AL  ■  ALM. 

relative-time  index  counters  used  for  identifying  and  printing  out 
the  columns  of  test-day  sound  pressure  levels,  perceived  noise  levels, 


■*  ■ 
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1 


IFA ,  IFP , 
ILA,  ILP 


INTEC 


tone-corrected  perceived  noise  levels,  and  FMX  values.  Initially, 
IS  —  1  and  IE  *  14  to  print  out  14  samples  at  0.5-second  intervals 
on  the  first  page  corresponding  to  relative  tesc  times  TT(1)  »  0.0 
to  TT(14)  ■  6.5.  IS  and  IE,  for  I  Start  and  I  End,  are  incremented 
by  appropriate  values  until  all  data  samples  have  been  listed, 
i . e . ,  at  IE  »  NS . 

re lative-time  index  counters  determined  by  subroutine  INTEG  for 
the  first  and  last  times  when  the  test-dav  A-weighted  sound  level 
AL(I)  is  closest  to  ALM  -  10.0  (l.e.,  IFA  and  ILA),  or  when  the 
test-dav  tone-correc ted  perceived  noise  level  PNLT<I)  is  closest 
to  PNLTM  -  10.0  (i.e.,  IFP  and  ILP).  The  index  counters  ire  used 
Co  denoce  the  relative  times  on  the  output  page  for  beginning  and 
ending  times  in  the  calculation  of  duration  correction  factors  far 
EPNL  and  SEL;  DCF  is  calculated  for  times  ranging  from  TT<  IFP)  to 
T T(ILP)  and  DSEL  is  for  times  from  TT(IFAj  to  TTMLA)  where  TT ;  I' 
the  array  of  relative  test  times  at  0.5-s  intervals. 

name  of  a  subroutine  that  performs  a  numerical  integration,  bv  a 
simple  summation  of  rectangular  ipproximat ions ,  of  the  square 
of  a  quantity,  derived  from  an  array  of  input  values  in  decibels, 

over  a  range  of  relative  test  times  which  are  the  first  md  the 

lasc  times  when  the  input  data  are  closest  to  DATA10.  The 

answer  returned  by  INTEG  is  called  SUMDAT  for  the  power  sum  of  the 

input  data. 


ft 


I  PM, 
IPTM 


ISBN 


Like  TAM  above,  except  the  reLative  time  indexes  for  that  value  or 
PNL(I)  which  is  set  equal  to  PNLM  and  for  chat  value  of  PNLTM)  which 
is  set  equal  Co  PNLTM 

international  standard  band  numbers  for  1 / 3-octave-band  filters, 
see  FC  above.  ISBN  is  defined  to  be  a  real  (not  in  integer) 
var lab le . 


JMX . 

JMXB, 

JMXC, 

JMXD, 
JMXE 


frequency  band  index  counters  returned  by  subroutine  P3bTC  in  a 
calculation  of  tone-corrected  perceived  noise  level  and  used  to 
determine  the  nominal  center  frequency  of  the  1/3-octave  band  that 
produced  the  largest  tone-correction  factor  1C(I)  or  CMXB  to  GMXE 
depending  on  which  set  of  sound  pressure  levels  were  being  examined). 
The  nominal  band  center  frequencies  associated  with  the  maximum 


1  3 


*  *•  '  O* 


tone-correction  factors  are  found  from  FMX  »  FREQ( JMX) ,  FMXB  = 
FREQ(JMXB),  FMXC  *  FREQ(JMXC) ,  and  so  on.  See  descriptions  of 
C(I),  FMX.  and  FREQ. 


KAA,  RAP  for  the  propagation  distances  at  the  time  of  maximum  test-day 

A-weighted  sound  level,  ALM,  and  at  the  time  of  maximum  test-day 
tone-corrected  perceived  noise  level,  PNLTM,  and  for  the  corresponding 
heights  of  the  aircraft  reference  point  above  the  microphone  height, 
AHSAM  and  AHSPTM,  index  counters  KAA  and  KAP  represent  the  heights 
from  the  array  HM(K)  of  heights  for  meteorological  data  which  are 
equal  to,  or  just  less  than,  the  corresponding  height  of  the  aircraft 
reference  point  above  the  microphone  height,  i.e.,  KAA  =  K  -  1  when 
[HM(K)  -  2M]  is  first  greater  than  AHSAM  in  a  search  from  K  *  1  to 
K  ■  NH  and  correspondingly  for  KAP  when  -  2M )  is  first  greater 

than  AHSPTM. 


KAA 1 ,  like  the  heighc-arrav  index  counters  KAA  and  KAP  defined  above, 

KAP  I 

except  that  KAAl  ■  KAA  -  1  and  KAP1  -  KAP  -1  to  specify  those 
heights  from  the  HM(K)  array  which  are  stored  in  the  locations  that 
are  1  less  than  the  locations  defined  by  KAA  and  KAP 

KK  index  of  the  height  array  HM(K)  when  HM(K)  is  first  greater  than 

the  microphone  height  ZM  in  a  search  from  1  to  NH.  To  find  the 
lengths  of  the  propagation-path  segments  in  arrays  DAM(K)  and  DPTM(K), 
index  KK  is  used  as  the  starting  index  for  calculations  backwards 
along  the  sound  propagation  path  from  the  microphone  to  the  aircraft 
reference  point. 

KLA,  KLP  final  indexes  of  the  propagation-pathlength-segment  arravs  DAM(K^ 
and  DPTM(K)  for  the  last  segment  of  the  path  from  the  microphone 
back  to  the  aircraft  reference  point:  KLA  ■  KAAl  -  KK  +  3  and 
KLP  »  KAP1  -  KK  +  3.  Indexes  KLA  and  KLP  are  used  in  the  lavered- 
atmosphere  procedures  of  methods  C,  D,  and  E  to  define  the  range  of 
steps  in  the  calculation  of  the  test-day  attenuation  over  the 
propagation  path  in  each  frequency  band. 


f 

I 

j 

I 


I 


I 

J, 

i; 

\ 


LPAGE, 

LPGTOT 


MIC 


MOLAR 


MONTH 


NH 


NPG 


NS 


NUMINT 


page  number  counters  used  to  label  the  pages  of  the  output  listing 
the  test-day  sound  pressure  levels  at  0.5-s  intervals  and  the 
associated  values  of  PNL,  PNLT,  AL,  and  FMX.  LPGTOT  is  the  total 
number  of  pages  of  input  data,  printed  with  a  maximum  of  14 
columns  of  SPL  per  page,  and  LPAGE  is  the  running  page  number 
counter  with  a  value  between  1  and  LPGTOT. 

abbreviation  for  microphone  number,  a  label  assigned  during  the 
flyover  noise  measurement  to  identify  the  microphone  and  its 
location.  MIC  is  read  in  from  the  input  data  file  and  is  only  used 
to  label  the  output  pages. 

name  of  a  subroutine  that  calculates  the  humidity  of  the  atmos¬ 
phere  as  the  molar  concentration  of  water  vapor,  in  percent,  given 
the  air  temperature,  in  degrees  kelvin,  relative  humidity,  in 
percent,  and  air  pressure,  in  standard  atmospheres 

integer  array  read  in  from  the  input  data  file  and  consisting  of 
one  4-character  computer  word  giving  the  month  of  the  year  when  the 
flyover  noise  data  were  measured  (e.g.,  OCT  or  JAN).  MONTH,  DATE, 
and  YEAR  are  used  as  label  information  for  the  output  pages 
containing  the  test-day  sound  pressure  levels. 

integer  constant  read  in  from  the  input  data  file  that  specifies  the 
total  number  of  heights  (in  array  HM)  where  the  test-day  meteorolo¬ 
gical  data  (in  arrays  TC  and  HR)  were  measured 

symbol  for  an  integer  variable  used  in  the  DO  loop  for  printing  out 
pages  2  and  3  of  the  summary  output  pages.  Information  on  page  2 
(the  EPNL  page)  is  printed  when  NPG  ■  1,  that  on  page  3  (the  SEL 
page)  when  NPG  »  2. 

integer  constant  read  from  the  input  data  file  to  specify  the  total 
number  of  0.5-s  samples  of  test-day  1/ 3-octave-band  sound  pressure 
levels  which  are  available  and  are  to  be  read  in  from  the  input  data 
file  as  the  array  SPL(I,J)  where  I  ranges  from  1  to  NS  for  the  time 
index  and  J  ranges  from  1  to  24  for  the  frequency-band  index 

name  of  a  subroutine  that  performs  the  numerical  Integrations  in 
method  D  for  determining  the  1 / 3-octave-band  adjustment  factors  for 
the  differences  in  attenuation  over  the  propagation  pathlength  as  a 
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I 


I 

I 


PA 


PAR 


PDAM, 

PDPTM 


PHIAM, 

PHIPTM 


PI 
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I 

i 

result  of  differences  in  atmospheric  absorption  under  test  and  , 

reference  meteorological  conditions.  NUMINT  returns  an  array  of  the 

band- adjustment  factors,  in  dB,  given:  (1)  an  array  of  24  test-day  | 

sound  pressure  levels  in  dB,  (2)  the  reference-day  uniform-atmosphere 

air  temperature,  TKR,  in  degrees  kelvin,  (3)  the  reference-day 

uniform-atmosphere  relative  humidity,  HRR,  in  percent,  (4)  the 

reference-day  uniform-atmosphere  air  pressure,  PAR,  in  standard 

atmospheres,  (5)  the  array  of  average  test-day  air  temperatures,  ATK, 

in  degrees  kelvin,  (6)  the  array  of  average  test-day  relative 

humidities,  AHR,  in  percent,  (7)  the  test-day  atmospheric  pressure, 

PA,  in  standard  atmospheres,  (8)  the  array  of  propagation-pathlength  i 

segments  along  the  sound  propagation  path  from  the  microphone  back  j- 

to  the  aircraft  reference  point,  and  (9)  the  index  counter  for  the 

last  segment  in  the  array  of  propagation  path  segments.  If  the  j 

array  of  test-day  sound  pressure  levels  is  given  by  SPLPTM(J) ,  with 

the  array  of  propagation  path  segments  given  by  DPTM(K)  and  the  last 

index  counter  specified  by  KLP,  then  NUMINT  returns  the  array  of  * 

band- adjustment  factors  BAPTD(J) .  Similarly,  SPLAM(J)  with  DAM(J) 

and  KLA  returns  the  array  BAAD(J) . 


test-day  atmospheric  pressure,  in  standard  atmospheres,  calculated 
from  BP/RBPO  and  assumed  to  be  uniform  over  the  sound  propagation 
path 

reference-day  air  pressure  for  a  uniform  atmosphere  and  set  equal  to 
1.0  standard  atmosphere 

total  propagation  distance  from  the  aircraft  reference  point  to  the 
microphone  at  the  relative  time  of  occurrence  of  ALM  with  sound 
directivity  angle  PSIAM  for  the  distance  PDAM,  and  at  the  relative 
time  of  occurrence  of  PNLTM  with  sound  directivity  angle  PSIPTM 
for  the  distance  PDPTM,  m 

angle  $  (or  PHI)  defines  the  location  of  the  aircraft  on  the  flight 
path  at  the  relative  time  of  occurrence  at  the  microphone  of  ALM  and 
PNLTM,  respectively.  Angle  ♦,  in  radians,  is  measured,  relative  to 
the  direction  of  flight,  between  the  flight  path  and  the  line  from 
and  aircraft  reference  point  to  the  microphone. 

3.141592654  radians  ( 360°) 


A 


PNL 


array  of  perceived  noise  levels  at  0.5-s  intervals  corresponding  to 
the  test-day  sound  pressure  levels  in  array  SPL,  dB 


PNLM 


PNLMB, 
PNL.MC , 
PNLMD, 
PNLME 


pnl: 


PNLTM 


?  NL  TMB , 
PNLTMC, 
PNL7MD, 
PNL  THE 


PS I AM. 
PS  r?TM 


P  36  TC 


maximum  value  of  test-day  perceived  noise  levels  in  arrav  PNL  (the 
last  one  if  there  is  more  than  one  PNL  with  the  same  maximum 
value) ,  dB 

similar  to  ALMB  to  ALME  above,  except  for  perceived  noise  levels 
instead  of  A-weighted  sound  levels  and  except  that  the  test-day 
spectrum  SPLPTM  at  the  time  of  occurrence  at  the  microphone  of 
PNLTM  is  adjusted  to  reference-day  conditions  by  methods  B,  C, 

D,  and  E  to  yield  the  spectra  SPLPTB  to  SPLPTE  from  which  the 
perceived  noise  levels  PNLMB  to  PNLME  are  calculated,  dB 

array  of  tone-corrected  perceived  noise  levels  corresponding  to  the 
array  of  test-day  sound  pressure  levels  in  array  SPL  and  perceived 
noise  levels  in  array  PNL,  dB 

maximum  value  of  test-day  tone-corrected  perceived  noise  levels  in 
array  PNLT  (the  last  one  if  there  is  more  than  one  PNLT  with  the 
same  maximum  value) ,  dB 

similar  to  PNLMB  to  PNLME  above  except  for  tone-corrected  per¬ 
ceived  noise  levels  in  dB.  The  values  are  found  from  PNLTMB  = 

PNLMB  +  CMXB ,  PNLTME  =  PNLMC  +  CMXC,  etc. 

at  the  same  relative  times  as  for  angles  PRIAM  and  PHIPTM  above, 
but  for  the  sound  directivity  angles  PSIAM  and  PSIPTM  between  the 
flight  path  and  the  line  to  the  microphone  from  the  point  on  the 
flight  path  where  the  aircraft  reference  point  was  when  the  aircraft 
emitted  the  sound  which  resulted  in  the  value  of  ALM  or  PNLTM, 
radians.  Note  that  PSIAM  «  PHIAM  -  sin’1 [AM  sin  (PRIAM)]  and 
PSIPTM  «  PHIPTM  -  sin"1 [AM  sin  (PHIPTM)]. 

name  of  a  subroutine  that  calculates  tone-correction  factors 
according  to  the  April  1978  version  of  Part  36  of  the  Federal 
Aviation  Regulations.  Given  a  spectrumof2i  1/ 3-octave-band  sound 
pressure  levels,  P36TC  returns  the  maximum  value  of  the  tone- 
correction  factors  for  that  spectrum  and  the  corresponding 
frequency-band  index  counter  [e.g.,  C ( I )  and  NMX  are  returned  for 
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I  . 


J 


t- 


RALPHA 


RBPO 


RHM 


RHR 


RHRL 


RHRO 

RPA 


RPAL 


RTK 


RTKL 


each  of  the  SPL(I,J)  spectra  [or  each  ASPL(J)],  CMXB  and  JMXB  for 
SPLPTB(J),  etc.] 

for  each  1/ 3-octave  band,  the  reference-day  sound  absorption 
coefficient  calculated  by  subroutine  ARP866A  for  uniform  reference 
meteorological  conditions  with  air  temperature  given  by  TCR  and 
relative  humidity  given  by  HRR  and  used  in  finding  the  band-adjust¬ 
ment  factors  in  methods  B  and  C,  dB/m 

reference,  sea-level,  standard-day  barometric  pressure  of  29.921 
in.  Hg  and  used  with  BP  to  determine  PA 

reference  height  above  ground  level  fo  10.0  m  used  to  define  the 
starting  height  for  the  calculations  of  reference  vertical  profiles 
of  reference  air  temperature  RTK,  reference  relative  humidity  RHR, 
and  reference  air  pressure  RPA 

array  for  the  reference  vertical  profile  of  relative  humidity 
starting  from  a  value  set  by  RHRO  at  the  reference  height  RHM  and 
decreasing  at  the  rate  specified  by  RHRL,  percent 

reference  lapse  rate  for  the  reference  relative  humidities  in  array 
RHR  with  the  value  of  -0.0065  percentage  points/m  (-6.5  percentage 
points/km) 

reference  relative  humidity,  at  the  reference  height  RHM,  with  the 
value  of  70.0  percent 

array  for  the  reference  vertical  profile  of  air  pressure  starting 
from  a  value  of  1.0  standard  atmosphere  at  the  reference  height 
RHM,  standard  atmospheres 

exponent  for  the  lapse  rate  of  air  pressure  in  standard  atmospheres 
or  a  range  of  heights  from  mean  sea  level  to  11  km  and  having  the 
value  of  -5.25588 


array  for  the  reference  vertical  profile  of  air  temperature 
starting  from  a  value  set  by  RTK0  at  the  reference  height  RHM  and 
decreasing  at  the  rate  specified  by  RTKL,  K 

reference  lapse  rate  for  the  reference  air  temperatures  in  array 
RTK  with  the  value  -0.0065  K/m  (-6.5  K/km) 
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KTKO 


aUN 


SEL 


SELB, 

SELC , 

SELD, 
SELE 


SPLAB, 

SPLAC, 

SPLAD, 

SPLAE, 

SPLAM 


SPLPTB, 

SPLPTC, 

SPLPTD, 

SPLPTE, 

SPLPTM 


SUMDAT 


reference  air  temperature  at  the  reference  height  RHM  with  the  value 
of  298.13  K  (i.e.,  25.0  C) 

integer  constant  that  is  read  from  the  input  data  file  to  specify 
the  run  number  assigned  at  the  time  of  rhe  measurements  to  the 
set  of  input  test-day  sound  pressure  levels  and  used  for  label 
information  on  the  output  pages 

sound  exposure  level,  in  dB,  as  calculated  by  subroutine  IN'TEG  from 
test-day  AL  data  with  the  corresponding  value  of  ALM  as  DATA10 
and  with  specified  values  for  IFA,  ILA,  and  NS 

sound  exposure  levels,  in  dB,  calculated  from  the  SEL  duration 
factor,  DSEL,  and  the  maximum  A-weighted  sound  levels  for  reference- 
day  meteorological  conditions  by  methods  B  to  E:  SELB  =  ALMB  +  DSEL, 
SELC  -  ALMC  +  DSEL,  etc. 

a  2-dimensional  array  SPL(I,J)  representing  the  test-day  sound  pressure 
levels  read  from  the  input  data  file  at  the  I- intervals  of  time 
during  the  flyover  noise  recording  for  the  J  1/ 3-octave  bands,  dB 

SPLAM  represents  the  array  of  24  1/ 3-octave-band  test-day  sound 
pressure  levels  that  are  associated  with  the  maximum  test-day 
A-weighted  sound  level,  ALM;  SPLAB  to  SPLAE  represent  the  corres¬ 
ponding  arrays  of  sound  pressure  levels  adjusted  to  uniform  refer¬ 
ence-day  conditions  by  adjustment  methods  B  to  E,  respectively,  dB 

SPLPTM  represents  the  array  of  24  1/ 3-octave-band  test-day  sound 
pressure  levels  that  are  associated  with  the  maximum  test-day  tone- 
corrected  perceived  noise  level,  PNLTM;  SPLPTB  to  SPLPTE  represent 
the  corresponding  arrays  of  sound  pressure  levels  adjusted  to 
uniform  reference-day  conditions  by  adjustment  methods  B  to  E, 
respectively,  dB 

SUMDAT  is  the  value  of  the  summation  returned  by  subroutine  INTEG 

.  ,  ,  _  _ [DATA( I) ]/10  ,  _  . 

as  the  nondimens ional  power  sum  110  between  the  values 

i 

of  the  test-day  times  when  DATA  is  closest  to  DATA10.  If  DATA  is 
the  array  of  PNLT  values  and  DATA10  *  PNLTM  -  10.0,  then  SUMDAT 
is  used  to  determine  EPNL;  if  DATA  is  the  array  of  AL  values  and 
DATA10  “  ALM  -  10.0,  then  SUMDAT  is  used  to  determine  SEL. 


TALM 


TALPHA 


TC 


TCR 


TC10 


TEXT, 

TEXTA 


TITLE1, 

TITLE2 


a  quantity  used  with  the  output  listing  of  test-day  sound  pressure 
levels  to  indicate  the  time  of  occurrence  of  ALM.  TALM,  in  seconds, 
is  determined  by  the  value  of  the  relative  test  time  from  array 
TT(I)  at  the  index  LAM  [i.e.,  TALM  =  TT(IAM) ] . 

for  band-adjustment  method  B,  an  absorption  coefficient  calculated 
by  subroutine  ARP866A  at  a  frequency  defined  in  the  subroutine  for 
the  air  temperature  and  relative  humidity  TC10  and  HR10  at  a  height 
of  10.0  m,  dB/m 

array  of  test-day  air  temperatures  measured  at  the  heights  defined 
by  array  HM  and  read  from  the  input  data  file,  C 

reference-day  air  temperature  for  a  uniform  atmosphere  and  set 
equal  to  25.0  C 

air  temperature  for  test-day  conditions  interpolated,  or  read 
directly,  from  the  measured  vertical-profile  data  for  a  height  of 
10.0  m,  C 

names  for  the  two  arrays  of  literal  characters  that  are  declared  to 
be  labeled  COMMON  and  which  are  defined  in  BLOCK  DATA  for  use  on 
summary  output  pages  2  and  3.  The  372-computer-word  array  TEXT  is 
a  12x31,  2-dimensional  array  of  12  computer  words  on  31  lines  of 
descriptive  text.  Each  line  has  a  maximum  of  48  literal  characters 
at  4  characters  per  computer  word.  There  are  31  possible  lines 
of  text  to  correspond  with  the  da "a  at  the  24  1/3-octave-band 
center  frequencies  with  an  allowance  for  7  blank  lines  that  are  used 
for  spacing  when  printing  the  output  listings.  TEXTA  is  a  12- 
computer-word,  one-dimensional  array  that  is  used  in  place  of 
TEXT(12,3)  when  printing  the  third  line  of  descriptive  text  on  page 
3  of  summary  output.  Page  3  contains  the  values  of  maximum  A- 
weighted  sound  level  and  sound  exposure  level  while  page  2  contains 
the  perceived  noise  levels  and  effective  perceived  noise  levels. 

arrays  of  20  4-character  computer  words  that  are  read  from  the  input 
data  files  and  used  for  titles  on  the  output  pages.  TITLE1  is 
used  on  the  pages  containing  the  input  test-day  sound  pressure 
levels  and  describes  general  characteristics  for  the  type  of  aircraft 
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« 


TKR 


TOH 


TPNLM, 

TPNLTM 


TRAM, 

TRPTM 


TS 


IT 


YEAR 


tested.  TITLE2  is  used  on  the  test-day  SPL  pages  and  also  on 
pages  2  and  3  of  the  summary  output  pages  to  describe  particular 
characteristics  of  the  aircraft,  engine,  and  test-site  location. 

reference-day  air  temperature  for  a  uniform  atmosphere  and  set 
equal  to  298.15  K 

time,  in  seconds,  on  the  day  of  the  test  when  the  aircraft  was 
directly  over  the  microphone  (i.e.,  the  time  at  overhead).  TOH 
is  read  from  the  input  data  file  as  the  time  after  midnight  on 
the  basis  of  86,400  seconds  per  day  (e.g.,  TOH  =  44184  s  means  12 
hours,  16  minutes,  and  24  seconds  after  midnight  or  1216:24  hours). 

similar  to  TALM  above,  except  here  the  symbols  indicate  the  relative 
times,  in  seconds,  from  array  TT  at  the  time  of  occurrence  of 
PNLM  [TPNLM  =  TT(IPM)]  and  of  PNLTM  [TPNLTM  =  TT(IPTM)]. 

relative  times  of  occurence,  in  seconds,  of  ALM  and  PNLTM.  The 
relative  times  are  determined  relative  to  the  time  of  day  TOH  when 
the  aircraft  was  directly  over  the  microphone  and  the  time  TS  of 
the  first  sample  of  test-day  sound  pressure  level.  A  0.23-s 
factor  is  included  in  the  calculation  of  TRAM  and  TRPTM  to  shift 
the  relative  time  scale  to  the  center  of  the  0.5-s  data  sample; 
thus  TRAM  =  TS  +  0.25  +  TT(IAM)  -  TOH  and  TRPTM  =  TS  +  0.25  + 
TT(IPTM)  -  TOH. 

time  of  day,  in  seconds,  read  from  the  input  data  file  and  defined 
as  the  time  at  the  beginning  of  the  first  0.5-s  sample  of  test-day 
sound  pressure  levels  [i.e.,  at  the  start  of  the  first  relative 
test  time  TT(1)],  and,  like  TOH,  on  a  86,400-s  daily  basis 

array  of  times  associated  with  each  set  of  test-day  sound  pressure 
levels,  seconds.  The  test  times  in  array  TT  are  determined 
relative  to  0.0  seconds  at  the  midpoint  of  the  averaging  period 
for  the  first  data  sample.  For  computational  purposes,  each  data 
sample  was  considered  to  have  been  averaged  for  exactly  0.5  seconds. 

a  number,  declared  to  be  an  integer  constant,  read  from  the  input 
data  file  (e.g.,  as  74  or  76)  and  used  with  DATE  and  MONTH  for 
specifying  when  the  flyover  noise  measurements  were  made 


ZM 


height  of  the  microphone  above  the  ground  plane  for  the  flyover 
noise  measurements  and  read  from  the  input  data  file,  m 
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3.  LISTING  OF  PROGRAM  STATEMENTS 

TESTREF  is  the  name  for  the  main,  or  driver,  program  that  calculates 
spectral  adjustment  factors  for  differences  in  atmospheric  absorption. 

The  program  is  specifically  applicable  to  aircraft  noise  measurements  made 
under  test  conditions  and  adjusted  to  reference  meteorological  conditions. 
Four  separate  adjustment  procedures  are  included. 

The  program  reads  in  1/3-octave-band  sound  pressure  levels  at  regular 
intervals  of  time  during  a  recording  of  aircraft  flyover  noise.  Meteorolo¬ 
gical  and  other  associated  test  data  are  also  read  in.  The  required  calcu¬ 
lations  are  performed  by  nine  subroutines  that  are  called  by  the  main 
program.  An  array  of  labeled  common  is  included  as  BLOCK  DATA  and  contains 
the  text  of  explanatory  legends  for  some  of  the  tabular  output  listings. 

Program  TESTREF  and  the  BLOCK  DATA  require  about  6825  words  or  28,000 
bytes  of  source-code  storage  capacity.  With  subroutines,  approximately 
37,500  bytes  are  required  for  storage  of  the  total  source  code. 

The  total  number  of  lines  of  computer  program  coding  are  substantially 
more  than  required  for  reading  in  the  data,  performing  the  calculations, 
and  printing  the  output  pages.  To  assist  the  user,  several  comments  have 
been  included  in  the  main  program  and  subroutines  to  explain  and  describe 
the  methods  used  ard  options  selected. 

The  program  was  written  in  the  FORTRAN  IV-extended  programming  language. 
It  was  written  for  a  compiler  and  central  processer  that  are  compatible 
with  an  IBM  System  370  computer.  The  program  makes  use  of  certain  features 
permitted  by  the  compiler  that  may  not  be  available  on  another  machine. 

For  the  user's  information,  those  features  that  may  require  special  atten¬ 
tion  during  implementation  of  the  program  are  listed  below. 

1.  Extended  assignment  of  a  value  to  several  variables  simultaneously 
is  permitted,  asA-B=C=D=E=  0.0. 


2.  GOTO  is  a  permitted  format;  GO  TO  is  not  required. 


3.  LOGIO  and  LOG  are  permitted  forms;  AL0G10  and  ALOG  are  not  required 
for  logarithms. 

4.  ASIN  is  permitted  for  the  arcsine;  ARSIN  is  not  required. 

5.  The  name  of  a  function  or  a  variable  can  have  as  many  as  8  charac¬ 

ters;  some  systems  have  a  limit  of  6  characters. 

6.  Arithmetic  operations  are  permitted  in  the  range  of  a  DO  loop 

and  in  the  specification  of  a  READ  or  WRITE  statement,  as  DO  20  J  *  1, 

3*NN  or  READ  (5,  15)  TC  +  273.15. 

7.  Two  or  more  statements  can  be  combined  on  one  line,  but  separated 
by  semicolons  as  PNM  *  -10.0;  PNT  «  0.0. 
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PROGRAM  TESTREF 


1.000 
2.000 
3.00  0 
4.000 
5.000 
6.  000 
7.000 
5.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
15.000 
19.000 
20.000 
21.000 
22.000 
23.000 
24.000 
25.000 
2  o. 000 
27.000 
26.000 
29.000 
30.000 
31.000 
32.000 
33.000 
34.000 
55.000 
56.000 
57.000 
33.000 
39. 000 
40.000 
41.000 
42.000 
43.000 
44.000 
45.000 
46.000 
47.000 
43.000 
49.000 
50.000 
51.000 
52.000 


C:::: 

C 

C  PROGRAM  TESTREF  WAS  PREPARED  BY  DYTEC  ENGINEERING,  INC.  OF 

C  LONG  BEACH,  CA .  VERSION  12  AUGUST  1974. 

C 

C  THE  PURPOSE  OF  THE  PROGRAM  IS  TO  ADJUST  MEASURED  AIRCRAFT 

C  FLYOVER  NOISE  LEVELS  FOR  DIFFERENCES  IN  ATMOSPHERIC  ABSORPTION 

C  LOSSES  UNDER  TEST  AND  REFERENCE  METEOROLOGICAL  CONDITIONS. 

C 

C  ADJUSTMENTS  are  DETERMINED  USING  PURE-TONE  ATMOSPHERIC 

C  ABSORPTION  CALCULATED  BY  THE  PROCEDURES  OF  SAE  ARo  ?ASa  (MARCH 

C  1975)  AND  ANS  SI. 26-1978.  ABSORPTION  LOSSES  AND  ADJUSTMENTS 

C  FROM  test  to  ACOUSTICAL  REFERENCE  METEOROLOGICAL  CONDITIONS 

C  ^  OR  1/3-OCTAVE-BAND  SOUND  PRESSURE  LEVElS  ARE  CALCULATED 

C  BY  A  NUMBER  OF  METHODS. 

C 

C  THE  EFFECT  OF  USING  the  VARIOUS  METHODS  IS  DETERMINED  FOR 

C  1/3-OCTAVE-BAND  SOUND  PRESSURE  LEVELS,  °E  R  C  E I VE  D  NOISE  LEVELS, 

C  TONE-CORRECTED  PERCEIVED  NOISE  LEVELS,  EFFECTIVE  PERCEIVED 

C  NOISE  LEVELS,  MAXIMUM  A-wEIGHT=D  SOUND  LEVELS,  AND  SOUND 

C  EXPOSURE  LEVELS.  VARIOUS  AIRCRAFT  NOISE  SPECTRA,  SOUND 

C  PROPAGATION  DISTANCES,  AND  TEST-DAY  MF tf OR OLOG I C AL  CONDITIONS 

C  ARE  EXAMINED  BY  VARYING  THE  INP'JT  DATA. 

C 


C::;:  INPUT  DATA  FILE: 

c 

C  for  each  SET  of  input  data,  the  INPUT  file  that  is  READ  ON 

C  DEVICE  5  CONTAINS:  RUN;  MIC,  NS,  DATE,  MONTH,  AND  YEAR; 

C  TS,  TOH,  AH,  ZM,  AS,  AM,  AND  GAMAD;  SPL:  NH  AND  HM ;  tq  AND  HR; 

C  BP;  AND  T I T  l  E  1 ,  T I T  L  E  2 ,  AND  GROUND  FOR  LABELS  ON  OUTPUT  DAGE  S . 

C 

C::::  SYMBOLS  AND  UNITS  FOR  PARAMETERS  and  VARIABLES  arc  DEFINED 

C  IN  COMMENTS  AND  IN  ACCOMPANYING  text  In  MORE  DETAIL. 

C 

C :s :s :::::::::: 

C 

DIMENSION  AHR (31 ) , AL C 1 0 0 ) , A S PL ( 2 4 ) , A TC ( 3 1  )  , A T< ( J 1  )  , A W( 2 4  )  , 
1SAAD(24),BAPTD(2  4),C(100),DAM(31),DPTM( 7 1  )  , RMX ( 1 0 n  )  , F R EP ( 2 4  )  , 
2GROUND( 5 ),HM( 31 ),HMEA(31 ),HM I (31 ),HR(31  )  ,HREF( 31 ),MOWTh( 1 ), 

5PNLC 100), PNLT(IOO), RHR(31 ) , R PA ( 3 1 ) , R T< ( 5 1  )  , 5  PL ( 1 0 0 ,  2 u ) , 
4SPLABC24), SPLAC(24), SPLAD(24), SPLAEC24), SPLAM( 25 ), SPLPTBC  24  ), 
5SPLPTC( 24), SPLPTDC24), SPLPTE(24), SPLPTM( 2  5), TC( 31 ), tt( i 00  ), 
6TITLE1(20),TITLE2(20) 

C 

C  THE  DIMENSIONS  OF  AHR  ,  A  TC  ,  A  T<  ,  DAM ,  DPTM  ,  HM  ,  HME  A  ,  HM  I  ,  HR  , 

C  HREF,RHR,RPA,RT<,  AND  TC  ARE  SPECIFIED  FOR  METEOROLOGICAL 
C  MEASUREMENTS  AT  31  HEIGHTS  ABOVE  GROUND.  IF  THERE  ARE  MORE 
C  HEIGHT  INTERVALS  (LAYERS)  THAN  31,  THE  DIMENSION  WILL  HAVE 
C  TO  BE  INCREASED  ACCORDINGLY,  HERE  AND  IN  CORRESPONDING  SUBROUTINES. 
C 
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53.000 

54.000 

55.000 

56.000 

c 

57.000 

c 

58.000 

c 

59.000 

60.000 

61.000 

c 

62.000 

c 

6  3.000 

c 

64.000 

c 

65.000 

66.000 

67.000 

68.000 

c 

69.000 

C:::! 

70.000 

c 

71.000 

c 

72.000 

c 

7  3.000 

c 

74.0  00 

75.000 

10 

76.000 

1 1 

77.000 

1  2 

78.000 

79.000 

13 

80.000 

c 

81.000 

c 

82.000 

c 

83.000 

c 

84.000 

c 

85.000 

86.000 

8  7  .  000 

14 

8  8.000 

C 

89.000 

90.000 

C 

91 .000 

C 

92.000 

C 

93.000 

9  4.000 

15 

9  5.0  00 

C 

96.000 

97.000 

1  7 

98.000 

C 

99.000 

C 

100.000 

C 

101.000 

102.000 

103.000 

C 

104.000 

C 

105.000 

c 

106.000 

107.000 

108.000 

c 

COMMON  /CTEXT/  TEX T ( 1 2 , 3 1 ) , T E XTA (  1  2  ) 

INTEGER  DATE, FREO, RUN, YEAR 
REAL  ISBN 

DATA  FREQ  =  24  NOMINAL  1/ 3-OCTAVE-BAND  CENTER  FREQUENCIES 

DATA  FREQ  / 5 0 , 6 3 , 8 0 , 1 0 0 , 1 2 5 , 1 6 0 , 2 0 0 , 2 5 0 , 3 1 5 , 4 0 0 , 5 0 0 , 6 3 0 , 8 0 0 , 
11000,1250,1600,2000,2500,3150,4000,5000,6300,8000,10000/ 

DATA  AW  =  A-SCALE  FREQUENCY  WEIGHTINGS  AT  THE  24  CENTER 
FREQUENCIES  OF  THE  1/3-OCTAVE  BANDS. 

DATA  AW  / - 3 0 . 2  , -2  6. 2  , -2  2. 5  , -  1  9. 1  ,-1  6. 1  ,-  1  3. 4, -  1  0. 9, -8. 6  , 

1-6. 6, -4. 8, -3. 2, -1.9, -0.8, 0.0, 0.6,  1.0,  1.2,  1.3,  1.2,  1.0, 

20. 5,-0. 1,-1.  1,-2. 5/ 

::::  READ  INPUT  FILE  FOR  AIRCRAFT  NOISE  LEVELS  AND  ASSOCIATED  DATA 

READ  RUN, MIC, NS, DATE, MONTH, YEAR, TS,TOH, AH, ZM, AS, AM,  AND 
GAMAD  . 

READ  (5,10)  RUN 
FORMAT  (13) 

READ  (5,12)  MIC, NS, DATE, MONTH, YEAR 
FORMAT  (I  1,  14,  I  3,  1A4,  13) 

READ  (5,13)  TS,TOH,AH,ZM,AS,AM,GAMAD 
FORMA  T(F  7 . 1,F8.1,F6.1,F4.1,F5.1,F6.3,F5.1) 

CALCULATE  TT ( I ) ,  THE  ARRAY  OF  TEST  Timfs  RELATIVE  TO  The 
TIME  of  the  first  DATA  SAMPLE  ASSUMING  Data  time  INTERVALS 
OF  EXACTLY  0.5  SECONDS. 

TT( 1 )=0 . 0 

DO  14  1=2, NS 

TT( I )  =  T  T ( I -1 )  +  0 . 5 

DO  25  1=1, NS 

READ  SP  L ( I , J ) ,  THE  24  1/ 5-OCTAVE-BAND  SPLS  AT  0.5-S  INTERVALS. 

READ  (5,15)  (SPL(I,J),J=1,24) 

FORMAT (  1  2F6  .  1  ) 

DO  17  J=l, 24 
ASPL(J)=SPL( I , J) 

CALCULATE  PERCEIVED  NOISE  LEVEL,  PNL(I) 

CALL  CPNL(ASPL,APNL) 

PNL ( I ) = APNL 

CALCULATE  TONE-CORRECTED  PERCEIVED  NOISE  LEVELS,  DNL  T (  I  ) 

CALL  P36TC(ASPL,C(I),UMX) 

PNL  T ( I )  =  PNL ( I )+C ( I ) 
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109.000  C  FIND  FMXCI),  THE  NOMINAL  CENTER  FREQUENCY  IN  KHZ  OF  THE 

110.000  C  BAND  THAT  PRODUCED  CMX  AND  HAD  FREQUENCY  BAND  INDEX  COUNTER 
111.000  C  JMX. 

112.000  C 

1  1  3.000  FMXCI )  =  (FLOAT(FREQ(JMX)))/l  000  .  0 

1 14. 000  C 

115.000  C  CALCULATE  A-WEIGHTEO  SOUND  PRESSURE  LEVEL,  AlCI). 

1  16.000  C  SET  AL  (0  =  0.0  DB  IF  ALL  SPL(J)=0.0  DB  . 

117.000  C 

1  1  8.000  AP  =  0  .  0 

119. 000  DO  20  U=1 , 24 

1  20  .  000  20  AP=AP  +  1Q.0::::((SPL([,J)  +  AW(J))/10.0) 

121.000  IFCAP.LT. 15. 0)GOTO  23 

1  22  .  000  AL(1  )=10  .  0::LOG1G(AP) 

123.000  GOTO  25 

124.000  23  AL(I)=0.0 

125.000  25  CONTINUE 

126.000  C 

127.000  C  FIND  PNLM,  THE  MAXIMUM  VALUE  OF  THE  PERCEIVED  NOISE  LEVELS 

1  28.000  C  PNLCO;  IPM,  THE  TIME  INDEX  COUNTER  WHEN  PNL(1)=PNLM;  AND 

1  29.000  C  TPNLM,  THE  RELATIVE  TEST  TIME  TT(IPM)  WHEN  I  =  I  PM . 

150.000  C 

131.000  C  FIND  PNLTM,  THE  MAXIMUM  VALUE  OF  THE  TONE-CORRECTED 

152.000  C  PERCEIVED  NOISE  LEVELS  PNLT(I);  IPTM,  THE  TIME  INDEX  COUNTER 
133.000  C  WHEN  PNLTCI )=PNLTM;  AND  TPNLTM,  THE  RELATIVE  TEST  TIME 
134.000  C  TT(IPTM)  WHEN  IrlPTM. 

1 35. 000  C 

136.000  C  ALSO,  FIND  ALM,  THE  MAXIMUM  VALUE  OF  THE  A-WEIGHTED  SOUND 

1  37  .  000  C  PRESSURE  LEVELS  AL(I);  I  AM ,  THE  TIME  INDEX  COUNTER  WHEN 

158.000  C  AL(I)=ALM;  AND  TALM,  THE  RELATIVE  TEST  TIME  TT(IAM)  WHEN 

1  39.000  C  I = I  AM . 

140 . 000  C 

141.000  PNLM=PNLTM=ALM=-1 000 . 0 

142.000  DO  30  1=1, NS 

143.000  IFCPNLCI ).LE.PNLM)GOTO  26 

144.000  PNLM=PNL(I) 

145.000  I PM= I 

146.000  TPNLM=TT(IPM) 

147.000  C 

148.000  26  IFCPNLTCI ) .LE . PNLTM)GOTO  27 

1  49.000  PNLTMrPNLT ( I  ) 

150.000  I PTM= I 

151.000  TPNLTM=TTCIPTM) 

152. 000  C 

1  53.000  27  IF  C  AL  C I ) . LE . ALM)GOTO  30 

154.000  ALM=AL( I ) 

155.000  I AM= I 

156.000  TALMrTT(IAM) 

157.000  30  CONTINUE 

158. 000  C 

159.000  C  LABEL  THE  1 /3-OCTAVE -BAND  SPLS,  SPLPTM(U),  CORRESPONDING 

160.000  C  TO  THE  TIME  OF  OCCURRENCE  OF  PNLTMrPNLT ( I PTM)  . 

161.000  C 

162.000  DO  33  J=1 , 24 

163.000  33  SPLPTM(U)=SPL(IPTM,J) 

164 . 000  C 


r 

c 

i 

1 

1  ; 

LABEL  THE  1/3-OCTAVE-BAND  SPLS  ,SPLAM(j),  j 

1  66  .  000 

c 

CORRESPONDING  TO  THE  TIME  OF  OCCURRENCE  l 

167. 000 

c 

OF  ALM  =  AL(IAM) .  j 

168. 000 

c 

I 

169.000 

DO  34  J= 1 , 24 

170.000 

34 

splam(j)=spl(iam,u) 

171.000 

C 

1  72  .  000 

C 

READ  NH,  THE  NUMBER  OF  HEIGHTS  ABOVE  GROUND  LEVEL  WHERE 

j  173.000 

c 

METEOROLOGICAL  DATA  WERE  MEASURED,  AND  HM(K),  THE  VALUES 

174.000 

c 

OF  THE  HEIGHTS. 

175.000 

c 

176. 000 

READ  C 5  ,  35)  NH 

1 77 . COO 

35 

FORMAT ( I 3) 

178.000 

READ  (5,36)  (HM(<) ,  K  =  1 ,  NH) 

179. COO 

36 

F  ORMAT (F  6 . 1)  '  ; 

180.000 

C 

>  --1 

181.000 

C 

READ  TCOO  AND  HR(K),  THE  AIR  TEMPERATURE  AND  RELATIVE 

182.000 

c 

HUMIDITY  AT  THE  TIME  OF  THE  FLYOVER  NOISE  MEASUREMENT  J 

1  83.  000 

c 

AND  AT  THE  VARIOUS  HEIGHTS  ABOVE  GROUND  LEVEL.  ^ 

184.000 

c 

,  i 

185.000 

c 

READ  BP,  THE  BAROMETRIC  PRESSURE  AT  THE  TIME  OF  THE 

186. 000 

c 

FLYOVER  NOISE  TEST.  THE  BAROMETRIC  PRESSURE  IS  USUALLY  i 

187.000 

c 

MEASURED  AT  A  HEIGHT  OF  10  M. 

188.000 

c 

189 . 000 

READ  (5,37)  ((TC(K),HR«)),K  =  1,NH)  f  j 

190. 000 

37 

FORMAT ( 2F  5 .  1 ) 

191  -  000 

READ  (5,38)  BP 

192 . 000 

38 

FORMAT (F6 . 3 )  ; 

1 93. 000 

C 

194.000 

c 

READ  INFORMATION  FOR  LABELING  OUTPUT  PAGES. 

195 . 000 

c 

196. 000 

READ  (5,39)  TITLE1  j 

197 . 000 

39 

F ORMAT (2  0A4 ) 

198. 000 

READ  (5,40)  TITLE2 

199. 000 

40 

FORMAT (2  0A4) 

200.000 

READ  (5,41)  GROUND 

201.000 

41 

FORMAT (5A4) 

202 . 000 

C 

203. 000 

C:: :: 

—  END  OF  READING  OF  INPUT  DATA  AND  OF  INITIAL  CALCULATIONS. 

204. 000 

c 

205 . 000 

c 

CALCULATE  THE  EFFECTIVE  PERCEIVED  NOISE  LEVEL,  EPNL ,  FROM 

206. 000 

c 

THE  INPUT  ARRAY  OF  PNLT(I)  VALUES.  ALSO,  FIND  INDICES  IFP  AND 

207.000 

c 

ILP  FOR  THE  RELATIVE  TIMES  AT  THE  START  AND  END  OF  THE  SUMMATION  u 

208 . 000 

c 

OVER  THE  1 0 -D9 -DOWN  TIME  OF  THE  SQUARE  OF  THE  TONE-CORRECTED  J, 

209. 000 

c 

PERCEIVED  NOISINESS.  | 

210.000 

c 

1  ; 

211 . 000 

c 

INTEGRATION  FOLLOWS  THE  RULE  OF  SECTION  636.9(E)  AND  IS 

212.000 

c 

JUST  BETWEEN  THE  RELATIVE  TIMES  OF  THE  0.5-SEC  DATA  SAMPLES  1 

213. 000 

c 

THAT  ARE  CLOSEST  TO  THE  VALUE  OF  PNLTM-10.0  DB .  INTEGRATION  i 

214.000 

c 

IS  PERFORMED  BY  SUBROUTINE  INTEG.  ■ 

215.000 

c 

,  „• 

216.000 

DATA10=PNLTM-1 0 . 0  - 

217.000 

CALL  INTEG(NS,PNLT,DATA10, IFP, ILP, SUMDAT) 

218. 000 

EPNL  =  1 0 . 0 "LOG 10(SUMDAT)-10. 0  "LOG  1  0  ( 2  . 0)  -10.0  : 

219.000 

c 

28 

4 

1 

4 

L 

.  j 

--'I'--  •*  iw  •-  v •  ^ •  "■ 

220.000  C 
221.000  C 
222. 000 
223.000  C 
224. 000  C 
225.000  C 
226.000  C 
227.000  C 
228.000  C 
229. 000  C 
230.000 
231.000 
252.000 
233.000  C 
234.000  C 
235.000  C 
236. 000 
237.000  C 
238.000  C 
2  39  .  000  C 
240.000  C 
241.000 
242.000 
243 . 000  C 
244.000  C 
245.000  C 
246. 000  C 
247. 000  C 
248. 000 
249.000 
2  50  .  0  00 
251.000 
252.000  C 
253.000  C 
254.000  C 
255.000  C 
256.000  C 
257.000  C 
258.000  C 
259.000 
2  60  .  0  00 
261.000  C 
262.000  C 
263.000  C 
264. 000  C 
265.000 
2  66 . 000  C 
267. 000  C 
268.000  C 
269. 000  C 
270.000  C 
271.000  C 
2  72 . 000  C 
273.000  C 
274.000  C 
275.000  C 
2  76  .  000  C 


FIND  duration  CORRECTION  FACTOR,  DCF 
DCFrEPNL-PNLTM 

CALCULATE  the  SOUND  EXPOSURE  LEVEL,  SEL,  from  the  array  of 
INPUT  A-WE1GHTED  SOUND  LEVELS,  AL ( ! ) .  ALSO  FIND  the  INDEXES 
I F  A  AND  I  LA  FOR  THE  RELATIVE  TIMES  AT  the  START  AND  END  OF 
SUMMATION  OVER  THE  10-DB-DOWN  TIME  FOR  THE  SQUARED  A-WEIGHTED 
SOUND  PRESSURE. 

DATA 10=ALM- 10.0 

CALL  INTEGCNS, AL, DATA10, IFA, ILA,SUMDAT) 

S  EL  =  1  0  .  0"LOG10(SUMDAT)-1  0  .  0  :;LOG  1  0  (  2  .  0  ) 

FIND  DSEL,  DURATION  OF  SEL 

DSEL=SEL-ALM 

DETERMINE  TEST-DAY  AIR  PRESSURE  IN  STANDARD  ATMOSPHERES 
USING  A  REFERENCE  ATMOSPHERIC  PRESSURE  OF  29.921  IN.  OF  HG . 

R8PO=29.921 

PA=BP/R8PO 

FOR  A  UNIFORM  REFERENCE  ATMOSPHERE,  DEFINE  REFERENCE  AIR 
TEMPERATURE  IN  DEGREES  CELSIUS  AND  KELVIN,  RELATIVE  HUMIDITY 
IN  PERCENT,  AND  AIR  PRESSURE  IN  STANDARD  ATMOSPHERES. 


TCR  =  2  5 . 0 
TKR  =  2  98  .  1  S 
H  R  R  =  7  0  .  0 
PAR=1  .  0 

IF  APPLICABLE,  USE  LINEAR  INTERPOLATION  ON  the  VALUES  OF 
AIR  TEMPERATURE  AND  RELATIVE  HUMIDITY  AT  The  FIRST  AND  SECOND 
HEIGHTS  REAO  FROM  THE  INPUT  FILE  TO  OBTAIN  CORRESPONDING  VALUES 
AT  A  HEIGHT  OF  10.0  M.  (HM(2)  IS  ASSUMED  T0  BE  >  10.0  M  AND 
HM(1)  IS  ASSUMED  TO  BE  <  OR  =  10.0  M.) 

TC  10  =  TCC  1  )-((TC(  1  )-TC(2))::((  1  0  .  0-HMC  1  ))  /CHMC2)-HM(  1  )))) 

HR  10=HR(1  )-C(HR(  1  )-HR(2))::CC  10  .  0-HMC  1  ))  /  CHMC  2)-HMC  1  )))) 

FIND  THE  MOLAR  CONCENTRATION  OF  WATER  VAPOR  FOR  THE 
TEST-DAY  CONDITIONS  AT  10  M. 

CALL  MOLAR(tc  10+2  73  .  1  5  ,  HR  1  0 ,  PA ,  H 1  0  ) 


DEFINE  THE  REFERENCE  HEIGHT,  RHM ,  OF  10.0  M  AS  THE  HEIGHT 
AT  WHICH  THE  REFERENCE  AIR  TEMPERATURE,  RT<0  OF  298.15  DEG  K, 
REFERENCE  RELATIVE  HUMIDITY,  RHRO  OF  70.0  PERCENT,  AMD 
REFERENCE  AIR  PRESSURE  OF  1.0  STANDARD  ATMOSPHERE  ARE  SPECIFIED. 

ALSO  DEFINE  THE  REFERENCE  LAPSE  RATES  AS  A  FUNCTION  OF  HEIGHT 
ABOVE  THE  REFERENCE  HEIGHT  OF  RT<l=  -0.00FF  DEG  K/M, 

RHRL=  -0.0065  PERCENTAGE  POINTS/M,  AND  RPAl=  -5.25FRR  FOR  the 
EXPONENT  OF  THE  EXPRESSION  FOR  LAPSE  RATE  OF  The  REFERENCE 
AIR  PRESSURE  IN  STANDARD  ATMOSPHERES. 
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2  7  7  .  0  C  0 
278.000 
279.000 
280.000 
281.000 
232.000 
23J.000 
234.000 
285.000 
296.00) 
287.00] 
288.000 
289.000 
290.000 
291.000 
292.000 
293.000 
294.000 
295.000 
2  96 . 0  00 
297.000 
298.000 
299.700 


C 

R  HM  =10.0 
R  1X0=  29  8  .  1  5 
RiXL  =  -0  .  0  06  5 
RHRO=  70  .  0 
R  HR  L  =  -  0 . 006  5 
RPAL=-5 . 25588 
C 

C  NOW  CALCULATE  the  VALUE  op  the  HEIGHTS  at  The  middle 

C  OF  each  height  interval  (layer). 

c 

C  cOR  the  REFERENCE  lapse  RATES,  CALCULATE  the  CORRESPONDING 

C  REFERENCE  TEMPERATURES,  RELATIVE  HUMIDITIES,  an  D  PRESSURES. 

C  ALSO  CALCULATE  The  AVERAGE  VALUES  OF  the  MEASURED  air 

C  TEMPERATURES  AND  RELATIVE  HUMIDITIES  AT  the  MIDPOINTS  OF 
C  THE  HEIGHT  INTERVALS. 

C 

C  MEASURED  TEMPERATURES  are  CONVERTED  from  DEGREES  CELSIUS 

C  TO  KELVINS. 

C 

C  ALSO  CALCULATE  THE  VERTICAL  PROFILES  OF  HM  E  A  AND  hREF  FOR  TME 

C  MOLAR  CONCENTRATION  OF  WATPR  VAPOR  UNDER  TEST-DAY  AND  REFERENCE- 
C  DA  s'  lAPSE-RATe  CONDITIONS. 


3  00.  ]  00 

301.000 

302.000 

303.000 

304.000 

505.000 

306.000 

307.000 

508.000 

509.000 

310.000 

311.000 

312.000 

313.000 

314.000 

315.000 

316.000 

317.000 

318.000 

319.000 

320.000 

321.000 

322.000 

323.000 

324.000 

325.000 

326.000 

327.000 

323.000 

329.000 

330.000 

331.000 

332.000 

333.000 

334.000 


C 

DO  42  K  =  1 , NH- 1 

HMI (k)=(hm(k+1 )+HM(K))/2 . 

RTK(K)=RTKO  +  RT<L-”(HMI  (k)-Rhm) 

RHR(K)  =  RHRO  +  RHRL::(HMI(K)-RHM) 

RPA(K)  =  ((RTKO-10.0)/((RTKO-10.0)  +  (RTKL::HM1(K))))::::RPAL 
ATC(K)=(TC(K+1)+tc(K))/2. 

At<(K)=ATC(K)  +  273.  1  5 

AHR(K)=(HR(<+1)+HR(K))/2. 

call  MOLARCATK(K), AHR(K), PA,HMEA(K)) 

42  CALL  MOLAR(RTK(K),RHR(K),RPA(K),HREF(K)) 

c 

C  IF  HM ( 1 ) ,  THE  FIRST  height  AT  WHICH  METEOROLOGICAL  DATA 

C  WERE  MEASURED,  IS  >  ZM,  THEN  THE  METEOROLOGICAL  PARAMETERS 
C  CALCULATED  ABOVE  WILL  NOT  BE  DEFINED  OVER  THE  LAYER  OF  THE 
C  ATMOSPHERE  THAT  INCLUDES  THE  MICROPHONE  AT  HEIGHT  ZM , 

C  IN  THIS  CASE,  RENUMBER  THE  INDEX  OF  THE  PARAMETERS 

C  (INCREASING  IT  BY  1)  AND  DEFINE  A  NEW  SET  OF  PARAMETERS  FOR 
C  INDEX  1  EQUAL  TO  THE  PREVIOUS  SET  FOR  OLD  INDEX  1  (OR  NEW 
C  INDEX  2),  THEREBY  MAKING  THE  PARAMETERS  IDENTICAL  FOR  ThE  FIRST 
C  TWO  LAYERS. 

C 

IF (HM( 1 ) . GT . ZM)GOTO  43 
GOTO  47 

43  DO  45  L  =  1 , NH 
LL  =  (NH  +  2  )-L 

HMI ( L  L ) =HM I ( LL - 1 ) 

RTK(LL)  =  RTK(LL-1  ) 

RHR(LL)=RHR(LL-1 ) 

RPA(LL)=RPA(LL-1 ) 

ATC(LL)=ATC(LL-1 ) 

AT<(LL)=ATK(LL-1 ) 
ahr(ll)=ahr(ll-i  ) 

HMEA(LL)=HMEA(LL-1 ) 

45  hREF(LL)=HREF(LL-1 ) 
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335 . 000 
336.000 
337 . 000 
338.000 
339.000 
3  40 . 0  JO 
341.000 
342.000 
343 . 000 
344 . 000 
345 . 000 
346.000 
347 . 000 
348.000 
349.000 
350. 000 
351.000 
352.000 
353. 000 
354. 000 
355.000 
356.000 
357.000 
358.000 
359  .  00  0 
360.000 
361.000 
362.000 
363.000 
364.  000 
36  5  .  000 
566  .  00  0 
367. 000 
368. 000 
369.000 
370.000 
371 . 000 
372.000 
373.000 
374. 000 
375-000 
376.000 
377.000 
378.000 
379-000 
380 . 000 
381.000 
382.000 
383 . 000 
384 . 000 
385 . 000 
386 . 000 
387.000 


C 

C  FOR  THE  RELATIVE  TEST  TIMES  TT(IPTM)  AND  TT(IAM) 

C  CORRESPONDING  TO  PNLTM  AND  ALM,  DETERMINE  THE  TOTAL 

C  LENGTHS  OF  THE  SOUND  PROPAGATION  PATHS,  PDPTM  AND  PDAM, 

C  AND  THE  LENGTHS  OF  THE  VARIOUS  SEGMENTS  DPTM(K)  AND 

C  DAMCO,  ALONG  THE  PATHS  THAT  CORRESPOND  TO  THE  INTERVALS 

C  OF  HEIGHT  FOR  THE  METEOROLOGICAL  MEASUREMENTS. 

C 

C  DETERMINE  THE  TIMES,  TRPTM  AND  TRAM,  FOR  THE  POSITIONS  ON 

C  THE  FLIGHT  PATH  WHERE  THE  AIRPLANE  WAS  WHEN  THE  SOUND  PRESSURE 
C  LEVELS  SPLPTM(J)  AND  SPLAM(j),  CORRESPONDING  TO  PNLTM  AND  ALM, 

C  WERE  RECEIVED  AT  THE  MICROPHONE. 

C 

C  TIMES  TRPTM  AND  TRAM  ARE  CALCULATED  RELATIVE  TO  THE 

C  TIME,  TOH,  WHEN  THE  AIRPLANE  WAS  DIRECTLY  OVER  THE  MICROPHONE 
C  AND  TO  THE  DATA  ANALYSIS  START  TIME,  TS. 

C 

C  A  FACTOR  OF  0.25  SECONDS  IS  INCLUDED  TO  SHIFT  the 

C  RELATIVE  TIME  SCALE  TO  THE  MIDDLE  OF  EACH  0.5-S 

C  SPL  DATA  SAMPLE. 

V.. 

47  TRPTMrT  S  +  0 .  25+TT ( 1 PTM)-TOH 
TRAM=TS+0. 25+TTC I AM3-TOH 

C 

C  find  the  DISTANCE  BETWEEN  THE  AIRCRAFT  AND  THE  MICROPHONE 

C  AT  THE  TIME  THE  AIRCRAFT  WAS  OVERHEAD. 

C 

AMHrAH-ZM 

C 

C  DEFINE  PI 

c 

PI  =3 .  14  1  592654 

C 

C  SPECIFY  THE  FLIGHT  PATH  ANGLE,  GAMAR,  IN  RAD  I ANS--POS I T 1 VE 

C  FOR  CLIMBING. 

C 

GAMAR zGAMAD”P 1/180.0 
C 

C  CALCULATE  DM,  THE  MINIMUM  DISTANCE  BETWEEN  THE  MICROPHONE 

C  AND  THE  FLIGHT  PATH;  DB ,  THE  DISTANCE  ALONG  THE  FLIGHT  PAT- 
C  FROM  THE  OVERHEAD  POINT  TO  THE  POINT  OF  CLOSEST  APPROACH;  AND 
C  DISTANCES  ALONG  THE  FLIGHT  PATH  FROM  THE  OVERHEAD  POINT  TO 
C  THE  AIRCRAFT  POSITION  AT  THE  TIME  THE  SOUND  WAS  RECEIVED  AT 
C  THE  MICROPHONE,  I.E.,  DISTANCES  DOPTM  AND  DOAM  AT  RELATIVE 
C  TIMES  TRPTM  AND  TRAM,  RE  SPEC T I VEL Y . 

DM  =AMH::COS  (GAMAR  ) 

DB  =AMH:5S  I  N  (GAMAR  ) 

C 

C  DETERMINE  SOUND  PROPAGATION  ANGLES,  PS  I ,  RELATIVE  TO  THE 

C  FLIGHT  PATH  AND  THE  DIRECTION  OF  FLIGHT,  AT  THE  TIME  WHEN  THE 
C  NOISE  WAS  EMITTED  WHICH  WAS  LATER  RECEIVED  AT  THE  MICROPHONE  AT 
C  RELATIVE  TIMES  OF  PNLTM  AND  ALM. 

C 
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388.000  DOPTM=AS::TRPTM 

3  89.0  00  IFCD8  +  OOPTM.LT.O . 0 ) PH  I PTM=( ATAN( -DM/ ( DB  +  DOPTM) ) ) 

390.000  IF(O8+DOPTM.EQ.0.0)PHIPTM=PI/2.0 

391.000  I F (DB+DOPTM. GT. 0 . 0 )PH I PTM=P I -(ATAN(DM/ (DB+DOPTM) ) ) 

3  92 . 0  U0  PS  1  PTMrPH  I  PTM-AS  I  N(  AM:JS  IN(PHIPTM)) 

593.000  C 

394.000  DOAM=AS"TRAM 

3  95.000  I F ( DB+DOAM . L  T . 0 . 0 )PH I AM  =  ( ATAN( -DM/ (DB+DOAM) ) ) 

396. 0  U0  I F (DB  +  DOAM . EQ . 0 . 0 )PH I AM= P  I  / 2  .  0 

397.000  I F (DB+DOAM. GT . 0 .  0  )PH I  AM=P l -(AT AN (DM/ (DB+DOAM) ) ) 

3  98.0  0  0  PS  I AM=  PH  I  AM- AS  I  N(  AM::  SIN(PHIAM)) 

399. 0u0  C 

4  u  0 . 0  0  0  C  NOW  FIND  THE  TOTAL  LENGTH  OF  THE  PROPAGATION  DISTANCE 

H01. 000  C  FROM  THE  SOURCE  TO  THE  RECEIVER  AT  THE  TIMES  OF 

4U2.UOO  C  PNLTM  AND  ALM. 

403.000  C 

404.000  PDPTM=DM/SIN(PSIPTM) 

+  0  5.0  00  PDAMrDM/ S I N( PS  I  AM) 

406.000  C 

+  0  7 . 0  U  0  C  DETERMINE,  AHS ,  THE  HEIGHT  OF  THE  AIRPLANE  ABOVE  THE 

+08.000  C  MICROPHONE  HEIGHT  AT  THE  TIME  THE  MAXIMUM  NOISE  WAS  RECEIVED 
409.000  C  AT  THE  MICROPHONE. 

410 . OuO  C 

4  1  1  .  0  00  AHSPTM=PDPTM”S 1N(PS I PTM-GAMAR) 

4  1  2.000  AHS  AM  =  PO  AM::  S  I  N(  PS  l  AM-GAMAR  ) 

+13.000  C 

410. 000  C  NOW  DETERMINE  THE  INDICES  FROM  THE  ARRAY  OF  HEIGHTS  FOR 

415.000  C  TmE  METEOROLOGICAL  MEASUREMENTS  THAT  ARE  ONE  AND  TWO  LESS  THAN 
416.000  C  The  HEIGHT  OF  THE  AIRCRAFT  A80VE  THE  MICROPHONE  HEIGHT  AT  THE 
417.000  C  TIMES  OF  MAXIMUM  NOISE. 

418.000  C 

419.000  DO  50  K=1,NH 

420. OuO  IF(hM(K)-ZM.GT. AHSPTM9GOTO  53 

+21.000  50  CONTINUE 

422.000  53  KAP=K-1 

4  2  3.0  00  KAPUKAP-1 

424.000  C 

+25 . U00  DO  55  K=1 ,NH 

426.000  IF(HM(<)-ZM.GT. AHS AM) GOTO  57 

427.000  55  CONTINUE 

426.000  57  <AA=K-1 

429.000  KAA1=KAA-1 

430.000  C 

431.000  C  NOW  FINO  THE  LENGTHS  OF  THE  SEGMENTS  OF  DISTANCE  ALONG 

432.000  C  THE  PROPAGATION  PATHS  CORRESPONDING  TO  THE  INTERVALS  OF 

+33-000  C  HEIGHT  FOR  THE  METEOROLOGICAL  MEASUREMENTS,  DPTM  AND  DAM. 

434.000  C 

435 . 000  DO  60  K=1 ,NH 

436.000  IF((HM(K)-ZM) .GT. 0 . 0)GOTO  65 

4  3  7.0  0  0  6  0  CONTINUE 

438.000  65  KK=K 

439.000  C 

4  +  0  .  O00  DP TM(1)  =  (HM(KK)-ZM) /SI N( PS  I PTM-GAMAR ) 

441.000  DO  70  <=<<,<AP1 

44  2.0  00  7  0  DPTM(K.-<K  +  2)  =  (HM(K+l)-hM(<))/SIN(PS  IRTM-GAMAR) 

443.000  DPTM( KAPl-<<+ 3 )=( AHSPTM-( HM(XAP )- ZM)) /S IN(PS I PTM-GAMAR) 

444.000  C 
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445.000  DAM(1 )=(HM(KK)-ZM)/ SIN  CPS IAM-GAMAR) 

446.000  DO  72  K=KK,KAA1 

447.000  72  DAMCK-KK+2 )=(HM(K+1 )-HM( K) ) / S I N(PS IAM-GAMAR) 

448.00  0  DAMCKAA1-KK+3 )  = C AHSAM- CHMCKAA )-ZM) ) / S I N (PS  I AM-GAMAR ) 

449.000  C 

450.000  C  NOW  COMPUTE  THE  1 / 3 -OC TAVE -B AND  ADJUSTMENT  FACTORS  FOR 

451.000  C  DIFFERENCES  IN  ATMOSPHERIC  ABSORPTION  UNDER  TEST  AND  REFERENCE 
452.000  C  METEOROLOGICAL  CONDITIONS  AND  APPLY  THE  ADJUSTMENTS  TO  THE  SPL 
453.000  C  SPECTRUM  AT  THE  TIME  OF  PNLTM  AND  ALM. 

454.000  C 

455.000  C  NOTE  THAT  IF  THE  TEST-DAY  BAND  LEVEL  IS  =  0.0  DB ,  THEN 

456.000  C  THE  REFERENCE-DAY  BAND  LEVEL  IS  ALSO  =  0.0  DB . 

457.000  C 

458.000  APMB=APMC=APMD=APME=0 . 0 

459.000  C 

460.000  DO  85  J=1 , 24 

461 . 000  C 

462.000  C  FIRST,  DETERMINE  ATMOSPHERIC  ABSORPTION  ADJUSTMENT  FACTORS 

463.000  C  USING  SAE  ARP  866A,  MARCH  1975,  AND  TEST-DAY  AIR  TEMPERATURE 
464.000  C  AND  RELATIVE  HUMIDITY  DATA  AT  A  HEIGHT  OF  10  M. 

465.000  C 

466.000  C  SOUND  ABSORPTION  COEFFICIENTS  FOR  REFERENCE-DAY  METEOROLOGICAL 

467.000  C  CONDITIONS  ARE  CALCULATED  DIRECTLY  FOR  THE  SPECIFIED  CONDITIONS. 

468 . 000  C 

469.000  CALL  ARP86 6A( J , TC 1 0 , HR  1 0 , TALPHA ) 

470.000  CALL  ARP866A(J, TCR, HRR, R ALPHA) 

471.000  IF(SPLPTM(J) . EQ. 0. 0 )GOTO  74 

472.000  SPLPTBC J)=S PLPTM(J)+ CTALPHA-R ALPHA )”PDPTM 

473.000  GOTO  75 

474.000  74  SPLPTB(J)  =  0  .  0 

475.000  75  I F C S PLAMC J).EQ.O.O)GOTO  76 

476. 000  SPLABC  J  )  =  S PLAMC  J  )+ CTALPHA-R ALPHA  )-”PDAM 

477.000  GOTO  77 

478.000  76  SPLABC  J  )  =  0  .  0 

479.000  77  APMB=APMB  +  10.0:s:;CCSPLABCJ)+AWCJ))/10.0) 

480 . 000  C 

481.000  C  SECOND,  DETERMINE  ATMOSPHERIC  ABSORPTION  ADJUSTMENT  FACTORS 

482.000  C  USING  SAE  ARP  866A  FOR  ATMOSPHERIC  ABSORPTION  COEFFICIENTS 
483.000  C  IN  DB/M  AND  AVERAGE  AIR  TEMPERATURE  AND  RELATIVE  HUMIDITY  ALONG 
484.000  C  SEGMENTS  OF  THE  SOUND  PROPAGATION  PATH  DEFINED  BY  HORIZONTAL 
485.000  C  LAYERS  OF  THE  ATMOSPHERE. 

486.000  C 

487.000  C  AN  AVERAGE  SOUND  ABSORPTION  COEFFICIENT  FOR  TEST  METEOROLOGICAL 

488. 0C0  C  CONDITIONS  IS  COMPUTED  FROM  THE  RATIO  OF  THE  TOTAL  ATTENUATION  OVER 
489-000  C  THE  PROPAGATION  PATH  TO  THE  LENGTH  OF  THE  PROPAGATION  PATH. 

490. 000  C 

491.000  C  DEFINE  AN  INDEX  COUNTER  FOR  THE  LAST  SEGMENT  OF 

492.000  C  DISTANCES  GOING  BACK  ALONG  THE  SOUND  PROPAGATION  PATH 
493.000  C  FROM  THE  MICROPHONE  TO  THE  AIRPLANE. 

494 . 000  C 

495.000  KLP=KAP 1 -KK  +  3 

496.000  KLA  =KA A  1 -KK  +  3 

497.000  C 

498.000  ATTPTC=ATTAC=0 . 0 

499 . 000  C 
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500.000  00  80  K=1,KLP 

501.000  CALL  ARP866A( J,ATC(K),AHR(K), TALPHA) 

5  02  .  000  80  ATTPTC=ATTPTC  +  (TALPHA”DPTM(lO) 

503.000  C 

504.000  DO  81  K= 1 , KLA 

5  05.000  CALL  ARP 86 6A( J , ATCC K ), AHR (K ), TALPHA ) 

5  06.000  8  1  ATTAC=ATTAC  +  CTALPHA::DAM(K )) 

507. 000  C 

508.000  AVALFPTC=ATTPTC/PDPTM 

509.000  AVALFAC=ATTAC/POAM 

510.000  IF(SPLPTM(J).EQ. 0. 0)GOTO  82 

5  1  1.000  SPLPTCC  J)  =  SPLP  TM(J)+(AVALF  PTC -R  ALPHA  )::PDPTM 

512.000  GOTO  83 

513.000  82  SPLPTCC J)  =  0 . 0 

514.000  83  IF(SPLAM(J).EQ.O.O)GOTO  84 

5  1  5  .  000  SPLAC  (  J  )  =  S  PLAMC  J)+CAVALFAC  -R  ALPHA  )::P  DAM 

516.000  GOTO  85 

5  1  7  .  000  84  SPLACC J)  =  0  .  0 

5  1  8  .  000  85  APMC=APMC+1  0  .  0 ::::  C  C  S  PLACC  J)+AWC  J))/10  .  0) 

519. 000  C 

520  .  000  C  THIRD,  DETERMINE  ATMOS PHER I C -AB SOR PT I  ON  BAND-ADJUSTMENT 

521.000  C  FACTORS  USING  ANS  SI. 26-1978  FOR  PURE-TONE  ATMOSPHERIC 
522.000  C  ABSORPTION  COEFFICIENTS  AND  COMPUTE  THE  ADJUSTMENT  FACTORS 
523.000  C  BY  INTEGRATING  OVER  THE  BANDWIDTH  OF  THE  1 / 3 -OCTAVE -BAND 
524.000  C  FILTERS  ASSUMING  EACH  FILTER'S  TRANSMISSION  RESPONSE  IS  THAT 
525.000  C  OF  AN  IDEAL  FILTER  WITH  THE  SAME  GEOMETRIC  MEAN  FREQUENCY. 
526. 000  C 

527.000  C  COMPUTE  THE  TOTAL  DIFFERENCE  IN  ATTENUATION  UNDER  TEST 

528.000  C  AND  REFERENCE  METEOROLOGICAL  CONDITIONS  AS  THE  SUM  OF  THE 
529.000  C  DIFFERENCES  IN  ATTENUATION  OVER  SEGMENTS  OF  THE  SOUND 
530.000  C  PROPAGATION  PATH.  THE  SEGMENTS  ARE  DEFINED  BY  HORIZONTAL 
531.000  C  LAYERS  OF  THE  ATMOSPHERE  OVER  WHICH  THE  AIR  TEMPERATURE 
532.000  C  AND  RELATIVE  HUMIDITY  ARE  ASSUMED  TO  HAVE  THE  AVERAGE 
535.000  C  VALUES  GIVEN  AT  THE  HEIGHTS  OF  THE  MIDPOINTS  OF  EACH  LAYER. 
534. 000  C 

535.000  C  CALL  SUBROUTINE  NUMINT  TO  GET  THE  BAND  ADJUSTMENT  FACTORS, 

536.000  C  BA,  FROM  TEST  TO  REFERENCE  METEOROLOGICAL  CONDITIONS  ALONG 
537.000  C  THE  SOUND  PROPAGATION  PATH  FOR  EACH  1/3-OCTAVE  BAND. 

538 . 000  C 

539.000  CALL  NUM INTCSPLPTM, TKR , HRR, PAR, ATK, AHR, PA, DPTM, KLP, BAPTD ) 

540 . 000  C 

541.000  CALL  NUM I NT CSPLAM, TKR , HRR , PAR , ATK , AHR , PA , DAM, KLA , B AAD ) 

542 . 000  C 

543.000  C  CALCULATE  SPLS  ADJUSTED  TO  REFERENCE  CONDITIONS 

544. 000  C 

545.000  DO  93  J=1 , 24 

546. 000  C 

547.000  IFCSPLPTMCJ).EQ. 0. 0)GOTO  90 

548.000  SPLPTDCJ9=SPLPTM( J)+BAPTDCJ) 

549.000  GOTO  91 

5  50.000  90  SPLPTDCJ )  =  0  .  0 

551.000  91  IF(SPLAM(J).EQ. 0. 0 )GOTO  92 

552.000  SPLAD(J)=S PLAMC J)+BAADCJ) 

553.000  GOTO  93 

5  54.000  92  SPLADC J  )  =  0  .  0 

5  55  .  000  93  APMDrAPMD  +  1  0.  0  ::::C  CSPLAD  C  J  )+A  WC  J  ))  / 1  0  .  0) 

556. 000  C 
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557.000  C  FOURTH,  DETERMINE  ATMOSPHERIC-ABSORPTION  BAND-ADJUSTMENT 

5  58  .  000  C  FACTORS,  BA PT E  AND  BAAE ,  USING  ANS  SI. 26-1  97  8  FOR  PURE-TONE 
550.000  C  ATMOSPHERIC  ABSORPTION  COEFFICIENTS. 

560.000  C 

561.000  C  APPROXIMATE  the  ATTENUATION  OVER  the  FREQUENCY  RANGE  OF 

5  62  .  0  00  C  EACH  1 / 3 -OC TA VE -0  AND  FILTER  BY  THE  ATTENUATION  AT  the 
563.000  C  GEOMETRIC  MEAN  FREOUENCY  ONLY.  COMPUTE  the  APPROXIMATE 
564.000  C  DIFFERENCE  IN  TOTAL  ATTENUATION  UNDER  TEST  AND  REFERENCE 
565.000  C  METEOROLOGICAL  CONDITIONS  AS  THE  SUM  OF  THE  DIFFERENCES  IN 
566.000  C  ATTENUATION  OF  THE  GEOMETRIC  MEAN  FREQUENCY  OVER  SEGMENTS 
567.000  C  OF  THE  SOUND  PROPAGATION  PATH  DEFINED  BY  HORIZONTAL  LAYERS 
568.000  C  OF  THE  ATMOSPHERE  OVER  WHICH  the  AIR  TEMPERATURE  AND 
569.000  C  RELATIVE  HUMIDITY  ARE  ASSUMED  TO  HAVE  The  AVERAGE  VALUES 
570.000  C  GIVEN  AT  THE  HEIGHTS  OF  THE  MIDPOINTS  OF  EACH  LAYER. 

571 .000  C 

572.000  DO  100  J=l,24 

5  7  '  . 0  00  C 

5  7  .  0  0  0  I  SBN  =  1  6  +  J 

5  7  5.0  0  0  FC=10.0::::  (ISBN/10.0) 

576.000  CALL  ANSAB( TKR, HRR, PAR, FC , AR ) 

5  7  ’  .  0  0  0  C 

5  7  8  .  0  0  0  BAPTE  =  BAAE  =  0  .  0 

5  7 a .0  0  0  DO  95  <  =  1,KLP 

588.000  CALL  ANSAB(ATK(<), AHR(K ), PA, FC, ATE) 

5  8  1.0  0  0  9  5  BAPTE  =  BAPTE  +  (ATE-AR  )-"DPTM(<  ) 

582.000  c 

58''.  000  DO  96  <  =  1,KLA 

5 . 0  0  0  CALL  ANSAB(ATkCK), AHR(<), PA, FC, ate) 

5^.0  00  96  BAAE  =  B  AAE  + (ATE-AR);:dAM(<  ) 

586.000  C 

58’. 000  IF(SPLPTM(J).EO.O. 0)GOTO  97 

533.000  SPLPTE(J)=SPLPTM(J)+8APTE 

5S9.000  GOTO  98 

590.000  97  SPLPTE(J)=0 . 0 

591.000  98  IF(SPLAM(j).EQ.O. 0)GOTO  99 

592.000  SPLAE( J )=SPLAM( J )+3AAE 

5  9  t  .  0  0  0  GOTO  100 

5  9U . 0  0  0  99  SPLAE(J)=0.0 

5  95.0  00  1  00  APME=APME  +  1  0  .  0::::((SPLA  E(J  )  +  AW(  J))/1  0  .  0  ) 

596.000  C 

597.000  C  CALCULATE  CORRESPONDING  PNL ,  PN  L T ,  AND  EPNL 

598 . 000  C 

599.000  CALL  CPNL(SPLPTB, PNLMB) 

600.000  CALL  P36TC(SPLPTB,CMXB, jmyb) 

601.000  FMXB=(FLOAT(FREQ(JMXB)))/inoO.O 

602.000  pnltmb=pnlmb+cmxb 

603.000  EPNLB=PNLTMB+DCF 

604.000  C 

605.000  CALL  CPNLCSPLPTC, PNLMC) 

606.000  CALL  P36TC(SPLPTC,CMXC, JMXC) 

607.000  FMXC=(FLOAT(FREO(JMXC)))/1000.0 

608.000  pnltmc=pnlmc+cmxc 

609.000  EPNLC=PNLTMC+DCF 

610.000  C 


611.000  call  cpnl( splptd, pnlmd) 

612.000  CALL  P36TC(SPIPTD,CMXD, JMXD) 

6  1  3  .  0  00  FMXDr(FLOAT(FREQ( JMXD)  ))  /  1  0  0  0  .  0 

614.000  pnltmd=pnlmd+cmxd 

615.000  EPNLDrPNLTMD+DCF 

616.000  C 

617.000  CALL  CPNLC SPLPTE, PNLME) 

".000  call  °3STC(SPLsTE,CMXE, JMXE) 

619.000  FMXE=(FLOAT(FREQ(JMXE)))/innO.O 

620.000  PNLTMEr PNLME+CMXE 

621.000  EPNLEr PNL TME+DCF 

622.000  C 

625.000  C  CALCULATE  CORRESPONDING  AL  AND  SEL 

624.000  C 

62  5  .  0  0  0  ALMB=  1  0  .  0::LOG1  0  (AP^B  ) 

626.000  SELB=ALMB+DSEL 

627.000  C 

62  8  .  000  ALMC  r 1 0 . 0  "‘LOG  I  0 ( APMC ) 

629.000  SELCrALMC+DSEL 

□30.000  C 

6  3  1.0  0  0  ALMD=  1  0  .  0  ::LOG  1  O(APMD) 

652.000  SELDrALMD+DScL 

633.000  C 

634.000  ALMErlO .  0::LOG10(APME) 

635.000  SELE=ALME+DSEL 

656.000  C 

□  3  7.0  0  0  END  OF  CALCULATIONS  :::: is :::::: 

638.000  C 

63}. ooo  c  set  sage  counters  and  initial  values  of  time  index  counters 
640.000  C 

64  1.0  00  LPGtOt r  INT((FLOAT(nS)  /  14 . 0  )  ■*•  0  .  99) 

6  42.0  00  I  Sr  1 

6  43.0  00  I Er  14 

6  44.0  30  L  PAGE  r  0 

6  45  .  0  00  IF(NS.GT.  150)  STOP  'NS  EXCEEDS  150' 

646.000  C 

64  7.000  C  LIST  THE  T£St-DAY  SPL  DATA  AND  THE  TE  ST -DA  f  NOISE  MEASURES 

643.000  C 

649.0  0  0  99  0  L  PAGE  r  L PAGE  -t- 1 

□50.000  WRITE  (6,995)  RUN, M I C , LPAGE, LPGTOT 

6  5  1.0  00  99  5  FORMA t (  1  i  1  , / , T 1 4 ,  'RUN  NO.:  ',13,';  MIC  NO.:  ',11, 

6  5  2.0  0  0  1  T  8  0 ,  'DATA  PAGE  ',12,'  OF  ’,12) 

653.000  WRITE  (6,1000)  T[TLE1 

5  5 4.  n  0  0  1  0  0  0  FORMAT^/, Tl,  '  |  , T  20  ,  2  0A4, T1  1  1,  '  |  ) 

6  5  5  .  0  00  WRITE  (6,1010)  t  i  tlE  2 

656.000  1010  FORMAT( T20 , 2 0A4) 

657.000  WRITE  (6,1020)  DA te , MONTH , YE AR , ZM, TS 

6  5  8.0  0  0  1  0  20  FORMAT(/,T27,  'TEST  date:  ' ,  12, 1  A 4 ,  ’  ’,12,';  ' 

6  5  9.0  00  l'MIC  HT;  • , F  3  -  1 ,  '  M;  START  TIME:  '  ,  F  7  .  1 ,  '  S') 

660.000  WRITE  (6,1030)  TOH , AH , A S , GR OUND 

6  6  1.0  0  0  1  0  30  FORMAT( T 1 7 ,  ' T I  ME  OH:  ',F7.1,'  S;  HT  OH :  * , F  6 .  1 ,  '  M;  ' 

66?  0  0  0  I  'AIRSPEED  OH:  ’,F4.1,'  M/S;  GR  ND  SURF:  >,5A4) 

663.000  WRITE  (6,1040) 

6  64.0  00  1040  FORMAT( / ; T 1 2 ,' CENTER ', 1 2X,  '  1/ 3-OCTAVE-BAND  SOUND  PRESSURE  ' 
6  65.0  00  1  'LEVEL,  DB  RE  20  MICROPASCALS') 

666.000  WRITE  (6,1050) 

6  6  7.0  0  0  1  0  5  0  FORMAT(Tl2,  'FREQ.,  '  ,27X,  'TIME  RELATIVE  TO  START  TIME,  S') 

668.000  WRITE  (6,1060)  (TT(I),[rIS,IE) 

6  69.0  0  0  10  60  FORMAT( T  14 ,  'HZ  '  , 2X, 14F6 . 1  ) 
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670 . 000  DO  1090  d  =  l, 24 

671.000  NJ=J-1 

672.000  IF(NJ. EO. 0)GOTO  1070 

673.000  IF(MOD(NJ, 3) . EQ. 0)GOTO  1070 

674.000  GOTO  1080 

675  .  000  1  0  70  WRITE  (6,  1075  )  ( FREQC J ) , ( S PL ( I , J  )  ,  I  =  I S ,  I E  )  ) 

6  76.000  1  0  75  F ORMAT( / , T 1 3 , I  5 , 14F 6 . 1 ) 

677.000  GOTO  1090 

678.000  1080  WRITE  (6  ,  1085  )  (FREQ( J  )  , ( S PL ( I , J  )  ,  I  =  I S ,  I  E  )  ) 

6  79  -  000  1  08  5  F ORMAT (T 1 3 ,  I  5 ,  1 4F 6  .  1  ) 

680.000  1090  CONTINUE 

681.000  WRITE  (6,1100)  (PNL ( I  )  ,  I  =  I S ,  I  E  ) 

682.000  1  1  00  FORMAT(/, T12,  ' PNL, DB  '  ,  14F6.  1) 

683.000  WRITE  (6,1110)  (PNLT (I), I = I S ,  IE) 

684.000  1  1  1  0  FORMAT(T  1 1,  ' PNLT, DB ' , 14F6 .  1 ) 

685.000  WRITE  (6,1120)  (AL( I ), I = I S , I E ) 

686  .  000  1  1  20  FORMAT(Tl,  '  I  ' ,T13,  'AL,DB ' ,  14F6.  1, Till,  1  |  '  ) 

687.000  WRITE(6,  1122)  (FMX( I ) , I  =  I  5, I E) 

688.000  1122  FORMAT(Tll, 'FMX,KHZ ' , 14F6. 3) 

689.000  WRITE(6, 1124) 

690.000  1124  FORMAT(T18, 'FMX  IS  THE  CENTER  FREQUENCY,  IN  KHZ,  OF  THE  ' 

691.00  0  1'1/3-OCTAVE  BAND  THAT  PRODUCES  C(  I  )  1  ) 

692  .  000  WRITE(6,  1126) 

693  .  000  1  1  26  F  ORMAT ( T 1 2 ,  ' C ( I )  IS  THE  LARGEST  OF  THE  TONE -COR REC T I  ON  1 

694.000  1  1  F  AC  T  OR  S  FOR  THE  CORRESPONDING  ' 

695.000  2 ' SPECTRUM  SPL(I,d)') 

696.000  IF(IE.NE.NS)GOTO  1150 

697.000  WRITE  (6,1130)  PNLM, TPNLM , PNLTM , TPNLTM, DC F , TT ( I F P ) , TT ( I LP ) 

698.000  1  1  30  FORMAT(/, Til,  ' PNLMr  •  , F  5 .  1 ,  '  DB  AT  ',F4.1,'  S;  PNL  TM  = 

699.000  1 F  5 .  1 ,  1  DB  AT  \F4.1,’  S;  DC  F  =  \F4.1,'  DB  FROM  «,F4.1, 

700.000  2  '  TO  '  ,F4.  1,  1  S') 

7  01  .  000  WRITE  (6,  1140)  E  PNL 

702.000  1  1  40  FORMAT(Tl,  •  I  1  , T 1 1 ,  1 E PNL =PNL TM+DC F =  * , F  5 .  1 ,  ’  DB ’ , T 1 1 1 ,  •  I  1  ) 

703  .  000  WRITE  (6,1145)  ALM , TALM , SE L , TT ( I F A ) , TT ( I  LA  ) 

7  04.000  1  145  FORMAT(Tll, ’ALM=  ’  , F  5 . 1 ,  '  DB  AT  • , F  4 . 1 ,  *  S;  SEL=  ’,F5.1, 

705.000  1’  DB ;  TSEL  FROM  ',F4.1,'  TO  ’,F4.1,'  S’) 

706. 000  C 

707.000  C  REITERATE  FOR  NEXT  PAGE  UNTIL  GET  TO  LAST  PAGE  FOR  INPUT  DATA 

708.000  C 

709.000  IF(IE.EQ.NS)GOTO  1800 

7  10.000  1  1  50  I S=M I N  0 ( l E  +  1 , NS  ) 

711.000  IE=MINQ(IE+14,NS) 

712.000  IF(IE.LE.NS)GOTO  990 

713.000  C 

714.000  C  TABULATE  THE  RESULTS  OF  APPLYING  VARIOUS  METHODS  FOR  ADJUSTING 

715.000  C  SOUND  PRESSURE  LEVELS  FOR  DIFFERENCES  IN  ATMOSPHERIC  ABSORPTION 

716.000  C  UNDER  TEST  AND  REFERENCE  METEOROLOGICAL  CONDITIONS 

717.000  C 
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718.000  C 
719. 000  C 
720.000  1800 
721.000  1810 
722.000 
723. 000 
724.000  1820 
725.000 
726 . 000 
727.000  1830 
728 . 000 
729.000 
730.000  1840 
731.000 
7  32  .  000 
733.000  1850 
734. 000 
735.000 
736.000  1855 
737. 000 
738.000 
739.000  1857 
740 . 000 
741.000 
742.000  1858 
743.000 
744.000  1860 
745.000 
746.000  1870 
747 . 000 
748.000 
749.000  1880 
750.000 
751.000 
752.000 
753.000  1890 
754. 000 
7  55 . 000  C 
756.000  C 
757.000  C 
758.000 
759. 000 
760.000 
761.000  2000 
762 . 000 
763.000  2010 
764. 000 
765.000 
766.000  2020 
767.000 
768. 000 
769.000  2030 
770 . 000 
771.000  2040 
772. 000 


SUMMARY  OUTPUT  PAGE  1  FOR  METEOROLOGICAL  DATA 
WRITE  (6, 1810)  RUN, MIC 

FORMATC  1  *,  5/,  Tl,  '  |  ’  ,T17,  'RUN  NO.:  ',13,';  MIC  NO.:  ',11, 

1T48,  'METEOROLOGICAL  DA TA 1 , T 76 ,  *  PAGE  1  OF  SUMMAR Y • , 1 8X ,  1  I  • ) 

WRITE  (6,1820) 

FORMATC//, T16, 'HEIGHT  AT  MEAS.  TEMP.  MEAS .  R.H.  MEAS.' 

1'  MOLAR  REF.  TEMP.  REF.  R.H.  REF.  MOLAR') 

WRITE  (6  ,  18  30  ) 

FORMAT(T 15, 'MIDPOINT  OF  OVER  HEIGHT  OVER  HEIGHT  CONC . ' 

1  '  OF  AT  MID.  OF  AT  MID.  OF  CONC.  OF') 

WRITE  (6,1840) 

FORMAT(Tl6, ' INTERVAL,  INTERVAL,  INTERVAL,  WATER  ' 

1 'VAPOR,  INTERVAL,  INTERVAL,  WATER  VAPOR,') 

WRITE  (6,1850) 

FORMATCT 20, ' M  C  PCT.  PC T . ' 

1  '  C  PCT.  PCT .  '  ) 

WRITE  (6,1855) 

FORMAT  (/) 

IF(HM( 1 ) .GT . ZM)GOTO  1857 
GOTO  1858 

WRITE(6,l860)((HMI(K),ATK(K)-273. 1 5 , AHR (K ) , HMEA(K ) , 

1 R  TK  (K  )-  2  7  3  . 1 5 , RHR (K ), HREF( <)),<= 1 , NH) 

GOTO  1860 

WRI TE(6, 1860 )((HMI (K),ATK(K)-273 . 1 5 , AHR (K) , HMEA( K) , 
lRTK(K)-273. 1 5 , RHR (K ) , HR EF ( K ) ) , K = 1 ,NH-1 ) 

FORMATCT 17, F 6. 1,9X,F4.1,9X,F4.1,8X,F5.3,9X,F4.1,9X,F4.1,8X,F5.3) 
WRITE  (6  ,  18  70  )  PA;:  1  0  1  .  3  2  5  ,  PA 

FORMAT  (/,T24, 'STATION  BAROMETRIC  PRESSURE  =  ',F7.3, 

1  '  KPA  (  ' ,F  5. 3,  '  STD.  ATM  ) ' ) 

WRITE  (6,1880)  TC10,HR10,H10 

FORMATC/, Tl, ' | ' ,T12, 'AT  10  M:  MEAS.  AIR  TEMP.=  ',F4.1, 

1'  DEG  C;  MEAS.  REL.  HUM . =  ',F4.1,'  PCT;  MEAS.  MOLAR  CONC.' 

2 '=  ' , F  5 . 3,  '  PCT'  , T 1 1 1 ,  '  I  ' ) 

WRITE  (6, 1890) 

FORMATCT23, 'REF .AIR  TEMP.=  25.0  DEG  C;  REF.  REL.  HUM.' 
l's  70.0  PCT;  REF.  MOLAR  CONC . =  2.18774  PCT') 

SUMMARY  OUTPUT  PAGES  2  AND  3  FOR  NOISE  DATA 

DO  2300  NPG  =  1,2 

IFCNPG.EQ. 1)WRITE(6, 2000 )R UN, MIC 
IF (NPG. EQ. 2)WRITE(6, 2010)RUN,MIC 

FORMATC ' 1  ' , 2/, Tl,  '  I  ' , T14,  'RUN  MO.:  ',13,';  MIC  NO.:  ',11, 

1 T  5 1 ,  ' PNL  AND  EPNL ’ , T 76 ,  ' PAGE  2  OF  SUMMAR Y ' , 1 8 X  ,  '  I  '  ) 

FORMATC ' 1  ' , 3/, Tl,  '  I  ' , T14,  'RUN  NO.:  ',13,';  MIC  NO.:  ',11, 

1  T  5  2 ,  ' AL  AND  SEL '  ,  T 76  ,  ' PAGE  3  OF  SUMMAR Y ', 1 8X,  '  I  ' ) 

WRITE  (6,2020) 

FORMATC/, T21, 'TEST-DAY  NOISE  DATA  ADJUSTED  TO  ' 

1 'ACOUSTICAL  REFERENCE-DAY  CONDITIONS') 

WRITE  (6,  2030  )  TITLE2 
FORMATCT20, 20A4) 

WRITE  (6,2040) 

FORMATCT19, 'REF . -DAY  COND.:  AIR  TEMP  =  25.0  C;  REL  HUM  =  ' 

1 '70.0  PCT;  PRESS  =  1.0  STD  ATM') 
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773.000  WRITE  (6,2050) 

7  7  4.0  00  20  50  FORMAT ( / , T  1  3 ,  ' BAND' ,3X,'(A)',4X,'(B)',4X/,CC),,4X,'(D)', 

7  7  5.0  00  14X,  '  (E)  '  ) 

776.000  WRITE  (6,2060) 

777.000  2060  FORMAT ( T\2,  ' CENTER1, 2X,  'TEST', 3X,  '  AD J.  '  , 3  X ,  ' AD J .  ' , 3X,  '  AD J. 
778.000  1  3  X ,  ' AD  J  '  ) 

7  79.0  00  WRITE  (  6,  20  70  ) 

780.000  2070  FORMA  T  (  T  I  2  ,  '  FR  EO  .  ,  '  ,  2  X  ,  '  SPL  ,  '  ,  3  X  ,  '  SPL  ,  '  ,  3  X  ,  '  SPI_  ,  '  ,  3  X  ,  '  SPL  , 
7  8  1.0  0  0  1  3X,  ' SPL,  '  ) 

732.000  WRITE  (6, 2080) 

783.000  2080  FORMA  T ( T  1 4 ,  '  HZ ' , 5 X,  ' OB  '  , 5 X,  • OB  '  , 5 X ,  '  DB  '  , 5 X,  '  OB  '  , 5 X ,  •  OB  '  ,  /  ) 

7  84.0  00  J  J  =0 

735.000  DO  2160  J=l,24 

736.000  JJ=JJ+1 

737.000  C 

733.000  C  PRINT  OUTPUT  DATA  AND  A  LINE  OF  text 

739.000  C 

790.0)0  IF(NPG.EO. 1 )GOTO  2110 

7  91  .  0  00  I F (u  .  EO  .  3 . AND .NPG .EO . 2 )GOTO  2130 

7  92  .  0  0  0  I  c (NPG . EO . 2 )GOTO  2120 

7  9  3.0  0  0  2  1  00  FORMAT( T13,  i 5, F6 .  1, 4F7 . 1 , T58, 1 2AA) 

794.000  2110  WRITE  (6,2100)  (FRE0(J),SPLPTM(J),SPLPT3(J),SPLPTC(J), 
795.00  0  1  SPLPTD( J), SPLPTE(u)), (TEXT( I  I , jj),  I  I  =  1,  12  ) 

796.000  GOTO  2140 

7  9  7  .  0  0  0  2  1  2  0  WR I TE ( 6 , 2 1 0 0 )( FR EO( J ) , SPLAM( J ) , S PLA8 ( J ) , $ pL AC ( J  )  , 

7  93.0  00  1  SPLAD( J) , SPLAE(J) ), (T£XT( I  I  , jj)  ,  I  I  =  1  ,  12  ) 

799.000  GOTO  2140 

300.000  2130  WRITE(6,2100)(FREQ(J),SPLAM(J),SPLAB(J),SPLAC(J), 

8  0  1.0  0  0  1  SPLAD(J),SPLAE(J)),(TEXTA( I  I ),  I  I =1, 12) 

802.000  2140  IF(MOD( J , 3 ) .NE . 0 )GOTO  2160 

3  0  3.0  0  0  I F ( j . EO . 2  4 )GOTO  2170 

804.000  C 

80  5  .  0  00  C  D  R I N  T  A  LINE  OF  TEXT  ONLY 

806.000  C 

807.000  JJ=JJ+1 

8  08.0  00  WRITE  (6,2150)  (TEXT( I  I ,JJ), I  I =1,  12) 

3  09.0  00  21  50  FORMAT( T58 , i 2A4  ) 

310.000  2160  CONTINUE 
311.000  2170  I F (NPG . EO • 1 )GOTO  2180 
812.000  IF(NPG.EQ.2)GOTO  2240 

815.000  C 

814.000  C  BOTTOM  OF  PAGE  2 

815.000  C 

316.000  2130  WRITE  (6,2190)  (PNLM, PNLMB , PNLMC , PNLMD, DNLME ) 

8  1  7.0  00  21  90  FORMAT(//, T  12,  'PNL,DB >,F6.  1,4F7.  1  ) 

8  1  8  .  0  00  WRITE  (  6,  2  20  0  )  (PNLTM, PNLTMB , PNLTMC , PNLTMD , PNL  THE  ) 

8  19.0  00  22  00  FORMAT(Tll/  >PNLT,DB',F6.  1,4F7.  1) 

820.000  WRITE  (6,2210)  ( EPNL , E PNLB , E PNLC , E PNL D , E PNL E ) 

8  2  1  .  0  00  22  1  0  FORMAT( T  1  1  ,  ' EPNL , DB  '  ,F6 .  1 , 4F7  .  1  ) 

8  22  .  000  WRITE  (6,2215)  (FMXB , FMXC , FM XD , F MX E  ) 

3  2  3.0  0  0  22  1  5  FORMAT ( T 1 1 ,  '  FMX, KHZ ' , T 24 , 4F 7 . 3  ) 

824.000  WRITE  (6,2220)  AM H 

8  25.0  00  22  20  FORMAT (//,T21,  'AIRPLANE  HEIGHT  OVER  MICROPHONE  ' 

3  26  .  0  0  0  1  '  =  '  ,F6  .  1 ,  '  M  '  ) 

827.000  WRITE  (6,2230)  PDPTM 

828.000  2230  FORMAT(Tl/ ' | ',T21, 'SOUND  PROPAGATION  PAThlENGTh' 

8  29.000  1  '  FOR  TEST  SPLS  AT  PNLTM  =  * , F  6 .  1  ,  '  KI,T111,,|') 

3  30.0  00  WRITE  (  6,  2  2  3  5  )  (PS  I P tm»  5  7  .  2  9  5  7  7 9  8  1  ) 

8  3  1  .  0  00  22  3  5  FORMAT( T2 1 ,' SOUND  PROPAGATION  ANGLE  RELATIVE  to  ' 

8  3  2  .  0  00  1'  DIRECTION  OF  FLIGHT  =  » , F  5 .  1 ,  ’  DEG  AT  PNLTM') 


GOTO  2500 


835.000 
8  34 . 0  00  C 
855.000  C 
8  36 . 000  C 
837.000  2240 
838.000  2250 
8  39 . 0  00 
840.000  2260 
841 .000 
842.000  2270 
843.000 
844.000 
845.000  2280 
346.000 
847.000 
848.000  2290 
849.000 
8  5  0.0  00  2  3  0  0 
851.000  C 
852.000  C 
353.000  C 
354.000 
355.000 
356.000 
857.000 
358.000  C 
359.000  C 


BOTTOM  OF  PAGE  3 

WR  1  TEC  6,  2  2  50  )C ALM, ALMB , ALMC, ALMD, AlME  ) 

FORMA  T(  //,  T  1  3,  '  A|_  ,  DB  '  ,F6.  1 ,4F7  .  1  ) 
WRITEC6,2260)CSEL,SELB,SELC,SELD,SELE) 

FORMA T( T 12, 'SEL,D8',F6.1,4F7.1) 

WR ! T£(  6,  2  2  70  )  AMH 

FORMAT(//^T21, 'AIRPLANE  HEIGHT  OVER  MICROPHONE  ' 
1 '=  ' ,F6 . 1, '  M' ) 

WR I te(  6,  2  2  80  )  POAM 

FORMATC T l ,  <  |  1 ' T 2 1  ,  ' SOUND  PROPAGATION  PAThLENGTh ' 
1*  FOR  TEST  SPLS  AT  ALM  =  ',F6.1,'  M',T11[/|') 

WRITE  C  6 ,  2  2  9  0  )  CPSI AM"  57. 29577951) 

FORMAT^ T2 1 SOUND  PROPAGATION  ANGLE  RELATIVE' 

1'  TO  DIRECTION  OF  FLIGHT  :  ',F5.1,'  DEG  AT  ALM' 

CONTINUE 

READ  NEXT  SET  OF  INPUT  DATA.  STOP  IF  RUN  NO.  = 

READ  C  5,  1 0 )  RUN 
IFCRUN.NE.999)GOTO  11 
STOP  'END  OF  INPUT  DATA' 

END 

END  OF  STATEMENTS  FOR  MAIN  PROGRAM. 


99  R  . 


-.0 


BLOCK  DATA 


The  labeled  common  statements  in  BLOCK  DATA  contain,  in  384  computer 
words,  the  text  for  31  lines  on  output  summary  pages  2  and  3.  Some  of  the 
lines  of  text  are  blank.  The  words  of  the  text  describe  the  legends  for 
the  columns  labeled  (A),  (B) ,  (C) ,  (D),  and  (E)  of  measured  test-day  and 
reference-day  1/3-octave-band  sound  pressure  levels. 
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LISTING  OF  BLOCK  OATA  FOR  CODES  ON  OUTPUT  SUMMARY  PAGES 


860.000  C 
861 . GOO  C 
86?. 000  C 
863.000  C 
864.000  C 
865.000 
866.000  C 
867.000 
868 . 000 
869.000 
870.000 
871 . 000 
872.000 
873.000 
874.000  C 
875.000 
876.000  C 
877.000 
878. 000  C 
879.000 
880 . 000 
881.000 
882.000 
883. 000 
884. 000 
885.000 
886.000 
887. 000 
888.000 
889.000 
890.000 
891 . 000 
892.000  C 
893.000 
894. 000  C 
895.000 


BLOCK  DATA 

COMMON  /CTEXT/  TEXT  1  Cl  2 ) , TEXT  2 C 1 2 ) , TEXT  3 ( 1 2 ) , TEXT 4( 1 2 ) , 

1TEXT 5(12), TEXT 6(12), TEXT7(12), TEXT  8(12), TEXT9(12), TEXT  10(12), 
2TEXT 11(12 ), TEXT  12 (12), TEXT  13(1 2), TEXT  14(1 2), TEXT  15(1 2), 

3TEXT 16(1 2), TEXT  17 (12), TEXT  18 (12), TEXT  19(1 2), TEXT  20 (12), 
4TEXT21(12), TEXT  22(12), TEXT23(12), TEXT  24(12), TEXT25(12), 

5TE XT  26(12), TEXT27(12), TEXT  28(12), TEXT  29(12), TEXT  30(12), 
6TEXT31(1 2), TEXTA(1 2) 

DIMENSION  TEXT (12, 31) 

EQUIVALENCE  (TEXT ( 1 , 1 ), TEXT  1 (1 )) 


DATA  TEXT/ 3  7 2 (4H  )/ 

DATA  TEXT  1 / '  COLUMN  CODES:  ’  / 

DATA  TEXT  3/' (A)TEST -DAY  SPLS  AT  TIME  OF  TEST-DAY  PNLTM  +  '/ 

DATA  TEXT4/ '  ASSOCIATED  INTEGRATED  NOISE  MEASURES.  '/ 

DATA  TEXT6/ 1  THEN,  TEST-DAY  SPLS  ADJUSTED  TO  ACOUSTICAL  '/ 

DATA  TEXT  7/ ' REFERENCE -DAY  CONDITIONS  +  ASSOCIATED  '/ 

DATA  TEXT 8/ ' MEASURES .  ADJUSTMENTS  MADE  USING:  '  / 

DATA  TEXT  1 0 / '  (B)10-M  MET.  DATA  AND  SAE  ARP  866A;  '/ 

DATA  TEXT  1 1 / 1  (C)LAYERED  MET.  DATA  ALOFT  AND  SA£  ARP  866A;  '/ 

DATA  TEXT  1 2 / '  (D)LAYERED  MET.  DATA  ALOFT  AND  ABSORPTION  BY  '/ 

DATA  TEXT  1 3 / 1  ANS  SI. 26-1978  INTEGRATED  OVER  PASSBAND;  '/ 

OATA  TEXT  14/ 1  (E)LAYERED  MET.  DATA  ALOFT  AND  ABSORPTION  BY  ■/ 

DATA  TEXT  1 5 / f  ANS  SI. 26-1978  AT  BAND  CENT.  FREQ.  ONLY.  '/ 

DATA  TEXTA/ ' (A)TEST -DAY  SPLS  AT  TIME  OF  TEST-DAY  ALM  +  '/ 


END 
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SUBROUTINE  CPNL(SPL,  PNL) 


Subroutine  CPNL  accepts  an  array  of  24  1/ 3-octave-band  sound  pressure 
levels  and  calculates  the  corresponding  perceived  noise  level,  PNL.  It 
requires  about  355  computer  words  or  approximately  1400  bytes  of  storage 
capacity. 
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SUBROUTINE  CPNL (S PL , PNL  ) 


896.000  C 
897 . 000 
898 . 000  C 
899. 000  C 
900.000  C 

901.000  C  SPL  IS  INPUT  ARRAY  OF  1 /3 -OCTAVE -BAND  SPLS  FROM  50  TO  10,000  HZ 
902.000  C  FOR  WHICH  PNL  IS  TO  BE  DETERMINED;  DB . 

905 . 000  C 

904.000  C  PERCEIVED  NOISINESS,  PN,  IN  NOYS  IS  CALCULATED  FOR  THE 

905.000  C  CORRESPONDING  1 /3-OCTAVE-BAND  SOUND  PRESSURE  LEVEL,  SPL,  AND 
906.000  C  BAND  INDEX  COUNTER,  J,  WHERE  U  RANGES  FROM  1  TO  24  FOR  THE 
907.000  C  NOMINAL  BAND  CENTER  FREQUENCIES  FROM  50  TO  10,000  HZ. 

908. 000  C 

909.000  C  THE  BAND  PERCEIVED  NOISINESS  IS  SET  EQUAL  TO  0.0  NOYS  FOR  ANY 

910.000  C  BAND  SPL  LESS  THAN  THE  CORRESPONDING  CONSTANTS  LI  IN  THE 

911.000  C  DATA  LIST. 

912. 000  C 

913.000  C  IF  ALL  SPL(U)  ARE  <  L1(J),  THEN  THE  PNL  IS  SET  =  0.0  DB . 

914.000  C 

915.000  C  THE  COMPUTATION  ASSUMES  THAT  NO  BAND  SOUND  PRESSURE  LEVEL 

916.000  C  IS  GREATER  THAN  150.0  DB . 

917.000  C 

918.000  C  DATA  CONSTANTS  ARE  FROM  SAE  ARP  865A,  15  AUGUST  1969. 

919.000  C 

920.000  C  PNL  IS  THE  PERCEIVED  NOISE  LEVEL  CORRESPONDING  TO  THE  ARRAY 
921.000  C  OF  SPLS;  DB . 

922. 000  C 

92  3.000  "• :: :: :: :: ::  5: :: :: :: :: :: ::  ”  -  »  » :: :: »  ~  "  -  -  » "•  -  ~  »  » :: -  »  »  -  "  ;t  ::  -  ”  -  -  :: - »  »  -  ~ 

924. 000  C 

925.000  REAL  L 1 , M 1 , L 2 , M 2 , L 3 , M 3 , LC , L4 , M4 

926.000  DIMENSION  L 1 ( 2 4 ) , M 1 ( 2 4 ) , L 2 (2 4 ) , M 2 ( 2 4 ) , L 3 C 2 4 ) , M 3 C 2 4 ) , 

927.000  1LC(24), L4(24),M4(24), SPL(24) 

928.000  C 

929.000  DATA  Ll/49.,44.,39.,34.,30.,27.,24.,21.,18.,16.,16.,16.,16.. 

930.000  116. ,15. ,12. ,9. ,5. ,4. ,5. ,6. ,10. ,17. ,21./ 

931.000  DATA  Ml/. 079520, .068160, .068160, .059640, .053013, .053013, 

932.000  1.053013, .053013, .053013, .053013,  .053013, .053013, .053013, 

933-000  2. 053013, .059640, .053013, .053013, .047712, .047712, .053013, 

934.000  3-053013, .068160, .079520, .059640/ 

935.000  DATA  L2/55., 51. ,  46  . ,  42  . ,39. ,36. ,33. ,30. ,27. ,25. ,25. .25. ,25., 

936.000  125. ,23. ,21. ,18. ,15. ,  14. ,14. ,15. ,17. ,23. ,29./ 

937.000  DATA  M2/. 058098, .058098, .052288, .047534, .043573, .04 3573, 

938.000  1.04 0221, .037349, .034859,  .034859,  .054859,  .034859,  .034859, 

939-000  2. 034859, .034859, .040221, .037349,  .034859, .034859, .034859, 

940.000  3.034859, .037349, .037349, .043573/ 

941.000  DATA  L3/64., 60. ,56. ,53. ,51. ,48. ,46. ,44. ,42. ,40. ,40. ,40. ,40., 

942.000  140. ,38. ,34. ,32. ,30. ,29. ,29. ,30. ,51. ,37.  ,41./ 

943.000  DATA  M3/.  043478,  .040570,  .036831,  .036831,  .035336,  .033333, 

944.000  1.033333, .032051, .030675, .030103, .030103,  .030103, .030103, 

945.000  2. 030103, .030103, .029960, .029960, .029960,  .029960, .029960, 

946.000  3.029960, .029960, .042285, .042285/ 

947.000  DATA  LC/91. 01,  85. 88,  87. 32,  79. 85,  79. 76,  75. 96, 73. 96, 74. 91, 94. 65, 

948.000  1100. ,100.  ,100. ,100. ,100. ,100. ,100. ,100.  ,100. ,100. ,100.  ,100., 

949  .000  2  100  . ,  44. 29  ,  50  .  72  / 
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950.000  DATA  M4/. 030103, .030103, .030103, .030103, .030103, .03010}, 

951.000  1.030103, .030103, .030103, .030103, .030103, .030103, .030103, 

952.000  2. 030103,. 030103,. 029960,. 029960,. 029960,. 029960,. 029960, 

953.000  3.029960, .029960, .029960, .029960/ 

954.000  DATA  L4/52., 51. ,49. ,47. ,46. ,45. ,43. ,42. ,41. ,40. ,40. ,40., 

9  5  5.000  140  . ,  40. ,  38. ,  34. ,  32. ,  30. ,  2  9  . ,  29. ,  30. ,  3  1 . ,  34. ,  3  7  .  / 

956.000  C 

957.000  PNSUM -0 . 0 ;  PNMX=-10.0 

958.000  DO  80  J=l,24 

959-000  IF(SPL(J).LT.Ll(J))GOTO  5 

960.000  GOTO  7 

961.000  5  PN=0.0 

962.000  GOTO  75 

963.000  7  IF(SPL(J).GE.L1(J). AND . SPL(U).LT.L2(J))  GOTO  10 

964.000  GOTO  20 

965.000  10  PN  =  0.1”(10.5!::((Ml(J))::(SPL(J)-Ll(J)))) 

966.000  GOTO  75 

967.000  20  IF(SPL(J).GE.L2(J).AND.SPL(J).LT.L3(J))  GOTO  30 

968. 000  GOTO  40 

969.000  30  PN  =  10."*::(CM2(J))“(SPLCJ)-L3CJ))) 

970.000  GOTO  75 

971.000  40  IF(SPL(J).GE.L3(J).AND.SPLCJ).LT.LCCj))  GOTO  50 

972.000  GOTO  60 

973.000  50  PN  =  10.““((M3(J))::(SPLCvJ)-L3CJ))) 

974.000  GOTO  75 

975.000  60  IFCSPL(J).GE.LCCJ).AND.SPlCO).LE.  150.)  GOTO  70 

976.000  70  PN=10.«"CCM4(J))»CSPLCJ)-L4(J))) 

977.000  75  PNMX=AMAX(PNMX,PN) 

978.000  PNSUMrPNSUM  +  PN 

979.000  80  CONTINUE 

980.000  PNT  =  0  .  85 ::PNMX  +  0.15::PNSUM 

981.000  IF(PNT)90, 90, 85 

982.000  85  PNL=40.0  +  33  .  22  :!LOG  1  0  (PNT  ) 

983.000  RETURN 

984.000  90  PNL  =0 . 0 

985.000  RETURN 

986.000  END 


SUBROUTINE  P36TC(SPL,  CMX,  JMX) 


Subroutine  P36TC  accepts  an  array  of  24  1/ 3-occave-band  sound  pressure 
levels  and  calculates  a  tone-correction  penalty,  CMX,  appropriate  for  that 
spectrum  using  a  procedure  consistent  with  the  requirements  of  the  April 
1978  version  of  Part  36  of  the  Federal  Aviation  Regulations.  The  index 
counter,  JMX,  for  the  1/ 3-octave-band  producing  the  largest  tone  correc¬ 
tion  factor  is  also  found  and  returned.  If  certain  requirements  are  not 
satisfied,  no  tone-correction  factor  is  calculated  (CMX  ■  0.0  d B ) 
and  a  value  of  JMX  *  1  is  returned  to  indicate  the  impossible,  by  defini¬ 
tion,  result  of  a  tone  in  the  band  at  50  Hz  as  a  warning  that  no  tone 
correction  was  calculated. 

P36TC  requires  about  400  words  or  about  1650  bytes  of  storage  for  the 


source  code. 


987.000  C 

988.000  SUBROUTINE  P 3 6 TC ( S PL , CMX , JM X ) 

989. 000  C 

99  0.000  C *'•' ::  •"  •" :: •" ::  •“ ::  •” ::  ” :: :: :: :: :: :: :: :: :: ::  •" :: :: :: :: :: :: :: :: :: ::  ::  ”• :: :: :: :: :: ::  "• ::  "• ::  ”  - ::  - :: ::  -  -  - 

991.000  C 

992.000  C  SPL  IS  THE  INPUT  ARRAY  OF  24  SPL  VALUES  FOR  1/3-OCTAVE  BANDS 
993.000  C  FROM  50  TO  10,000  HZ  AT  THE  I-TH  INSTANT  OF  TIME  DURING 

994.000  C  A  FLYOVER;  DB . 

995.000  C 

996.000  C  CMX  IS  THE  LARGEST  VALUE  OF  THE  TONE -C ORREC T I  ON  FACTORS 
997.000  C  CALCULATED  BY  SUBROUTINE  P56TC;  DB 

998.000  C 

999.000  C  JM X  IS  THE  FREQUENCY  BAND  INDEX  COUNTER  FOR  THE  1/3-OCTAVE 

1000.000  C  BAND  PRODUCING  THE  VALUE  OF  C  IDENTIFIED  AS  CMX  BY  THE 

1001 . 000  C  SUBROUTINE. 

1002 . 000  C 

1  0  0  3  0  0  0  C:: :::: :: :: :: :: :: :: :: :: :: :: :: :: :: :: :c :: :: :::::: :: :: :: :: 
1004.000  C 

1005.000  C  SUBROUTINE  P 3 6TC( S PL , CM X , JM X )  EXAMINES  A  SET  OF  1/3-OCTAVE- 

1  006.00  0  C  BAND  SOUND  PRESSURE  LEVELS  AND  DETERMINES  A  TONE -CORREC T I  ON 
1007.000  C  PENALTY  WHICH  IS  TO  BE  ADOED  TO  THE  PERCEIVED  NOISE  LEVEL  FOR  THE 
1008.000  C  SPL  SPECTRUM  TO  GIVE  THE  TONE -CORREC TED  PNL ,  OR  PNLT. 

1009.000  C 

1010. 000  C  TONE-CORRECTION  FACTORS  FOR  THE  PRESENCE  OF  PRONOUNCED 

1011.000  C  IRREGULARITIES  IN  THE  SPL  SPECTRUM  ARE  CALCULATED  BY  THE 
1012.000  C  PROCEDURE  GIVEN  IN  SECTION  836.5  OF  APPENDIX  B  TO  FAR  PART  36, 

1013.000  C  CHANGE  8,  EFFECTIVE  3  APRIL  1978,  WITH  INTERPRETATIONS  DESCRIBED 
1014. 000  C  BEuOW . 

1 0 1 5 . COO  C 

1015.  COO  C  THE  REQUIREMENT  OF  A36.  5(D)(3)  OF  CHANGE  8  OF  PART  36  THAT  NO 

1017.000  C  E  PNL  BE  CALCULATED  OR  REPORTED  IF  ANY  SPECTRUM  WITHIN  THE 
1018.000  C  1 0 -DB -DOWN  DURATION  HAS  MORE  THAN  FOUR  MISSING  BAND  LEVELS  OUT 
1019.000  C  OF  24  IN  THE  RANGE  OF  BAND  CENTER  FREQUENCIES  FROM  50  TO  10,000  HZ 
1020.000  C  IS  SPECIFICALLY  NOT  INCORPORATED  HERE  BECAUSE  TO  HAVE  DONE  SO  WOULD 
1021.000  C  HAVE  ELIMINATED  THE  PNLT  AND  EPNL  CALCULATIONS  FOR  SOME  OF 
1022.000  C  THE  AIRCRAFT  NOISE  MEASUREMENTS  AVAILABLE  FOR  THE  STUDY  OF 
1023.000  C  ATMOSPHERIC  ABSORPTION  ADJUSTMENTS  FOR  AIRCRAFT  NOISE  DATA. 

1 02  A . 000  C 

1025.000  C  PARAGRAPH  B36.5(M)  OF  PART  36,  CHANGE  8,  PERMITS  THE  EXCLUSION 

1  0  2  6.00  0  C  OF  TONE -CORRECT  I  ON  PENALTIES  RESULTING  FROM  SPECTRAL  IRREGULARITIES 
1  0  2  7.  00  J  C  (OR  PSEUDOTONES)  CAUSED  BY  GROUND-REFLECTION  EFFECTS  IN  THE  FREQUENCY 
1028.000  C  RANGE  COVERED  BY  THE  1/3-OCTAVE  BANDS  WITH  CENTER  FREQUENCIES  OF 
1029.000  C  800  HZ  AND  LOWER  (I.E.,  BAND  NO.  13  AND  LOWER). 

1050.000  C 

1051.0ii0  C  SINCE  MOST  AVAILABLE  AIRCRAFT  FLYOVER  NOISE  DATA  WERE  MEASURED 

1052.000  C  By  MICROPHONES  LOCATED  1.2  M  ABOVE  THE  GROUND  PLANE,  MOST  OF 
1055. OuO  C  Tut  MEASURED  SOUND  PRESSURE  LEVELS  WERE  KNOWN  TO  CONTAIN  LOW- 
1058.000  C  FREQUENCY  SPECTRAL  IRREGULARITIES  CAUSED  BY  GR OUND -R EF L EC T I  ON 

1 0  3  F . 0  G  0  C  EFFECTS. 

1056. 000  C 

1  0  5  7  .  0  0  0  C  BECAUSE  THERE  WAS  NO  CONVENIENT  WAY  TO  ELIMINATE  THE 

1058.000  SPECTRAL  IRREGULARITIES  CAUSED  BY  GROUND  REFLECTION  EFFECTS 
10  5  1.  0 0  0  C  AND  BECAUSE  IT  WAS  NECESSARY  TO  HAVE  A  CONSISTENT  DETERMINATION 
1)80.00  1  OF  The  r  F  f  E T  S  OF  DIFFERENT  METHODS  OF  ACCOUNTING  FOR 

om  ;  a-mcspmlric  absorption  on  tone -corrected  pnls  and  epnls, 


1 042 . 000 
1043.000 
1 044. 000 
1043.000 
1046.000 
1 047 . 000 
1 048 . 0O0 
1 U49 ■ 0  0  0 
1050.000 
1 051 . 000 
1052.000 
1053. 000 
1 054. 000 
1055.000 
1056. OuO 
1057.000 
1058.000 
1059. OOU 
1060.000 
1061.000 
1062.000 
i 063 . 000 
1064.000 
1065.000 
1 o 06  .OuO 
1 0u7 . UOO 
1 068 . 000 
1069.000 
1070.000 
1071 .000 
1072.000 
1073.000 
1 074 . UOO 
1 075 . OuO 
1076.000 
1077.000 
1078.000 
1079.000 
1 080 . OUO 
1 08 1 . UOO 
1082.000 
1083.000 
1 U84 . O00 
1 085 . 000 
1086.000 
1  0  8  7  .  0  u  0 
1 088 . UOO 
1 u89 . 000 
1  0  9  0 . 0  0  0 
1  0  9  1  .  u  0  0 
1 092 . 000 
1 093 . UOO 
1  0  94 .  o  0  0 
1  0  9  5  .  0  0  o 


C  SUBROUTINE  P36TC  INCORPORATES  THE  OPTION  OF  836. 5(M)  AND 
C  CALCULATES  TONE  CORRECTIONS  ONLY  OVER  THE  RANGE  OF  BAND 
C  CENTER  FREQUENCIES  FROM  1000  HZ  (BAND  14)  TO  10,000  HZ 
C  CBAND  24). 

C 

C  THE  AIRCRAFT  NOISE  SPECTRA  WHICH  ARE  BEING  EXAMINED  ARE  KNOWN 

C  TO  HAVE  BAND  LEVELS  *HICH  HAVE  BEEN  SET  ARBITRARILY  TQ  0.0  DB 
C  BECAUSE  THE  SPL  SIGNAL  FROM  THE  AIRCRAFT  In  THAT  BAND  AT  SOME 
C  PARTICULAR  INSTANT  OF  TIME  WAS  NOT  MORE  THAN  5  DB  GREATER  THAN 
C  THE  CORRESPONDING  AMBIENT  NOISE  LEVEL  IN  THAT  BAND.  STRAIGHT- 
C  FORWARD  APPLICATION  OF  ThE  FAR  PART  36  TONE-CORRECTION  PROCEDURE 
C  WOULD  INDICATE  THAT  SPECTRA  WITH  MISSING  BAND  LEVELS  WOULD  HAVE 
C  PRONOUNCED  SPECTRAL  IRREGULARITIES,  EVEN  THOUGH  IT  *AS  NOT  BECAUSE 
C  OF  A  PURE-TONE  NOISE  SOURCE  ON  The  AIRCRAFT. 

C 

C  A  PROCEDURE,  ADOPTED  FROM  ThE  USAF  OMEGA  5.5  COMPUTER 

C  PROGRAM  FOR  ANALYZING  AIRCRAFT  FLYOVER  NOISE  MEASUREMENTS, 

C  IS  USED  TO  AVOID  SPURIOUS  TONE  CORRECTIONS  FOR  THOSE  SPECTRA 
C  WITH  MISSING  BAND  LEVELS.  BY  THIS  PROCEDURE,  ThE  SPL  SPECTRUM 
C  IS  SEARCHED  FOR  SPECTRAL  IRREGULARITIES  ONLY  OVER  A  LIMITED 
C  RANGE  OF  1/3-OCTAVE  BANDS.  THE  LIMITED  RANGE  IS  DETERMINED  BY 
C  STARTING  AT  THE  BAND  CONTAINING  THE  MAXIMUM  SPL  FOR  THE  SPECTRUM 
C  AND  WORKING  TOWARD  LOWER  FREQUENCIES  UNTIL  THE  SPL  BECOMES  LESS 
C  THAN  20.0  DB  AND  TriEN  WORKING  TOWARD  HIGHER  FREQUENCIES  UNTIL 
C  THE  SPL  AGAIN  BECOMES  LESS  THAN  20.0  DB .  THERE  MUST  BE  AT  LEAST 
C  6  CONSECUTIVE  SANDS  CONTAINING  SPLS  That  ARE  >  OR  =  20.0  DB  FOR 
C  ThE  TONE-CORRECTION  CALCULATION  TO  PROCEED.  IF  TnER E  ARE  FEWER 
C  THAN  6  BANDS  IN  This  RANGE,  THEN  THE  TONE -C OR RE C T I  ON  FACTOR  C 

C  IS  SET  =  0.0  DB  AND  RETURNED  TO  THE  CALLING  DR OG R AM  ALONG  WlTH 

C  J M X  = 1  FOR  THE  50-hZ  BAND  AS  A  WARNING. 

C 

C  USE  OF  A  LIMITED  FREQUENCY  RANGE  TO  DETERMINE  TONE-CGRRECT I  ON 

C  FACTORS  FOR  THE  SPL  SPECTRUM  AT  ANY  SPECIFIC  INSTANT  OF  TIME  IS 
C  DERIVED  FROM  SECTION  B36.5CL)  OF  PART  36  WHERE  PERMISSION  IS  GIVEN 
C  TO  USE  A  NARROW  BAND  ANALYSIS  TO  PROVE  That  A  CALCULATED  TONE 

C  CORRECTION  IS  NOT  THE  RESULT  OF  AN  AIRCRAFT  NOISE  SOURCE  IF  IT 

C  IS  SUSPECTED  That  DISCRETE-FREQUENCY  AIRCRAFT  NOISES  ARE  NOT 
C  ACT(jALLv  PRESENT.  THE  INTENT  OF  PARAGRAPH  CD  IS  ASSUMED  T0 
C  BE  TO  CREATE  TONE -CORREC T  I  ON  PENALTIES  ONLY  FOR  SPECTRAL 
C  IRREGULARITIES  RESULTING  FROM  D I S C R E T E -F RE Q uENC Y  AIRCRAFT 
C  NOISE  SOURCES. 

C 

C  NO  ATTEMPT  IS  MADE  TO  SYNTHESIZE  A  COMPLETE  24-BAND 

C  SPECTRUM  BY  ESTIMATING  SPLS  FOR  THOSE  BANDS  CONTAMINATED  BY 
C  HIGH  BACKGROUND  NOISE  LEVELS. 

C 

c 

DIMENSION  F(24),JCT(24),S(24),SP(25),SPL(24),S?LPC24),SPLPPC24) 

r 

C  SET  INITIAL  VALUES  FOR  The  TONE -C GR R EC t I  ON  FACTORS. 

C 

C  =  0  .  0  ;  CMX  =  0.0 

C 
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1096.000  C  SEARCH  THE  SPL  SPECTRUM  FROM  BAND  14  (1000  HZ)  TO  BAND  24 

1097.000  C  (10,000  HZ)  AND  IDENTIFY  THE  MAXIMUM  SPL  VALUE  AND  THE 

1098.000  C  CORRESPONDING  BAND  NUMBER  JM . 

1099.000  C 

1100.000  SPLMr-1 000 . 0 

1101. 000  DO  20  J - 1 4 , 24 

1  1  0  2.000  IF(SPL(J)-SPLM)20 ,  10,  10 

1103.000  10  JM  =  J 

1 IOh . 000  5PLM=SPL(U) 

1105.000  20  CONTINUE 

1106.000  C 

1107.000  C  NOW  SEARCH  TmE  SPL  SPECTRUM  FROM  SPLM  TO  HIGHER  FREQUENCIES 

1108.000  C  AND  TO  LOWER  FREQUENCIES  TO  FIND  THE  BANDS,  IF  ANY,  WHERE  THE 

1109.000  C  SPL  FIRST  BECOMES  LESS  THAN  20.0  DB  AND  NOTE  THE  CORRESPONDING 

1110.000  C  BAND  NUMBERS  UL1  AND  JH 1 . 

1111.000  C 

1112.000  DO  30  J  =  JM, 24 

1113.000  JH 1  =  J 

1  1  1  4.00  0  IF(SPL(J  )  .  LT  .  20 .  0  )GOTO  40 

1115.000  30  CONTINUE 

1115.000  GOTO  50 

1117.000  C 

1118.000  C  IF  A  BAND  LEVEL  HAS  BEEN  FOUND  THAT  IS  <  20.0  DB ,  THEN  RE-SET 

1119.000  C  THE  JH 1  COUNTER  TO  THE  NEXT  LOWER  VALUE  TO  IDENTIFY  THE  BAND 

1120.000  C  BEFORE  THE  BAND  WHERE  THE  LEVEL  IS  <  20.0  DB . 

1121.000  C 

1  1  2  2.00  0  40  JH  1  r  JH  1-1 

1123-000  C 

1124.000  5  C  DO  60  U  =  14,JM 

1125.000  C 

1125.000  C  CHANGE  TO  A  NEW  COUNTER  TO  SEARCH  THE  SPLS  FROM  THE  BAND  AT 

1127.000  C  JM  BACKWARDS  TO  LOWER  FREQUENCIES. 

1128.000  C 

1  1  29.00  0  JLUJM-J  +  14 

1  1  30.000  IF(SPL(JL  1  )  .  LT  .  20 . 0)GOTO  70 

1131.000  60  CONTINUE 

1132.000  GOTO  80 

1133.000  C 

1  1  34.00  0  C  IF  A  LOW-FREQUENCY  BAND  LEVEL  HAS  BEEN  FOUND  THAT  IS  <  20.0  DB  , 

1  1  3  5.00  0  C  RE-SET  THE  JL  1  COUNTER  TO  THE  NEXT  HIGHER  VALUE. 

1136.000  C 

1137.000  70  JL1=JL1+1 

1133.000  C 

1  1  39  -  00  0  C  J  L 1  AND  JH 1  ARE  INDICES  FOR  THE  RANGE  OF  1/3-OCTAVE  BANDS 

1140.000  C  WHERE  THE  SPLS  ARE  ALL  >  OR  =  20.0  DB .  TO  PROCEED  WITH  THE 

1  1  4  1.000  C  CALCULATION  OF  A  TONE -CORRECT  I  ON  F  A  C  TOR ,  THERE  MUST  BE  AT 

1  1  42.00  0  C  -EAST  6  3  AND  S  IN  THE  RANGE  FROM  JL 1  TO  JH  1  .  IF  (JH1-JL1-4) 

1  1  4  3.00  0  C  IS  <  1,  T-1  E  N  THERE  ARE  AT  MOST  ONLY  5  BAND  LEVELS  IN  THE 

1  1  44.00  0  C  RANGE  WHICH  ARE  >  OR  =  20.0  DB  . 

1145.000  C 

1  1  46.00  0  C  IN  THIS  CASE,.  A  TONE -CORREC T I  ON  FACTOR  OF  CMX  =  0 .  0  DB  IS 

1 147. 000  C  DECLARED  and  RETURNED  ALONG  WITH  THE  INDEX  JMX=1  (50  HZ) 

1148.000  C  AS  A  WARNING. 

1149.000  C 

1  1  5  0.0  0  0  3  0  I  f  (.  J  “  1  -  „  _  1  -  a  ,2  70,270,90 

1151.000  C 
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1152. 000 
1153. 000 
1154.000 
1155. 000 
1156.000 
1157.000 
1158. 000 
1159.000 
1160.000 
1161.000 
1162.000 
1163.000 
1164.000 
1165.000 
1166.000 
1167. 000 
1168. 000 
1169. 000 
1170. 000 
1171.000 
1172.000 
1173. 000 
1174.000 
1175.000 
1176.000 
1177.000 
1178. 000 
1179. 000 
1180.000 
1181. 000 
1182. 000 
1183.000 
1184. 000 
1 185. 000 
1186. 000 
1187.000 
1188. 000 
1189. 000 
1190.000 
1191. 000 
1192.000 
1193. 000 
1194.000 
1195. 000 
1196. 000 
1197.000 
1198. 000 
1199. 000 
1200. 000 
1201 . 000 
1202. 000 
1203. 000 
1204. 000 
1205. 000 
1206. 000 


C  FROM  HERE  TO  STATEMENT  NUMBER  120  CORRESPONDS  TO  STEPS  1,2,  AND 

C  3  OF  SECTION  836.5  OF  FAR  PART  36. 

C 

C  SET  INITIAL  VALUES  FOR  THE  FIRST  CALCULATION  OF  BAND-LEVEL 

C  SLOPE  S(J)  AND  OF  THE  FLAG  FOR  THE  BAND  INDEX  COUNTER  JCT(j). 

C 

90  S(JL1)=0.0 

JC  T ( JL  1  )=  0 
C 

C  SET  BAND  INDEX  COUNTERS  FOR  1  AND  2  MORE  THAN  JL1. 

C 

JL2=JL 1+1 
JL  3  =  UL  2  +  1 
JC  T ( JL  2 )  =  0 
C 

C  COMPUTE  VALUE  OF  SECOND  BAND  SLOPE. 

C 

S(JL2)=SPL(JL2)-SPL(JL1  ) 

C 

C  COMPUTE  BAND-LEVEL  CHANGES  (SLOPES)  FOR  REMAINDER  OF  THE 

C  FREQUENCY  RANGE  TO  BAND  INDEX  JH 1 .  SET  THE  JCT(J)  FLAG  =  1 

C  FOR  THOSE  SPLS  WHICH  ARE  TO  BE  ENCIRCLED  ACCORDING  TO  RULES 

C  OF  STEPS  2  AND  3  OF  B36.5  OF  FAR  PART  36. 

C 

DO  120  J=JL3,JH1 
JC  T (  J  )  =  0 

S(J)  =  SPL(J)-SPL(J-1  ) 

IF(ABS(S(J)-S(J-l))-5. 0)120,  120,  100 
100  IF(SCJ) .GT. 0. 0.AND. S(J).GT. S(J-1 ))GOTO  110 
IF(S (J ) . LE . 0 . 0. AND. SCJ-l )  .  GT  .  0  .  0  )  JC  T  (  J -1  )=  1 
GOTO  120 
110  JC  T ( J )=  1 
120  CONTINUE 
C 

C  FROM  HERE  TO  TWO  LINES  AFTER  STATEMENT  NUMBER  160  CORRESPONDS 

C  TO  STEPS  4  AND  5  OF  B36.5  OF  FAR  PART  36. 

C 

C  STEP  4-COMPUTE  THE  INITIAL  ESTIMATE,  SPLP(J),  FOR  THE  BROADBAND 

C  PART  OF  THE  NOISE  SPECTRUM. 

C 

SPLP(JL1)=SPL(JL1) 

DO  160  J=JL2, JH  1 
C 

C  TEST  THE  FLAG  FOR  THE  ENCIRCLED  BAND  LEVELS. 

C 

IF( JCT(J  ))1  30,  1  30,  140 
130  SPLP ( J )  =  S PL ( J  ) 

GOTO  160 
C 

C  SPECIAL  CALCULATION  FOR  SPLP  IF  THE  FLAG  IS  SET  FOR  THE 

C  LAST  BAND. 

C 

140  IFCJ.EQ. JHl)GOTO  150 
C 

C  COMPUTE  SPLP(J)  FOR  ENCIRCLED  BAND  LEVELS. 


51 


1207.000 
1208.000 
1209.000 
1210.000 
1211.000 
1212.000 
1213.000 
1219.000 
1215. 000 
1216. 000 
1217 . 000 
1218. 000 
1219-000 
1220 . 000 
1221 . 000 
1222.000 
1223 . 000 
1224. 000 
1225.000 
1226.000 
1227. 000 
1228.000 
1229-000 
1230.000 
1231 . 000 
1232.000 
1233-000 
1234.000 
1235.000 
1236.000 
1237.000 
1238. 000 
1239.000 
1240 . 000 
1241.000 
1242. 000 
1 243. 000 
1244. 000 
1245 . 000 
1246. 000 
1247. 000 
1 248 . 000 
1249 . 000 
1250 . 000 
1251.000 
1252.000 
1253.000 
1254. 000 
1255.000 
1256 . 000 
1257.000 
1258.000 
1259.000 
1260 . 000 
1261. 000 


C 

SPLP(J)  =  0  .  5;:(SPL(J-1  )+SPL(J  +  l  )) 

GOTO  160 

150  SPLP(J)=SPL(J-1 )+S(J-l ) 

C 

C  STEP  5-COMPUTE  NEW  BAND  LEVEL  SLOPES,  SP(J),  FOR  THE  SPLP(J) 

C  VALUES. 

C 

160  SP(J)  =  SPLP(J)-SPLP(J-1  ) 

C 

C  COMPUTE  VALUE  OF  SP  FOR  FIRST  BAND  AND  FOR  FICTITIOUS  BAND 

C  AFTER  THE  BAND  AT  JH  1  . 

C 

SP(JL1)=SP(JL2) 

SPCJH 1+1 )=SP(JH 1 ) 

C 

C  FROM  HERE  TO  STATEMENT  NUMBER  260  CORRESPONDS  TO  STEPS  6  TO  10 

C  IN  B  36 . 5  OF  FAR  PART  36. 

C 

C  COMPUTE  AN  INITIAL  VALUE  FOR  THE  FINAL  ESTIMATE,  SPLPP(J),  OF 

C  THE  BROADBAND  COMPONENT  OF  THE  NOISE  SPECTRUM. 

C 

SPLPP(JL1)=SPL(JL1) 

DO  260  J  =  JL  1 , JH  1 
C 

C  TEST  J  TO  SEE  IF  J=JHl.  USE  SPECIAL  RULE  TO  FIND  F(J)  FACTOR 

C  WHEN  J  =  JH  1  . 

C 

IF  C  J -JH  1  )1  70,  180,  180 

C 

C  STEP  6-COMPUTE  THE  ARITHMETIC  AVERAGE,  SA,  OF  THE  3  ADJACENT 

C  SP(J)  BAND  LEVEL  SLOPES. 

C 

170  SA=(SP(J)+SP(J+1 )+SP(J+2 ))/3  .  0 
C 

C  STEP  7-COMPUTE  REMAINDER  OF  VALUES  OF  SPLPP(J). 

C 

SPLPP ( J  +  l  )=SPLPP ( J )  +  SA 
C 

C  STEP  8-CALCULATE  F(J)  FACTORS  AS  THE  DIFFERENCE  BETWEEN  THE 

C  ORIGINAL  BAND  LEVEL,  SPL(J),  AND  THE  FINAL  ESTIMATE  OF  THE 
C  BROADBAND  COMPONENT  OF  THE  NOISE  SPECTRUM,  SPLPP(J). 

C 

180  F(J)=SPL(J)-SPLPP(J) 

C 

C  NOTE  THOSE  F(J)  VALUES  >  OR  =  0.0  DB .  NEGLECT  ANY  F(J)  <  0.0  DB . 

C 

IF(F( J))260, 190, 190 
C 

C  STEP  9-FIND  TONE-CORRECTION  FACTORS  C  WITH  SPECIAL  RULE  (TABLE 

C  B  2  OF  B36.5  OF  FAR  PART  36)  FOR  FREQUENCY  BANDS  14  TO  21  (1000 
C  TO  5000  HZ). 

C 

190  IF( J . GE . 14 . AND. J . LE . 21 )GOTO  220 
C 


I 
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1262.000  C  COMPUTE  TONE -CORRECT  I  ON  FACTOR  C  FOR  BAND  CENTER  FREQUENCIES 

1263.000  C  BETWEEN  6300  AND  10,000  HZ  AS  APPROPRIATE  TO  THE  RANGE 
1264.000  C  FROM  JL 1  TO  JH  1 . 

1265. 000  C 

1266.000  C  NOTE  SPECIAL  RULE  TO  LIMIT  C  IF  F(J)  >  20.0  DB . 

1267. 000  C 

1268.000  IF(F(J)-20. 0)200, 210, 210 

1269-000  200  C=F(d)/6.0 

1270.000  GOTO  250 

1271.000  210  C=3. 3333333 

1272.000  GOTO  250 

1273.000  C 

1274.000  C  COMPUTE  C  FOR  BAND  CENTER  FREQUENCIES  BETWEEN  1000  AND  5000  HZ. 

1275.000  C 

1276.000  220  IF ( F( J )-2 0 . 0)2 3 0 , 240 , 240 

1277.000  230  C=F(J)/3.0 

1278.000  GOTO  250 

1279.000  240  C=6. 6666667 

1280.000  C 

1281.000  C  STEP  10-COMPUTE  LARGEST  OF  C  VALUES  AND  STORE  IN  CMX. 

1282. 000  C 

1283.000  250  CMX=AMAX(CMX, C) 

1284. 000  C 

1285.000  C  TEST  CMX  AGAINST  C  TO  SET  THE  BAND  INDEX  JMX=J  WHEN  CMX=C . 

1286. 000  C 

1  28  7.000  IF  C  CMX . EQ . C )  JMX= J 

1288.000  260  CONTINUE 

1289.000  RETURN 

1290.000  270  CMX=0.0 

1291 . 000  C 

1292.000  C  IF  THE  CALCULATION  WAS  SENT  TO  STATEMENT  NO.  270,  SET 

1293.000  C  JMX=1  AS  A  WARNING  FLAG  THAT  THERE  WERE  NOT  ENOUGH  BANDS 
1294.000  C  WITH  SPL  >  20.0  DB . 

1295.000  C 

1296.000  JMX=1 

1297.000  RETURN 

1298.000  END 
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SUBROUTINE  INTEG(NS,  DATA,  DATAIO,  IF,  IL,  SUMDAT) 


Subroutine  INTEG  takes  the  number  of  samples,  NS,  the  array  of  input 
data  as  a  function  of  time  (in  dummy  array  DATA),  and  the  quantity  DATA10 
(a  number  which  is  10  dB  less  than  the  maximum  value  of  DATA)  and  deter¬ 
mines  the  energy  summation,  SUMDAT,  of  the  quantities  specified  by  DATA 
from  the  first,  IF,  to  the  last,  IL,  time  that  DATA  ■  DATA10.  The  DATA 
array  can  be  tone-corrected  perceived  noise  levels,  PNLT,  or  A-weighted 
sound  pressure  levels,  AL,  or  any  other  quantity  which  is  to  be  integrated 
over  defined  limits.  If  the  input  data  are  PNLT  or  AL  values,  then  DATA10 
is  PNLTM-10.0  or  ALM-10.0,  and  SUMDAT  is  used  'to  calculate  effective 
perceived  noise  level,  EPNL,  or  sound  exposure  level,  SEL. 

INTEG  requires  about  90  words  or  about  350  bytes  for  storage  of 
the  source  code. 


SUBROUTINE  I NTEGCNS , DA TA, DA TA  1 0 ,  IF,  I L,  SUMDAT) 


1299.000 
1300 . 000 
1  301 . 000 
1302 . 000 
1303. 000 
1304.000 
1305. 000 
1 306. 000 
1307. 000 
1308.000 
1309. 000 
1310. 000 
1311.000 
1312.000 
1313.000 
1314.000 
1315.000 
1516.000 
1317. 000 
1318.000 
1319.000 
1320. 000 
1321.000 
1322. 000 
1323. 000 
1324.000 
1325.000 
1326.000 
1327.000 
1328.000 
1329. 000 
1330. 000 
1331 . 000 
1332. 000 
1333.000 
1334.000 
1335. 000 
1336.000 
1337.000 
1338 . 000 
1339.000 
1 340 . 000 
1  341 . 000 
1342.000 
1343. 000 
1344.000 
1  345 . 000 
1346.000 
1347.000 
1348. 000 
1349. 000 
1350.000 
1351.000 
1352.000 
1353. 000 


C 


C 

C:i 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C:i 

c 

c 

c 

c 

c 


10 


c 

c 

c 

c 

c 


20 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  INTEG  DETERMINES  THE  INDEXES  IF  AND  IL  FOR  THE  FIRST 
AND  LAST  OF  THE  RELATIVE  TIMES  WHEN  THE  INPUT  DATA  ARE  EITHER 
EQUAL  OR  CLOSEST  TO  10-DB  LESS  THAN  THE  MAXIMUM  VALUE  OF  THE 
DATA  C I . E . ,  EQUAL  OR  CLOSEST  TO  DA  T  A  1 0  )  . 

SUMDAT  IS  THE  ENERGY  SUMMATION  OF  THE  DECIBEL  VALUES 
OF  THE  DATA  OVER  THE  TIME  INDEXES  FROM  IF  TO  IL. 

DATA  IS  A  DUMMY  ARRAY  CONTAINING  THE  INPUT  VALUES  IN 
DECIBELS  OF  PNLTCI)  OR  ALCl)  USED  FOR  CALCULATIONS 
OF  EPNL  OR  SEL,  RESPECTIVELY. 

NS  IS  THE  NUMBER  OF  DATA  SAMPLES  PASSED  8Y  THE  CALL  STATEMENT 
IN  THE  MAIN  PROGRAM. 


:;::c :::::: 


DIMENSION  DATA(l) 

DETERMINE  IF,  THE  RELATIVE  TIME  INDEX  COUNTER  WHERE  DATA(I) 

IS  FOR  THE  FIRST  TIME  GREATER  THAN,  OR  IS  CLOSEST  TO,  DA T A 1 0 . 

DO  10  1=1, NS 
IF  =  I 

IFCDATACI ).GT.DATA10)GOTO  20 
CONTINUE 


TEST  THE  VALUE  OF  IF,  OBTAINED  FROM  THE  ABOVE  SEARCH  FOR 
THE  FIRST  TIME  INDEX  COUNTER,  AGAINST  THE  INITIAL  TIME 
INDEX  1=1. 


I F C  IF -1  )40, 40, 30 


WHEN  IF  IS  >  l,  THEN  TEST  THE  VALUE 
SEE  IF  IT  IS  LESS  THAN  THE  VALUE  OF 
IF  IT  IS  NOT,  THEN  RE-SET  THE  VALUE 
VALUE  C I F  =  I F -1 )  TO  FIX  THE  STARTING 


OF  CDATA(IF)-DATAIO)  TO 
(DATA10-DATAC IF-1  ))  AND, 
OF  IF  TO  THE  NEXT  SMALLER 
TIME  FOR  THE  SUMMATION. 


THERE  ARE  THREE  POSSIBLE  CASES  WHEN  IF  IS  RE-SET  TO  IF-1: 

Cl)  WHEN  DATAC IF-1 )=DATA10;  (2)  WHEN  DATA(IF)  AND  DATA(IF-l) 
ARE  MORE  AND  LESS  THAN  DATA10  BY  EQUAL  AMOUNTS;  AND  (3)  WHEN 
DATAC IF-1 )  IS  LESS  THAN  DATA  10  BY  AN  AMOUNT  THAT  IS  LESS 
THAN  THE  AMOUNT  BY  WHICH  DATAC IF)  IS  GREATER  THAN  DATA  1 0 . 


THERE  ARE  TWO  POSSIBLE  CASES  WHEN  IF  IS  NOT  RE-SET  TO  IF-1: 
Cl)  WHEN  DATAC IF-1)  IS  LESS  THAN  DATA  10  BY  AN  AMOUNT  THAT 
IS  GREATER  THAN  THE  AMOUNT  BY  WHICH  DATAClF)  IS  GREATER 
THAN  DATA  10  Cl.E.,  WHEN  DATAClF)  IS  CLOSER  TO  DATA10 
THAN  IS  DATAC IF-1 );  AND  C2)  WHEN  IF=1  Cl.E.,  WHEN  THE 
FIRST  VALUE  OF  DATAC I )  =  DATAC 1 )  I S  =  OR  >  DATA  1 0 ) . 
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1354.000  C 

1355.000  30  IF((DATA(IF)-DATA10).LT.(DATA1G-DATA(1F-1 )))GOTO  40 

1  356.000  IF  =  I F -1 

1357.000  C 

1358.000  C  NOW  DETERMINE  IL,  THE  RELAT I VE -T I  ME  INDEX  COUNTER  WHERE 

1359-000  C  DATA(I)  IS  FOR  THE  LAST  TIME  GREATER  THAN,  OR  CLOSEST 

1360.000  C  TO,  DA T A 1 0 . 

1361 . 000  C 

1362.000  40  DO  50  1=1, NS 

1363.000  I L=NS-I +1 

1364. 0b0  IF(DATA(IL).GT.DATA10)GOTO  60 

1365.000  50  CONTINUE 

1366.000  C 

1367.000  C  TEST  THE  VALUE  OF  IL,  OBTAINED  FROM  THE  ABOVE  SEARCH  FOR 

1368.000  C  THE  LAST  TIME  INDEX  COUNTER,  AGAINST  THE  LAST  TIME 

1369.000  C  INDEX  WHEN  I=NS. 

1370.000  C 

1371.000  60  IF(IL-NS)70, 80, 80 

1372.000  C 

1373.000  C  WHEN  IL  IS  =  OR  >  NS,  THEN  USE  THE  VALUE  OF  IL  DETERMINED 

1374.000  C  BY  THE  ABOVE  SEARCH  SINCE  THIS  RESULT  MEANS  THAT  THE  LAST 

1  375.000  C  VALUE  OF  DATA( I )  =  DATA(NS )  I S  =  OR  >  DAT  A 1 0 . 

1376.000  C 

1377.000  C  WHEN  IL  IS  <  NS,  THEN  RE-SET  IL  TO  IL+1  TO  DETERMINE  THE 

1378.000  C  RELATIVE -TIME  INDEX  COUNTER  FOR  THE  VALUE  OF  DATA(I) 

1  379.000  C  WHICH  IS  EITHER  EQUAL,  OR  CLOSEST  TO,  THE  VALUE  OF  DATA  1 0 . 

1380.000  C 

1381.000  70  IFCCDATA(IL)-DATA10).LT.(DATA10-DATACIL+1 )))GOTO  80 

1382.000  I L= I L+l 

1383.000  C 

1384.000  C  NOW  DETERMINE  THE  SUMMATION,  ON  AN  ENERGY  BASIS,  OF  THE 

1385.000  C  INPUT  VALUES  OF  DATA(I)  IN  DECIBELS  BETWEEN  THE  10-DB-DOWN 

1386.000  C  TIMES  IF  TO  IL. 

1387.000  C 

1388  .  000  80  SUMDAT  =  0.0 

1  389.000  DO  90  I = I F ,  I L 

1  390.000  90  SUMDA  T  =  5  UMDA  T+l  0  .  0  ::;:(DATA(  I  )  / 1  0  .  0  ) 

1391.000  C 

1392.000  C  RETURN  TO  THE  MAIN  PROGRAM  WITH  THE  VALUES  OF  IF,  IL, 

1393.000  C  AND  SUMDAT. 

1394.000  C 

1395.000  RETURN 

1396.000  END 


This  page  intentionally  blank 


SUBROUTINE  MOLAR(TK,  HR,  PA,  WV) 

Subroutine  MOLAR  calculates  the  molar  concentration,  WV,  of  water  vapor 
in  the  air,  in  percent,  for  given  air  temperature,  TK,  relative  humidity, 

HR,  and  air  pressure,  PA. 

MOLAR  requires  about  85  words  or  about  350  bytes  for  storage. 


r 


SUBROUTINE  MOLAR  (TIC,  HR,  PA,  WV) 


1397 . 000  C 
1598.000 
1399. 000  C 
1400.000  C:: 

1401 . 000  C 

1402.000  C  SUBROUTINE  MOLAR  CALCULATES  WV,  THE  MOLAR  CONCENTRATION  OF 

1403.000  C  WATER  VAPOR  IN  THE  AIR  IN  PERCENT.  A  WV  VALUE  IS  CALCULATED 
1404.000  C  FOR  THE  AIR  TEMPERATURE,  RELATIVE  HUMIDITY,  AND  AIR  PRESSURE 
1405.000  C  SPECIFIED  BY  THE  CALLING  ROUTINE. 

1406. 000  C 

1407.000  C  TK  IS  AIR  TEMPERATURE,  DEG  KELVIN 

1408.000  C 

1409.000  C  HR  IS  RELATIVE  HUMIDITY,  PERCENT 
1410. 000  C 

1411.000  C  PA  IS  AIR  PRESSURE,  STANDARD  ATMOSPHERES 

1412.000  C 

1413.000  C  WV  IS  MOLAR  CONCENTRATION  OF  WATER  VAPOR,  PERCENT 
1414.000  C 

141  5.000  C :: :: :: J! :s !s « « JS !: 5! :: St s: :: :s :s :c ;j :s 
1416. 000  C 

1417.000  REAL  LPSOPO 

1  418.000  T  0  =  293 .  1  5 

1  419.000  T 0 1 =2  7 3 .  16 

1  420.000  LPSOPOzl  0  .  79586::(1  ■  -(T01  /TK))  -  5  .  02  8  0  8  -LOG  1  0  (  TK  /  T  0 1  ) 

142  1.000  1  +1  .  50474E-4”(1  . -1  0  .  - :s  C  -  8  .  29  69  2 ::(  (TK  /  T  0 1  )-l  .  ))) 

1  422.000  2  +0 . 428  7 3E -35{(  1 0  .  ::;:(4  .  7695  5::(  1  .  -(T01 /TK)))-1  .  ) 

1423.000  3  -2.2195983 

1  424.000  PSOP0  =  10.;:"  LPSOPO 

1  425.000  WV=HR::PSOPQ/PA 

1426.000  RETURN 

1427.000  END 


SUBROUTINE  ARP866A(J,  TC,  HR,  ALPHA) 


Subroutine  ARP866A  uses  the  procedure  of  SAE  AKP866A  (March  1975)  to 
calculate  the  pure-tone  atmospheric  absorption  coefficient,  ALPHA  in  dB/m, 
applicable  to  a  1/3-octave  band  specified  by  index  counter  J  and  for  given 
temperature,  TC,  and  relative  humidity  HR.  TC  is  in  degrees  Celsius  and 
HR  is  in  percent. 

ARP866A  requires  about  280  words  or  about  1100  bytes  for  storage. 


r 

1428.000 
1429 . 000 

c 

SUBROUTINE  ARP866A  ( J , TC , HR , ALPHA ) 

1430. 000 

c 

1431.000 

c 

1432. 000 

c 

* 

1433.000 

c 

J  IS  THE  INDEX  COUNTER  FOR  THE  GEOMETRIC  MEAN  FREQUENCY  OF 

1434. 000 

c 

THE  1/3-OCTAVE  BANDS  FROM  50  TO  10,000  HZ. 

- - - 

1435 . 000 

c 

1436. 000 

c 

TC  IS  AIR  TEMPERATURE,  DEGREES  CELSIUS. 

1437. 000 

c 

1458. 000 

c 

HR  IS  RELATIVE  HUMIDITY  OF  THE  AIR,  PERCENT. 

1439- 000 

c 

1440 . 000 

c 

ALPHA  IS  THE  VALUE  OF  THE  ATMOSPHERIC  ABSORPTION  COEFFICIENT  IN 

1441 . 000 

c 

DB/M  CALCULATED  BY  THE  METHOD  OF  SAE  ARP  866A,  MARCH  1975. 

1442 . 000 

c 

1443.000 

o 

j 

1444. 000 
1445 . 000 

c 

DIMENSION  FREQ(24),TA8LE(58) 

1446. 000 

DATA  TABLE  / 0 . 0 0 ,  0  .  0 0 ,  0 . 2 5 0  ,  0 . 3 1  5 ,  0 . 50  ,  0 . 70 ,  0 . 6 ,  0 . 84 ,  0 . 7 ,  0 . 9 3 , 

1447. 000 

10. 8, 0.975, 0.9, 0.996, 1.0, 1.0,  1.1, 0.97, 1.2, 0.9, 1.3, 0.84, 1.5, 0.75, 

1448. 000 

21. 7, 0.67, 2. 0,0. 57, 2. 3, 0.49  5, 2. 5, 0.45, 2. 8, 0.4, 3. 0,0. 37, 3. 30, 0.33, 

1449. 000 

33. 6, 0.3, 4.  15, 0.26, 4. 45, 0.24 5, 4. 80, 0.23, 5. 25, 0.22, 5. 7, 0.21, 

1450. 000 

46. 05, 0.20  5, 6. 5, 0.2, 7. 0,0. 2, 10. 0,0. 2/ 

1451 . 000 

DATA  FREQ  / 5  0 . ,  6 3  . ,  8 0 . ,  1 0 0  . ,  1 2 5 . ,  1  6 0 . ,  20 0 . ,  2 5 0  . ,  3 1  5 . ,  40 0 . ,  5 0 0 . , 

1452 . 000 

1630. ,800. ,1000. ,1250. ,1600. ,20 00. ,2500. ,3150. ,4000. ,4500., 

1453. 000 

25600 . , 7100. , 9000. / 

1454.000 

B  =  l. 328924- 3.  1  7  9  7  6  8E -0  2  ::TC +2  .  1 7  3  7  1  6E -0 4 "TC:! ;:2  .  -1  .  7496E-0  6:tTC::::3  . 

1455. 000 

ABSHUMzlO.  0:::S(LOG10(HR)-B) 

1456. 000 

AMOLMXrl  0.  0:!»(LOG10(FREQ(J))  +  8. 42 9 94E -0 3 “TC -2  .  755624) 

1457. 000 

HMOLMX=SQRT(FREQ(J)/1010. ) 

1458. 000 

HUMRAT  =ABSHUM/ HMOLMX 

1459- 000 

L  =2 

1460 . 000 

I  F  C HUMRAT . LE . TABLE ( 1 ))GOTO  40 

1461 . 000 

DO  10  1=3,57,2 

1462. 000 

L  =  I  +1 

1463. 000 

IF (TABLE  Cl  )-HUMRAT ) 1 0 , 40 , 20 

1464. 000 

10 

CONTINUE 

1465. 000 

GOTO  40 

1466. 000 

20 

1  =  1-2 

1467. 000 

IF  C I  .GE. 3)GOTO  30 

1468. 000 

1=1+2 

1 469 . 000 

30 

XAl=HUMRAT-TABLE(l ) 

1470. 000 

XA0=HUMRAT-TA8LE(I -2 ) 

1471.000 

XA2=HUMRAT-TABLE(I +2 ) 

1472. 000 

X0  1  =T  ABLE ( I -2 )-TABLE( I ) 

1473. 000 

X02=TABLE(I -2 )-TABLE(I +2 ) 

!  1474.000 

X 1 2  =TA8LE ( I )- TABLE (I +2 ) 

1475.000 

ALPRAT  =  TABLE(I  -1  ):!(XA1/X01  ):{(XA2/X0  2) 

1476.000 

1  -TABLE  (  I +1  ):-'(X  A0/X01  )5:(XA2/X12) 

|  1477.000 

2  +TABLE  (  I  +3  ):: ( XA  0 /X 0  2  )i;( XA  1  /X  1  2  ) 

1478. 000 

GOTO  50 

1479. 000 

40 

ALPRAT=TABLE(L) 

1480.000 

50 

alpmol=alprat«amolmx 

1481 . 000 

ALPCLA  =  1  0  .  Q::::(2  .  05  "LOG  10(FREQ(J)/1  000  .  )  +  1  .  1  3  9  4E  -  0  3  ::TC 

1  , 

1482. 000 

1  -  1.916984) 

;  1483.000 

ALPHA  =  0  .  01  ::(ALPCLA  +  ALPMOL) 

1  484 . 000 

RETURN 

1485. 000 

62 

END 

L 

m 

SUBROUTINE  ANSAB(TK,  HR,  PA,  F,  A) 

Subroutine  ANSAB  uses  the  procedure  of  American  National  Standard 
ANS  SI. 26-1978  to  calculate  the  pure-tone  atmospheric  absorption  of 
sound  A,  in  dB/m,  at  frequency  F,  in  Hz,  for  given  air  temperature  TK,  in 
degrees  kelvin,  relative  humidity  HR,  in  percent,  and  air  pressure  PA, 
in  standard  atmospheres. 

ANSAB  requires  about  200  words  or  about  850  bytes  for  storage. 


14-66.  000  C 

148  7.000  SUBROUTINE  ANSAB  (TK,HR,PA,F,A) 

1488. 000  C 

14  8  9  0  0  0  Q  <•  J*  I«  •»  <*  «*  <•  **  **  **  <>  '*  at  at  St  J»  ft  i(  at  at  at  at  at  ft  at  at  at  at  at  at  a!  at  at  at  at  at  JC  at  at  at  at  at  St  St  St  St  St  at  St  St  St  St  St  SC  St  5!  St  SC  SC  SC  SC  SC  SC  St  SC  SC 

1490. 000  C 

1491.000  C  SUBROUTINE  ANSAB  CALCULATES  PURE-TONE  SOUND  ABSORPTION 

1492.000  C  COEFFICIENTS  ACCORDING  TO  THE  PROCEDURE  OF  AMERICAN  NATIONAL 
1493.000  C  STANDARD  ANS  SI. 26-1978  FOR  THE  ATMOSPHERIC  ABSORPTION  OF 
1  494.  000  C  SOUND . 

1495.000  C 

1496.000  C  TK  IS  AIR  TEMPERATURE,  DEGREES  KELVIN. 

1497 . CFO  C 

1498.000  C  HR  IS  THE  RELATIVE  HUMIDITY  OF  THE  AIR,  PERCENT. 

1499. 000  C 

1500.000  C  PA  IS  THE  PRESSURE  OF  THE  AIR,  STANDARD  ATMOSPHERES. 

1501.000  C 

1502.000  C  F  IS  THE  FREQUENCY  OF  THE  SOUND,  HZ. 

1503.000  C 

1504.000  C  A  IS  THE  PURE-TONE  ATMOSPHERIC  SOUND  ABSORPTION  COEFFICIENT, 

1505.000  C  D8/M. 

1506.000  C 

1507.000  C  TO  AND  T01  ARE  REFERENCE  TEMPERATURES,  DEGREES  KELVIN. 

1508.000  C 

1  5  0  9.00  0  C  L  PS  OP  0  IS  THE  LOGARITHM  OF  THE  RATIO  OF  THE  VAPOR  PRESSURE 

1510.000  C  AT  SATURATION  TO  THE  STANDARD  ATMOSPHERIC  PRESSURE,  I.E.,  THE 
1511.000  C  LOGARITHM  OF  PSOPO. 

1512.000  C 

1513.000  C  H  IS  THE  MOLAR  CONCENTRATION  OF  WATER  VAPOR,  PERCENT. 

1514.000  C 

1  5  1  5.000  C  FR02  AND  FRN  2  ARE  THE  RELAXATION  FREQUENCIES  OF  OXYGEN 

1516.000  C  AND  NITROGEN,  RESPECTIVELY,  HERTZ. 

1517.000  C 

1518.000  C  ALPHA  IS  THE  ATMOSPHERIC  ABSORPTION  COEFFICIENT, 

1519. 000  C  NEPERS/M. 

1520.000  C 

1  5  2  1.000  C "  " "  ■"  - ::  "  "  " "  -  " :: :: "  " " ::  -  " ::  " "  »  "  "  -  "  "  "  -  "  -  -  -  -  »  "  "  " "  "  "  -  "  "  "  "  » » "  -  - -  " ::  - 
1522.000  C 

1523. 000  REAL  LPSOPO 

1  5  24.000  T  0  =  2  9  3 .  15 

1525.000  101=273.16 

1  526  .  00  0  LPSOP0  =  10.  79586::  Cl  .-(T01/TK))  -  5  .  02808::LOG1  0CTK/T01  ) 

1  52  7.000  1  +1  .  50474E-4::(1  . -1  0  .  8  .  29  69  2  ::C  (  TK  /  T  0 1  )- 1  .))) 

1  5  28.000  2  +  0, 42873E-3"(10."»(4.  76955::(  1  . -(T01/TK)))-1  .) 

1529.000  3  -2.2195983 

1  5  30.000  PSOP0  =  l  0  .:::;LPSOP0 

1  53  1.000  H  =HR  ::PS  OP  0  /  PA 

1  5  3  2.000  FR02  =  PA"C24.  +4. 41 E  04  ::H  "(  (  0  .  0  5 +H  )  /  (  0  .  39  1 +H  )  )  ) 

1533.000  FRN2=(PA/SQRT(TK/T0))"(9.+ 

1534.000  1  350."H"EXPC-6. 142"CC(TK/T0)""(-1 ./3.))-l .))) 

1535.000  ALPHA=(F""2 . ) - ( C ( 1 . 8 4E -1 1 )"(1 . / PA ) "S QR T C TK / T 0 ) ) 

1  5  36.000  1  +CCTK/T0)""C-5./2.))"CCl.  278E-25: 

1537.000  2  (EXPC-2239. 1/TK))/(FR02+CCF""2 . )/FR02))) 

1538.000  3  +(0. 1068-CEXPC-3352 . / TK ) ) / ( FRN 2+ ( ( F- "2 ,)/FRN2))))) 

1  5  39.000  A=  8 . 68  6 "ALPHA 

1540.000  RETURN 

1541.000  END 
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SUBROUTINE  NUMINT  (SPL,  TKR,  HRR,  PAR,  ATK,  AHR,  PA,  D,  KL,  BA) 


Subroutine  NUMINT  calculates  the  array  of  band  adjustment  factors  BA,  in 
decibels,  appropriate  for  the  set  of  measured  test-day  1/ 3- octave- band  sound 
pressure  levels  specified  in  dummy  array  SPL.  The  reference  atmospheric 
conditions  (given  by  parameters  TKR,  HRR,  and  PAR)  and  the  test-day  atmos¬ 
pheric  conditions  as  a  function  of  height  above  ground  level  (given  by 
arrays  ATK  and  AHR)  and  the  uniform  pressure  PA  are  specified  by  the  calling 
program  as  is  the  array  D  of  pathlength  segments  along  the  sound  propaga¬ 
tion  path  from  the  microphone  back  to  the  aircraft  at  its  position  on  the 
flight  path  at  the  time  it  emitted  the  sound  measured  by  the  microphone 
as  the  sound  pressure  levels  SPL. 

Subroutine  NUMINT  uses  a  numerical  integration  technique  to  integrate 
the  product  of  an  equivalent,  normalized  sound  pressure  spectral  density 
function  and  an  absorption  function  over  the  sound  propagation  path  for 
the  two  atmospheric  conditions  and  thereby  calculate  a  band  adjustment 
factor  for  ideal  1/3-octave-band  filters. 

Subroutine  DIFFS  Is  called  to  calculate  a  set  of  band-level  differences, 
or  slopes,  for  the  input  SPL  array.  The  differences  are  used  in  the 
approximation  of  the  sound  pressure  spectral  density  function  over  the 
bandwidth  of  each  filter. 

Subroutine  ANSAB  is  called  to  determine  a  pure-tone  sound  absorption 
coefficient  at  specified  frequencies  over  the  filter  bandwidth  and  over  the 
sound  propagation  path. 

Subroutine  QSF  is  calleu  to  perform  the  actual  numerical  integration, 
at  equally  spaced  frequency  intervals,  of  the  calculated  integrand  values 
over  the  bandwidth  of  the  ideal  filter. 

NUMINT  requires  about  470  computer  words  or  about  1900  bytes  for 
storage. 
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1342.000  C 

1543.000  SUBROUTINE  NUM1NT(SPL,TKR,HRR,PAR,ATK,AHR,PA,D,KL,BA) 

1544. 000  C 

1 546 .000  C 

1547.000  C  SUBROUTINE  NUMINT  USES  A  NUMERICAL  INTEGRATION  TECHNIQUE 

1548.000  C  TO  CALCULATE  AN  ADJUSTMENT  FACTOR  TO  BE  APPLIED  TO  1/3- 

1549.000  C  OCTAVE-BAND  SOUND  PRESSURE  LEVELS  TO  ACCOUNT  FOR  DIFFERENCES 

1 550 . 00C  C  IN  THE  ATMOSPHERIC  ABSORPTION  OF  SOUND  OVER  A  SEGMENTED  SOUND 
1551.000  C  PROPAGATION  PATH  UNDER  TWO  DIFFERENT  SETS  OF  METEOROLOGICAL 

1552.000  C  CONDITIONS,  I.E.,  TEST  AND  REFERENCE  CONDITIONS. 

1553.000  C 

1554.000  C  SPL  --  THE  ARRAY  OF  24  MEASURED  1 / 3 -OC TA VE -BAND  SPLS 
1555.000  C  WITH  NOMINAL  BAND  CENTER  FREQUENCIES  BETWEEN  50  AND 

1556.000  C  10,000  HZ,  DB. 

1557.000  C 

1558.000  C  TKR  --  THE  REFERENCE  AIR  TEMPERATURE,  DEG  KELVIN. 

1559.000  C 

1560.000  C  HRR  --  THE  REFERENCE  RELATIVE  HUMIDITY  OF  THE  AIR,  PERCENT. 
1561.000  C 

1562.000  C  PAR  --  THE  REFERENCE  ATMOSPHERIC  PRESSURE,  STANDARD  ATMO- 
1563.000  C  SPHERES. 

1564.000  C 

1565.000  C  ATK  --  THE  ARRAY  OF  AVERAGE  AIR  TEMPERATURES,  OVER  EACH 

1566.000  C  SEGMENT  OF  THE  SOUND  PROPAGATION  PATH,  MEASURED  AT 

1567.000  C  THE  TIME  THE  MEASURED  SPLS  WERE  OBTAINED,  DEG  KELVIN. 

1568.000  C 

1569.000  C  AHR  --  THE  ARRAY  OF  AVERAGE  RELATIVE  HUMIDITIES,  OVER  EACH 

1570.000  C  SEGMENT  OF  THE  SOUND  PROPAGATION  PATH,  MEASURED  AT  THE 

1571.000  C  TIME  THE  MEASURED  SPLS  WERE  OBTAINED,  PERCENT. 

1572.000  C 

1573.000  C  PA  --  THE  ATMOSPHERIC  PRESSURE  AT  THE  TIME  THE  MEASURED 
1574.000  C  SPLS  WERE  OBTAINED,  STANDARD  ATMOSPHERES. 

1 575. 000  C 

1576.000  J  D  --  THE  ARRAY  OF  DISTANCES  ALONG  THE  SEGMENTS  OF  THE  SOUND 

1  5  7  7  .  00  0  C  PROPAGATION  PATH  CORRESPOND  I NG  TO  THE  HORIZONTAL  LAYERS 

1578.000  C  DEFINED  BY  THE  INTERVALS  OF  HEIGHT  WHERE  THE  AIR 

1579.000  C  TEMPERATURE  AND  RELATIVE  HUMIDITY  WERE  MEASURED,  METERS. 

1580.000  C 

1581.000  C  KL  --  THE  INDEX  COUNTER  FOR  THE  ARRAYS  ATK,  AHR,  AND  D  THAT 
1582.000  C  INDICATES  THE  LAST  INCREMENT  ALONG  THE  PATH  FROM 

1583.000  C  THE  MICROPHONE  TO  THE  AIRPLANE. 

1584.000  C 

1585.000  C  BA  --  THE  BAND-ADJUSTMENT  FACTOR  TO  BE  ADDED  TO  THE 
1586.000  C  MEASURED  SPL  TO  DETERMINE  THE  SPL  THAT  WOULD  HAVE  BEEN 

1587.000  C  MEASURED  UNDER  REFERENCE  METEOROLOGICAL  CONDITIONS,  DB . 

1588.000  C 

1  589.000  C::::::  SUBROUTINE  NUMINT  CALLS  SUBROUTINE  DIFFS  TO  DETERMINE  THE  BAND 

1590.000  C  LEVEL  DIFFERENCES,  OR  DATA  SLOPES,  IN  DB/BAND  OVER  THE 

1591.000  C  LOWER  AND  UPPER  HALVES  OF  EACH  OF  THE  24  FILTER  BANDS 

1592.000  C  CORRESPONDING  TO  THE  BAND  LEVELS  IN  THE  ARRAY  SPL. 

1593.000  C 

1  594.000  C:::: ::  SUBROUTINE  NUMINT  ALSO  CALLS  SUBROUTINE  ANSAB  TO  DETERMINE 
1595.000  C  PURE-TONE  ATMOSPHERIC  ABSORPTION  COEFFICIENTS,  IN  DB/M, 

1596.000  C  AS  A  FUNCTION  OF  FREQUENCY  FOR  SPECIFIED  AIR  TEMPERATURE, 

1597.000  C  RELATIVE  HUMIDITY,  AND  AIR  PRESSURE. 
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1598.000  C 

1599.000  C!!5!:s  THE  ACTUAL  PROCESS  OF  NUMERICALLY  INTEGRATING  THE  ARRAY 
1500.000  C  OF  CALCULATED  INTEGRAND  VALUES  OVER  THE  BANDWIDTH 

1601.000  C  OF  THE  1/3-OCTAVE-BAND  FILTERS  (ASSUMED  TO  BE  IDEAL 

1602.000  C  FILTERS)  IS  PERFORMED  BY  A  STANDARD  ROUTINE  DESCRIBED 

1603.000  C  IN  SUBROUTINE  QSF  WHICH  IS  ALSO  CALLED  BY  SUBROUTINE 

1604. 000  C  NUMINT. 

1605.000  C 

1606.000  C::;::J  NOTE  THAT  A  VARIETY  OF  COMPARABLE  NUMERICAL  INTEGRATION 


1607.000  C  METHODS  COULD  BE  SUBSTITUTED  FOR  SUBROUTINE  QSF  WITH 
1608.000  C  APPROPRIATE  MODIFICATIONS  TO  THE  PROGRAM  STATEMENTS. 

1609. 000  C 

16  1  0.000  c ::  55 :: :t  5: ::  5! :: :: :: :: :: :: :: :: ::  5: 5: :: :: :: :: :: ::  :s :: :: :: ::  :s ::  5:  :: :: :: :: " ::  5: :: « :s « 5: ::  ”  »  -  -  - 

1611.000  C 

161  2.000  DIMENSION  AHR  (  3  1  )  ,  ATK  (  3  1 ),  BA(  2  4  )  ,  D(  3  1  )  ,  NSTEP  (  2  4  )  , 

1613.000  IS LOPE (24, 2), SPL(2 5 ) , YNL( 1 6 ) , YNU( 1 6 ) , ZNL ( 1 6 ) , ZNU ( 1 6 ) 

1614.000  C 

161  5  .  000  DATA  NSTEP  / 1  3  ::4,  5 ,  6 ,  7 ,  8 ,  9 ,  1 0 ,  1 1 ,  1 2 ,  1 3 ,  1 4,  1 5  / 

1616. 000  C 

1617.000  C  THE  DIMENSIONS  OF  THE  ARRAYS  HAVE  BEEN  SET  TO  SPECIFIC 

1618.000  C  VALUES  TO  AGREE  WITH  THE  CORRESPONDING  ARRAYS  IN  THE 
1619.000  C  CALLING  ROUTINE,  OR,  IN  THE  CASE  OF  THE  Y  AND  Z  ARRAYS, 
1620.000  C  TO  1  MORE  THAN  THE  LARGEST  NUMBER  OF  FREQUENCY  STEPS,  NSTEP, 
1621.000  C  OVER  THE  WIDEST  BANDWIDTH  OF  EACH  HALF  OF  THE  IDEAL 
1622.000  C  FILTER  AT  BAND  24  (10,000  HZ). 

1623.000  C 

1624.000  C  DEFINE  RF,  THE  EXACT  VALUE  OF  THE  FREQUENCY  RATIO 

1625.000  C  FOR  1/3-OCTAVE  BANDS. 

1626.000  C 

1527.000  RF  =1  0 . 05t::0  .  1 

1628 . 000  C 

1629.000  C  DETERMINE  FC,  THE  EXACT  BAND  CENTER  FREQUENCY,  IN  HZ, 

1630.000  C  FROM  THE  DEFINITION  BASED  ON  THE  SET  OF  INTERNATIONAL 
1631.000  C  STANDARD  BAND  NUMBERS,  ISBN,  STARTING  WITH  ISBNrl7  FOR 
1632.000  C  THE  FIRST  BAND,  J=l,  WITH  NOMINAL  CENTER  FREQUENCY  OF  50  HZ. 
1633.000  C 

1634.000  DO  30  J=l,24 

1655 . 000  C 

1636.000  I SBN=1 6+U 

1  637.000  FC  =  10  .  05:::(I  SBN/10.  0) 

1638.000  C 

1639.000  C  CALCULATE  THE  EXACT  LOWER,  FL,  AND  UPPER,  FU,  BANDEDGE 

1640.000  C  FREQUENCIES  FOR  THE  FILTER  BAND. 

1641 . 000  C 

1642.000  C  CALCULATE  THE  FREQUENCY  INTERVALS  FOR  THE  LOWER,  DELFL, 

1643.000  C  AND  UPPER,  DELFU,  HALVES  OF  THE  FILTER  BAND. 

1644. 000  C 

1645.000  FL  rFC / SQR  T(RF) 

1646.000  FU=F  C::SQR  T(RF  ) 

1647. 000  C 

1648.000  DELFL=(FC-FL)/NSTEP(J) 

1649.000  DELFU=(FU-FC)/NSTEP(J) 

1650. 000  C 

1651.000  C  ASSIGN  THE  VALUES  FOR  THE  FREQUENCIES  FI  AND  F2  BEFORE 

1652.000  C  THE  START  OF  THE  NUMERICAL  INTEGRATION  OVER  THE  LOWER 
1653.000  C  AND  UPPER  HALVES  OF  THE  FILTER  BAND. 
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165^.000  C 

1  65  5.00  0  FUFL-DELFL 

1656.000  F2=FC-OELFU 

1657. 000  C 

1658.000  C  CALCULATE  NN,  THE  TOTAL  NUMBER  OF  STEPS  (OR  FREQUENCY 

1659.000  C  INTERVALS)  IN  THE  NUMERICAL  INTEGRATION  OVER  EACH  HALF 
1660.000  C  OF  A  FILTER  BAND. 

1661.000  C 

1662.000  NN=NSTEP(J)+1 

1665. 000  C 

1664.000  C  FIND  THE  BAND-LEVEL  SLOPES,  IN  DB/BAND,  CORRESPONDING 

1665.000  C  TO  THE  INPUT  SPLS. 

1666. 000  C 

1667.000  CALL  DIFFSCSPL, SLOPE  ) 

1668.000  C 

1669.000  C  CALCULATE  SPSL  AND  SPSU  FOR  THE  ESTIMATED  SLOPES  OF  THE 

1670.000  C  NORMALIZED  SOUND  PRESSURE  SPECTRAL  DENSITY  FUNCTIONS  OVER 
1671.000  C  THE  LOWER  AND  UPPER  HALVES  OF  THE  FILTER  BAND  BY  USING  THE 
1672.000  C  CORRESPONDING  BAND-LEVEL  SLOPES. 

1673. 000  C 

16  7  4.00  0  SPSL  =  S  LOPE  (  J ,  0/(10  .  0  -"LOG  1  0  (R  F  )  ) 

1  67  5  .  00  0  SPS  U  =S  LOPE ( J , 2)/ (10. 0 -LOG  1 0 (RF ) ) 

1676. 000  C 

1677.000  C  TEST  THE  VALUE  OF  THE  SLOPES  TO  DETERMINE  WHICH  FORM  TO 

1678.000  C  USE  FOR  CALCULATING  THE  DENOMINATOR  TERMS  FOR  THE  LOWER  HALF 
1  6  79.000  C  OF  THE  BAND,  DENI,  AND  THE  UPPER  HALF  OF  THE  BAND,  DEN  2 . 

1680. 000  C 

1  68  1.000  IF(SPSL.NE.  0.  0)DEN1  =  (FC/SPSL)::(1  .  -  (RF:: ::  (  -  S  PS  L  /  2  .))) 

1  68  2.00  0  IF  (SPSL.  EQ.  0.  0  )D  EN  1  =  (  FC  /  2  .  )::LOG(RF) 

1685.000  C 

1684.000  IF(SPSU.NE.  0.  0)DEN2  =  (FC/SPSU)::((RF””(SPSU/2  .  ))-l  .  ) 

1685.000  IF(SPSU.  EQ.  0.  0)DEN2=(FC/2  .  )::LOG(RF) 

1686.000  C 

1687.000  C  DETERMINE  THE  ARRAY  OF  INTEGRAND  VALUES  FOR  THE  NUMERATOR 

1688.000  C  TERMS  OVER  THE  LOWER  AND  UPPER  HALVES  OF  A  FILTER  BAND. 

1689 . 000  C 

1690.000  DO  20  JJ=1,NN 

1691.000  C 

1692.000  C  DETERMINE  THE  FREQUENCY  AT  EACH  STEP  IN  THE  LOWER,  FI, 

1693.000  C  AND  UPPER,  F2,  HALF  OF  THE  BANDWIDTH. 

1694. 000  C 

1695.000  F1=F1+DELFL 

1  696.000  F  2=F  2+DELF  U 

1597.000  C 

1698.000  C  DETERMINE  THE  VALUE  OF  THE  NORMALIZED  SOUND  PRESSURE 

1699.000  C  SPECTRAL  DENSITY  FUNCTION  AT  FREQUENCIES  FI  AND  F2. 

1700. 000  C 

17  0  1.00  0  XL  =  (F1/FC)::::(-1  .+SPSL) 

1  702  .  0  00  XU  =  (F2/FC):;::(-1  . +SPSU ) 

1703.000  C 

1704.000  C  DETERMINE  THE  ATMOSPHERIC  ABSORPTION  COEFFICIENTS  FOR 

1705.000  C  THE  REFERENCE  METEOROLOGICAL  CONDITIONS  AT  FREQUENCIES 
1706.000  C  FI  AND  F2  WITH  THE  ASSUMPTION  OF  A  UNIFORM  REFERENCE  ATMOSPHERE. 


1707. 000  C 

1708.000  CALL  ANSABC TKR , HRR , PAR , F  1 , AR  1 ) 

1709.000  CALL  AN SA8( TKR , HRR , PAR , F 2 , AR 2 ) 

1710. 000  C 

17  1  1.000  C  DETERMINE  THE  EXPONENT,  AE,  OF  THE  ATMOSPHE R I C -AB SOR PT I  ON 

1712.000  C  ATTENUATION  FACTOR,  AF ,  BY  SUMMING  THE  DIFFERENCES  IN 
1713.000  C  ATTENUATION  UNDER  TEST  AND  REFERENCE  METEOROLOGICAL  CONDITIONS 
1714.000  C  OVER  THE  SEGMENTS  ALONG  THE  LENGTH  OF  THE  SOUND  PROPAGATION 
1 715. 000  C  PATH. 

1716. 000  C 

1717.000  AEL=AEU=0 . 0 

1718.000  DO  10  K=1,KL 

1  7  1  9  .  000  CALL  ANSABC ATK(K ) , AHR (K ) , PA , F 1 , A 1 ) 

1  720  .  000  CALL  ANSABC ATK « ) , AHR (K ) , PA , F 2, A  2 ) 

1721.000  AEL=AEL+CA1-AR1)“DCK) 

1  722  .  000  1  0  AEU=AEU  +  CA2-AR2):!DCK) 

1723. 000  C 

1  724.000  AFL  r  1  0  .  0  ::::  CAEL/10.  0) 

1  72  5  .  000  AFUrlO.  0:”"CAEU/10.  0) 

1726. 000  C 

1727.000  C  DETERMINE  THE  VALUES  OF  THE  NUMERATOR  INTEGRANDS  AS  THE 

1728.000  C  PROOUCTS  OF  THE  NORMALIZED  EQUIVALENT  SOUND  PRESSURE  SPECTRAL 
1  729  .  000  C  DENSITY  FUNCTIONS  AND  THE  ATMOSPHER I C -ABSORPT  I  ON  FACTORS. 
1730.000  C  STORE  THE  RESULTS  IN  ARRAYS  YNL  AND  YNU . 

1731.000  C 

1  7  32  .  000  YNLCJJ)  =  XL::AFL 

1  733.000  20  YNUCJ«J)=XU"AFU 

1734.000  C 

1735.000  C  CALL  THE  NUMERICAL  INTEGRATION  SUBROUTINE  QSF  TO 

1736.000  C  INTEGRATE  THE  INTEGRAND  ARRAYS  OVER  THE  LOWER  AND  UPPER 
1737.000  C  HALVES  OF  THE  FILTER  BANDWIDTH. 

1738. 000  C 

1739.000  CALL  QSFCDELFL, YNL, ZNL, NN ) 

1  740.000  CALL  QS FCDELF U , YNU , ZNU , NN  ) 

1741 . 000  C 

1742.000  C  THE  LAST  VALUE  AT  NN  OF  THE  OUTPUT  ARRAYS  ZNL  AND  ZNU 

1743.000  C  IS  THE  RESULT  OF  THE  NUMERICAL  INTEGRATION  OVER 
1  744.000  C  THE  RANGE  OF  FREQUENCIES  IN  EACH  HALF  OF  THE  FILTER  E AND  . 

1745. 000  C 

1746.000  C  DETERMINE  THE  BAND  ADJUSTMENT  FACTOR,  BA,  FROM 

1747.000  C  THE  RESULTS  OF  THE  NUMERICAL  INTEGRATIONS. 

1748. 000  C 

1749.000  30  BACJ)=10.  05:LOG1  0  CC  ZNLCNN  )  +  Z  NUCNN  )  )  /  CDEN  1 +DEN  2  )  ) 

1750. 000  C 

1751.000  RETURN 

1752.000  END 
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SUBROUTINE  DIFF$(SPL,  SLOPE) 


Subroutine  DIFFS  accepts  an  array  of  24  1/ 3-octave-band  sound  pressure 
levels  and  calculates  the  spectral  slopes  (SLOPE)  or  band-level  differences 
over  the  lower  and  upper  halves  of  each  band.  Special  rules  are  defined 
for  the  first  and  last  bands  and  for  handling  bands  where  the  indicated 
sound  pressure  level  has  been  set  equal  to  0.0  dB  because  the  signal  level 
was  not  more  than  5.0  dB  greater  than  the  corresponding  background  noise 
level . 

DIFFS  requires  about  170  words  or  about  700  bytes  for  storage  of  the 
source  code. 
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SUBROUTINE  DIFFSCSPL, SLOPE) 


17  53. 

.  000 

1754. 

,  000 

1  7  55  . 

000 

1  756. 

000 

1  757  . 

000 

1758. 

.  000 

1  759. 

000 

1  760. 

000 

1761. 

000 

1  762  . 

000 

1  763. 

000 

1764. 

000 

1  765. 

000 

1766. 

000 

1  767. 

000 

1  768. 

000 

1  769. 

000 

1  770. 

000 

1771. 

000 

17  72  . 

000 

1  773. 

000 

1  774. 

000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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000 
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SUBROUTINE  DIFFS  COMPUTES  THE  DIFFERENCES  IN  THE  LEVEL  OF 
THE  1 /3-OCTAVE -BAND  SOUND  PRESSURE  LEVELS  OVER  THE  LOWER  AND 
UPPER  HALVES  OF  EACH  OF  THE  24  BANDS  WITH  CENTER  FREQUENCIES 
FROM  50  HZ  TO  10,000  HZ.  THE  LOWER  HALF  OF  A  BAND  IS 
DEFINED  TO  BE  FROM  THE  LOWER  BANDEDGE  FREQUENCY  TO  THE 
EXACT  GEOMETRIC  MEAN  FREQUENCY;  THE  UPPER  HALF  IS  FROM  THE 
GEOMETRIC  MEAN  FREQUENCY  TO  THE  UPPER  BANDEDGE  FREQUENCY. 

THE  FILTER  BANDS  ARE  ASSUMED  TO  HAVE  IDEAL  TRANSMISSION 
RESPONSE  CHARACTERISTICS. 

SPL  THE  ARRAY  OF  24  SOUND  PRESSURE  LEVELS,  IN  DB , 

PASSED  BY  THE  CALLING  PROGRAM  PLUS  A  LEVEL 
OF  0.0  DB  FOR  A  FICTITIOUS  BAND  25. 

SLOPE  THE  ARRAY  OF  LOWER  AND  UPPER  BAND-LEVEL 

DIFFERENCES,  OR  SLOPES,  CALCULATED  FOR  EACH 
OF  THE  24  BANDS,  IN  DB/BAND. 


DIMENSION  SLOPE ( 24, 2),SPL(25),SS(24) 

DATA  SS/3::l.,13::0.,-l.,-2.,-3.,-4.,-5.,-7.,-10.,-15./ 

SPL(25)=0  .  0 

DO  70  J  =  1 ,23 
IFCU.LT. 23)GOTO  10 


C 

C  SPECIAL  RULE  FOR  BAND  24  (1  0,  000  HZ)  WHEN  SPL(24)  IS  >  0.0  DB , 

C  BUT  SPL (2  3 )  IS  <  OR  =  0.0  DB . 

C 

IF(SPL(24).GT.0. 0 . AND . SPL (2 3 ) . LE . 0. 0)GOTO  80 
C 

C  ANOTHER  SPECIAL  RULE  FOR  BAND  24  WHEN  SPL(23)  IS  >  0.0  DB 

C  BUT  SPLC24)  IS  <  OR  =  0.0  DB . 

C 

IF(SPL(23) .GT. 0. 0.AND. SPL(24) .LE. 0. 0)GO  TO  90 
C 

C  ANOTHER  SPECIAL  RULE  FOR  BAND  24  WHEN  SPL(23)  AND  SPL (24) 

C  ARE  BOTH  <  OR  =  0.0  DB . 

C 

IF(SPL(23).LE. 0. 0. AND. SPL (2 4) .LE. 0. 0)GOTO  100 
C 

C  IF  SPL(J)  IS  <  OR  =  0.0  DB,  THEN  IN  ORDER  TO  HAVE  SOME 

C  SPECIFIC  VALUE  TO  USE  IN  SUBSEQUENT  CALCULATIONS  SET  THE  SLOPE 
C  FOR  THIS  BAND  ARBITRARILY  TO  +1  DB/BAND. 

C 

C  A  SLOPE  OF  +1  DB/BAND  REPRESENTS  A  WHITE-NOISE  SPECTRUM  AND, 

C  WITH  SUBROUTINE  NUMINT,  YIELDS  A  MEASURE  OF  THE  "AVERAGE" 

C  ATTENUATION  OVER  THE  FREQUENCY  RANGE  OF  THE  1/3-OCTAVE  BAND. 
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1809. 000  C 
1810. 000 
1811.000 
1812. 000  C 
1813.000  C 
1814. 000  C 
1815. 000  C 
1816.000 
1817.000 
1818.000  C 
1819.000  C 
1820. 000  C 
1821.000  C 
1822.000  C 
1823.000 
1824. 000  C 
1825.000  C 
1826. 000  C 
1827. 000  C 
1828 . 000 
1829. 000  C 
1830. 000  C 
1831.000  C 
1832. 000  C 
1833. 000  C 
1834. 000  C 
1835. 000  C 
1836.000  C 
1837. 000  C 
1838.000  C 
1839 . 000 
1840. 000  C 
1841 . 000  C 
1842. 000  C 
1343.000 
1844.000  C 
1845 . 000  C 
1846. 000  C 
1847. 000  C 
1848 . 000 
1849. 000 
1850. 000  C 
1851.000  C 
1852. 000  C 
1853. 000  C 
1854. 000 
1855. 000 
1856.000  C 
1857.000  C 
1858. 000  C 
1859- 000  C 
1860. 000  C 
1861.000  C 
1862. 000  C 
1863.000 


10  IF(SPL(J).LE. 0. 0)GOTO  60 
IFCJ.GT.  DGOTO  20 

SPECIAL  RULE  FOR  BAND  1  (50 
BUT  SPLC2)  IS  <  OR  =  0.0  DB  . 

IF(SPL(2).LE. 0. 0)GOTO  50 
IFCJ.EQ.  DGOTO  30 


HZ)  WHEN  SPL(i)  IS  >  0.0  DB , 


SPECIAL  RULE  FOR  THE  CASE  WHEN  J  IS  BETWEEN  2  AND  23 
AND  SPL(U)  IS  >  0.0  DB,  BUT  SPL(J-l)  AND  SPL(J+1)  ARE 
BOTH  <  OR  =  0.0  DB . 

20  IFCSPLCJ-l ).LE. 0. O.AND. SPLCJ+l ).LE. 0. 0)GOTO  50 

PROCEED  WITH  NORMAL  CALCULATION  OF  BAND-LEVEL  SLOPES 
WHEN  SPL(J)  AND  SPL(J+1)  ARE  BOTH  >  0.0  DB . 


IF(SPL(J+1 ) .GT. 0. 0)GOTO  30 

IF  SPL(J)  >  0.0  DB,  BUT  SPL(J+1)  I S  <  OR  =  0.0  DB,  THEN 
SKIP  SLOPE  CALCULATION  BY  NORMAL  PROCEDURE  BECAUSE  THIS  CASE 
IS  TREATED  AS  THOUGH  THE  BAND  WHERE  SPL(J)  I S  >  0 . 0  DB  I S 
THE  LAST  BAND  WITH  VALID  DATA  AND  THE  SEARCH  STARTS  OVER 
AGAIN  AT  SPLCJ+2).  THE  SLOPE  FOR  THE  UPPER  HALF  OF  THE  BAND 
AT  J  IS  EXTRAPOLATED  FROM  THE  SLOPE  CALCULATED  FOR  THE 
LOWER  HALF.  THE  ARBITRARY  SLOPE  OF  +1  DB/BAND  IS  CALCULATED 
FOR  THE  BAND  AT  J  +  l  IN  THIS  CASE. 

GOTO  70 

FIRST,  FIND  SLOPE  OVER  UPPER  HALF  OF  BAND  AT  J. 

30  SLOPECU, 2)=SPL(J+1 )-SPL(J) 

SLOPE  OVER  LOWER  HALF  OF  BAND  AT  J+l=SLOPE  OVER  UPPER 
HALF  OF  BAND  AT  d. 


SLOPE (J+l , 1 )= SLOPE (J, 2) 
I F ( J . GT . l)GOTO  40 


SPECIAL  RULE  FOR  SLOPE  OVER  LOWER  HALF  OF  BAND  1  WHEN 
SPL(l)  AND  SPL (2 )  ARE  BOTH  >  0.0  DB . 

SLOPE ( 1 ,  l)=SLOPE(l, 2) 

GOTO  70 


SPECIAL  RULE  FOR  THE  SLOPE  OVER  THE  UPPER  HALF  OF  THE 
BAND  AT  J+l  WHEN  SPL(J)  AND  SPL(J+1)  ARE  BOTH  >0.0  D3 , 

BUT  THE  SPL(J+2 )  I S  <  OR  =  0.0  DB .  THIS  RULE  ALSO  APPLIES 
TO  BAND  24  WHEN  SPL(23)  AND  SPL(24)  ARE  BOTH  >  0 . 0  DB 
SINCE  SPL (25)  IS  ARBITRARILY  SET  =  0.0  DB . 

40  IF(SPL(J+2).LE. 0. 0)S LOPE (J+l, 2 )= SLOPE (J+l , 1) 
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1864. 000  C 

1865.000  C  SPECIAL  RULE  FOR  THE  SLOPE  OVER  THE  LOWER  HALF  OF  THE  BAND 

1866.000  C  AT  J  IF  SPL(J-l)  I S  <  OR  =  0.0  DB  AND  J  IS  >  1. 

1867.000  C 

1868.00  0  IFCSPLCJ-O.LE.  0. 0)SLOPE(J,  1  )=  SLOPE  (J,  2) 

1869.000  GOTO  70 

1870.000  C 

1871.000  C  FOR  THOSE  BANDS  WHERE  SPL  IS  >  0.0  DB  BUT  SPL(J-l)  ANO 

1872.000  C  SPLCJ+l)  ARE  BOTH  <  OR  =  0.0  DB ,  USE  THE  ARBITRARY  SUBSTITUTE 
1873.000  C  BAND-LEVEL  SLOPES  SS(J)  FROM  THE  DATA  STATEMENT. 

1874.000  C 

1875.000  50  SLOPECJ, OrSLOPECU, 2)=SSCU) 

1876.000  GOTO  70 

1877.000  C 

1878.000  C  SPECIAL  RULE  FOR  BANDS  WHERE  SPL(U)  I S  <  OR  =  0.0  DB . 

1879. 000  C 

1880.000  60  SLOPECJ, 1 )=SLOPE (U , 2)=1 . 0 

1881.000  C 

1882.000  70  CONTINUE 

1883.000  C 

1884.000  RETURN 

1385.000  C 

1886.000  C  SPECIAL  RULE  FOR  BAND  24  WHEN  SPL(24)  IS  >  0.0  DB 

1  887.000  C  BUT  SPLC2 3 )  I S  <  OR  =  0.0  DB . 

1888. 000  C 

1889.000  80  SLOPEC23, O=SL0PEC23, 2)=1 . 0 

1890.000  SLOPEC24, l)=SLOPE(24, 2)=SS(24) 

1891.000  RETURN 

1892.000  C 

1893.000  C  SPECIAL  RULE  FOR  BAND  24  WHEN  SPL(2  3)  IS  >  0.0  DB , 

1894.000  C  BUT  SPL(24)  I S  <  OR  =  0.0  DB . 

1895.000  C 

1896.000  90  SLOPEC24, l)=SLOPE(24, 2)=1 . 0 

1897.000  RETURN 

1898. 000  C 

1899.000  C  SPECIAL  RULE  FOR  BANDS  23  AND  24  WHEN  BOTH  SPL( 2  3 )  AND 

1900.000  C  SPLC24)  ARE  <  OR  =  0.0  DB . 

1901.000  C 

1  90  2.00  0  100  SLOPEC2  3, l)=SLOPE(23, 2)=SLOPEC24, l)=SLOPE(24, 2)=1 . 0 

1903.000  RETURN 

1 904. 000  C 
1905.000  END 
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SUBROUTINE  QSF ( H ,  Y,  Z,  NDIM) 


Subroutine  QSF  numerically  integrates  the  array  of  function  values 
Y  to  determine  the  array  of  integral  values  Z  at  the  specified  increment 
H  for  the  argument  of  Y.  NDIM  is  the  dimension  of  arrays  Y  and  Z.  A 
combination  of  Simpson's  and  Newton's  rules  is  used  to  perform  the  calcu¬ 
lations.  The  QSF  program  is  taken  from  the  IBM  Scientific  Subroutine 
Package.  It  is  included  here  for  completeness.  Other  numerical  integra¬ 
tion  methods  could  be  used  in  place  of  QSF  with  appropriate  modifications 
to  Subroutine  NUMINT. 

QSF  requires  about  310  words  or  about  1250  bytes  for  storage. 
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1906. 000 
1907. 000 
1908 . 000 
1909. 000 
1910. 000 
1911.000 
1912. 000 
1913. 000 
1914. 000 
1915.000 
1916. 000 
1917 . 000 
1918.000 
1919.000 
1920.000 
1921.000 
1922.000 
1923. 000 
1924. 000 
1925 . 000 
1926.000 
1927.000 
1928. 000 
1  929  -  000 
1950.000 
1931.000 
1932.000 
1933.000 
1934.000 
1935.000 
1936. 000 
1957. 000 
1 938 . 000 
1939. 000 
1940.000 
1941.000 
1 942 . 000 
1943.000 
1 944. 000 
1945. 000 
1946 . 000 
1947. 000 
1948. 000 
1 949. 000 
1950. 000 
1951.000 
1952. 000 
1953.000 
1954. 000 
1955. 000 
1956. 000 
1957. 000 
1958.000 
1959. 000 
1960.000 
1961.000 
1962.000 
1963.000 
1  964.  00  0 


C 

c 
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c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 
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SUBROUTINE  QS F (H , Y , Z , NO  I M) 

»» *►  *»  •  *  <* »» !» • » <•  * » Ji »» •»  it  •«  St  i!  it  it  <C  !t  n  »t  It  SC  SC  St  «C  St  «C  SC  SC  SC  SC  SC  «C  SC  St  SC  i%  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  SC  £ 

SUBROUTINE  QS F  FROM  IBM  SCIENTIFIC  SUBROUTINE  PACKAGE  (SSP), 
FIFTH  EDITION, AUGUST  1970. 

PURPOSE 

TO  COMPUTE  THE  ARRAY  OF  INTEGRAL  VALUES  FOR  A  GIVEN 
EQUIDISTANT  TABLE  OF  FUNCTION  VALUES. 

USAGE 

CALL  QS  F  (H, Y ,  Z  , ND I M) 

DESCRIPTION  OF  PARAMETERS 

H  -THE  INCREMENT  OF  ARGUMENT  VALUES. 

Y  -THE  INPUT  ARRAY  OF  FUNCTION  VALUES. 

2  -THE  RESULTING  ARRAY  OF  INTEGRAL  VALUES.  7  MAY  BE 
IDENTICAL  WITH  Y. 

ND I M  -THE  DIMENSION  OF  ARRAYS  Y  AND  Z. 

REMARKS 

NO  ACTION  IN  CASE  ND I M  LESS  THAN  3. 

NO  SUBROUTINES  OR  FUNCTION  SUBPROGRAMS  ARE  REQUIRED. 

METHOD 

BEGINNING  WITH  ZCl)=0,  EVALUATION  OF  ARRAY  Z  IS  DONE  BY 
MEANS  OF  SIMPSON'S  RULE  TOGETHER  WITH  NEWTON'S  3/8  RULE  OR 
A  COMBINATION  OF  THESE  TWO  RULES.  TRUNCATION  ERROR  IS  OF 
ORDER  H::::5  (I.  E., FOURTH  ORDER  METHOD).  ONLY  IN  CASE  NDIM=3 
IS  THE  TRUNCATION  ERROR  OF  Z( 2)  OF  ORDER  H«*4. 


C 


DIMENSION  Y(1),Z(1) 

H T -0  .  33  3 3 3 3 3  ::H 
lF(NDIM-5)7, 8, 1 

ND I M  IS  GREATER  THAN  5.  PREPARATION  OF  INTEGRATION  LOOP. 

1  SUM1=Y(2)+Y(2) 

SUM  1 =SUM 1+SUM 1 

SUM 1=HT“(Y(1 )+SUMl+Y(3)) 

AUX 1 =Y (4 )  + Y (4  ) 

AUX1=AUX1+AUX1 

AUX1  =SUM1+HT”(Y(3 )  +  AUXl+Y(5 )) 

AUX2=HT::(Y  Cl  )  +  3  .  875:c(Y(2)  +  Y(5  ))  +  2  .  62  5  ::C  Y  (  3  )  +  Y  C4  ))+Y  (6  )  ) 
SUM  2  =  Y ( 5 )  +  Y C  5  ) 

S  UM  2  =  S  UM  2  +  S  UM  2 

S  UM  2  =  AUX  2  -HT  ::C  Y  C4  )+SUM  2+Y  C6  )  ) 

Z( 1 )=0  .  0 
AUX  =  Y  (.  3  )+Y  (  3  ) 

AU  X  =A  UX  +AU  X 

ZC  2  )  =  SUM2-HT-"(Y(2)+AUX+Y(4)) 

ZC  39-SUM1 
Z( 4 )  =  S UM  2 
|F(NDIM-6)5, 5, 2 
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INTEGRATION  LOOP. 


1965. 000  C 
1966. 000  C 
1967.000  C 
1968.000  2 

1969.000 
1970. 000 
1971 . 000 
1972 .000 
1973 . 000 
1974.000 
1975 . 000 
1976.000  3 

1977. 000 
1978 . 000 
1979.000  4 

1980.000  5 

1981 . 000 
1982 . 000 
1983.000  6 

1984. 000 
1985. 000 
1986.000  C 
1987 . 000  C 


DO  4  1=7,NDIM, 2 
SUM  1 =AUX 1 
SUM2=AUX2 

AUX1=Y(1-1)+Y(I-1) 

AUX1=AUX1+AUX1 

AUXlrSUMl+HT"CY(I -29+AUX1+YCI )) 

ZC 1 -2)=SUM1 
I F ( I -ND I M) 3 , 6, 6 
AUX2=Y(I )+Y(I ) 

AUX2=AUX2+AUX2 

AUX2  =  SUM  2+HT  ::(  Y  (  I  -1  )+AUX2+Y C I +1 ) ) 

Z( I -1 )=SUM2 
ZCNDIM-1 )=AUX 1 
Z( ND I M) =AUX2 
RETURN 

ZCNDIM-1 )  =  S  UM  2 

ZCNDIM)=AUX1 

RETURN 

END  OF  INTEGRATION  LOOP  FOR  ND I M  GREATER  THAN  5. 


1988.000 

C 

1989 • 000 

7 

IF  (ND I M-3 ) 1 2 , 11,8 

1990 . 000 

C 

1991.000 

c 

NDIM  IS  EQUAL  TO  4  OR  5. 

1992 . 000 

c 

1993.000 

8 

SUM  2  =  1  .  1  2  5  ::H T  Y  (  1  )+Y(2)  +  Y(2)  +  Y(2)  +  Y(3)  +  Y(3)  +  Y(3)  +  Y(4)) 

1994. 000 

SUM  1  =  Y(2 )+Y(2 ) 

1995.000 

SUMUSUM1+SUM1 

1996.000 

SUM1=HT::(Y(1  )  +  SUMl+Y(3  )) 

1997. 000 

ZC 1 )  =  0  .  0 

1998.000 

AUX 1 =Y  (3  )  +  Y  (3  ) 

1 999 .000 

AUX1=AUX1+AUX1 

2000.000 

Z(2)  =  SUM2-HT::(Y(2  9+AUX1 +Y (4 ) ) 

2001.000 

IF(NDIM-5)10,9,9 

2002.000 

C 

2003 . 000 

c 

HERE  NDIM  IS  EQUAL  TO  5. 

2004. 000 

c 

2005.000 

9 

AUX1=Y(4)+Y(4) 

2006.000 

AUX1=AUX1+AUX1 

2007.000 

Z(5)  =  SUM1+HT::(Y(3)+AUX1+Y(5)) 

2008.000 

C 

2009 . 000 

C 

HERE  NDIM  IS  EQUAL  TO  4. 

2010.000 

C 

2011 . 000 

10 

Z( 3 )  =  SUM  1 

2012.000 

Z(  4  }  =  S  UM  2 

2013-000 

RETURN 

2014. 000 

C 

2015. 000 

C 

HERE  NDIM  IS  EQUAL  TO  3. 

2016.000 

c 

2017.000 

1 1 

SUM  1  =HT::(  1  .  2 5  "Y  C  1  )  +  Y  (2  )  +  Y  (2  )- 0  .  25:;Y(3)) 

2018. 000 

S  UM  2=Y(2)  +  Y(29 

2019.000 

S  UM  2  =  S  UM  2+S  UM  2 

2020.000 

Z  (  3  )  =  HT;:(Y  C  1  9  +  SUM2  +  Y  (3  )) 

2021 . 000 

ZC1)=0 • 0 

2022.000 

Z(2  )  =  SUM1 

2023. 000 

12 

RETURN 

2024.000 

END 

4.  SAMPLE  INPUT  DATA  FILE 


The  next  four  pages  present  an  input  data  file  in  punched-card- image  format 
for  use  as  a  test  case  with  the  program  statements  of  Section  3.  The  test  case 
selected  was  run  378  from  the  DC-9  data  analyzed  in  Volume  I. 

The  data  are  stored  in  the  order  and  with  the  format  that  they  are  read  by 
program  TESTREF.  For  the  test  case  shown,  there  are  211  lines  of  data  arranged 
as  follows: 

line  nos.  description 

1  RUN  stored  as  13 

2  MIC,  NS,  DATE,  MONTH,  YEAR  and  stored  as  II,  14,  13,  1A4 , 

13 

3  TS,  TOH,  AH,  ZM,  AS,  AM,  GAMAD  and  stored  as  F7.1,  F8.1, 
F6.1,  F4.1,  F5.1,  F6.3,  F5 . 1 

4  to  143  SPL ( I , J)  where  I  -  1,  NS  and  J  »  1,24  and  stored  two  lines 

at  a  time  as  12F6.1  for  NS  data  samples  (here  NS  ■  70  for 
34.5  seconds  of  data  starting  at  0.0  seconds) 

144  NH  stored  as  13 

145  to  175  HM(K)  where  K  ■  1,  NH  and  stored  as  F6.1  for  NH  data 
samples  (here  NH  -  31) 

176  to  206  TC(K) ,  HR(K)  where  K  -  1,  NH  and  stored  as  2F5.1  for  the  NH 
data  samples 

207  BP  stored  as  F6.3 

208  TITLE1  stored  as  20A4 

209  TITLE2  stored  as  20A4 

210  GROUND  stored  as  5A4 

211  999  for  the  run  number  used  to  indicate  the  end  of  the 
data  and  to  signal  the  program  to  STOP  execution  and  END 
calculations 

The  symbols  for  the  various  data  items  described  above  are  those  used 
in  the  program  statements  and  are  defined  in  Section  2. 


EXAMPLE  INPUT  DATA  FILE 


line 

number  Data 

1.000  378 

2.000  4  70  25  OCT  74 


3 

.  000 

44184.0 

44199 

.  7 

158 

.  0 

1 .  2 

72 

.  6 

0.212 

0.  i 

0 

4 

.  000 

,  0 

.0 

.  0 

.  0 

.  0 

.  0 

.  0 

.0 

50  , 

.  3 

51 

.  5 

.  0 

.0 

5 

.  000 

.  0 

51 

.  5 

52 

.  0 

49 

.  1 

47 

8 

46. 

1 

43. 

,  8 

44. 

8 

44. 

3 

.  0 

,  0 

.  0 

6  . 

.  000 

0 

.  0 

.0 

.0 

.  0 

.  0 

.  0 

.  0 

50. 

3 

.  0 

48 

.  0 

50 

.  6 

7  , 

,  000 

0 

50 

.  6 

50 

.6 

.  0 

46 

8 

44. 

6 

43. 

.  8 

44. 

8 

44. 

3 

0 

.  0 

.0 

8. 

,  000 

0 

.0 

.  0 

.  0 

,  0 

.  0 

.  0 

,  0 

.0 

,  0 

50 

1 

52 

.  5 

9. 

,  000 

48 

8 

50 

.  3 

51 

.  3 

48 

,  1 

48 

,  1 

43. 

,  3 

43. 

.  8 

44. 

.  8 

44  . 

3 

.  0 

.0 

.  0 

10. 

000 

0 

.  0 

.  0 

.  0 

.  0 

.  0 

.  0 

.  0 

.  0 

52 

.  6 

55 

.  1 

55 

.  1 

1 1 . 

,  000 

48 

8 

50 

,  6 

52 

.  0 

53 

0 

51 

1 

44. 

1 

43. 

,  8 

44. 

8 

44  . 

3 

,  0 

0 

.  0 

12  . 

000 

0 

.0 

.  0 

.  0 

,  0 

0 

.  0 

52. 

8 

51 . 

5 

52 

.  6 

54 

.6 

53 

.  8 

13. 

.  000 

0 

51 

.  3 

51 

.  3 

53 

0 

51 

.  1 

46. 

6 

43. 

,  8 

44. 

8 

44. 

3 

.  0 

.  0 

.  0 

14. 

.  000 

0 

.  0 

.  0 

.  0 

.  0 

.  0 

.  0 

52. 

.  6 

52. 

6 

52 

.6 

53 

.  0 

52 

.  1 

15. 

.  000 

0 

49 

.  8 

52 

.  0 

54 

.  3 

55 

.  8 

50. 

.  5 

46. 

.  1 

44  . 

,  8 

44  . 

3 

.  0 

.  0 

.  0 

1  8  . 

,  000 

0 

.  0 

.  0 

.0 

.  0 

.  0 

54, 

.  8 

55  , 

.  8 

52  . 

8 

52 

.  6 

52 

.  8 

51 

.  1 

17. 

.  000 

51 

8 

53 

.  0 

55 

.  8 

57 

.  0 

56 

.  8 

53. 

.  8 

46 

.8 

44  , 

.  8 

44  . 

.  3 

.  0 

.  0 

.  0 

13. 

.  000 

0 

.0 

.  0 

.  0 

.  0 

.  0 

59 

.  1 

59  . 

.  6 

57. 

.  1 

54 

.6 

49 

.  8 

49 

.  8 

19. 

.  000 

51 

1 

54 

.6 

55 

.  3 

55 

,  1 

56 

3 

50. 

.  8 

46, 

.  3 

44  . 

.  8 

44  . 

3 

.  0 

.  0 

.  0 

20, 

.  000 

0 

.  0 

.  0 

.  0 

.  0 

.  0 

59. 

.  3 

61  , 

.  0 

59. 

.  0 

54 

.  0 

51 

.  1 

.  0 

21. 

.  000 

0 

54 

.  3 

53 

.6 

54 

,  1 

55 

.  8 

50. 

,  5 

46, 

.  1 

44  , 

.8 

44. 

3 

.  0 

.  0 

.  0 

22. 

.  000 

0 

.  0 

.0 

.  0 

.  0 

59. 

,  3 

61 , 

.  8 

59. 

.  8 

58. 

8 

56 

.  6 

52 

.  3 

50 

.  1 

23. 

.  000 

51 

8 

59 

•  5 

59 

.  1 

53 

.6 

55 

.  1 

51  . 

.3 

47 

.6 

44  , 

.  8 

44  . 

.  3 

.0 

.0 

.0 

24. 

.  000 

0 

.  0 

.  0 

.  0 

.  0 

59. 

.  8 

63 

.  0 

62  , 

.  8 

61. 

,  8 

58 

.  5 

54 

.6 

51 

.  5 

25. 

.  000 

53 

8 

60 

.  8 

62 

.  3 

58 

,  0 

58 

.  3 

54. 

.0 

51 

.  0 

46 

3 

44. 

3 

.0 

.0 

.0 

26  , 

.  000 

0 

.  0 

.0 

.  0 

.  0 

60. 

,  8 

64 

.5 

65 

.3 

62  . 

.  8 

59 

.  1 

56 
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5.  SAMPLE  OUTPUT  LISTING 


The  next  eigh'-  pages  contain  the  listings  of  output  data  that  resulted 
from  running  program  TESTREF  in  Section  3  with  the  input  data  in  Section  4. 
The  output  data  are  in  the  order  and  with  the  format  specified  by  the  WRITE 
and  FORMAT  statements  in  TESTREF. 

The  first  five  of  the  eight  pages  are  called  the  data  pages.  They 
contain  the  test-day  sound  pressure  levels  [read  from  input-data  array 
SPL(I,J)]  at  0.5-s  intervals  for  the  1/3-octave  bands  with  center  frequencies 
from  50  to  10,000  Hz.  They  also  contain  the  test-day  PNL,  PNLT,  and  AL 
values.  The  center  frequency,  FMX  in  kHz,  of  the  1/ 3-oc cave-band  sound 
pressure  level  that  produced  the  largest  tone-correction  factor  is  also 
listed  on  the  data  pages.  The  magnitude  of  the  tone-correction  factor 
can  be  seen  by  subtracting  the  PNL  from  the  PNLT  value.  A  value  of  .0  for 
a  sound  pressure  level  means  that  the  aircraft  noise  signal  at  that  time 
for  chat  band  was  not  more  than  5.0  dB  greater  than  the  corresponding 
background  noise  level  (ambient  noise  +  instrument  noise) . 

The  fifth,  or  final,  data  page  also  contains  the  values  of  the  tesc-dav 
PNLM,  PNLTM,  and  ALM  along  with  their  relative  times  of  occurrence.  The 
value  of  the  duration-correction  factor,  DCF,  for  the  test-day  EPNL  is 
given  with  the  relative  times  over  which  the  PNLT  integration  was  performed 
to  find  EPNL.  Test-day  EPNL  is  listed  on  the  line  below  that  with  PNLM, 
PNLTM,  and  DCF.  The  last  line  on  the  fifth  page  contains  ALM,  SEL, 
and  the  relative  times  over  which  the  AL  values  were  integrated  to  find 
SEL. 


The  sixth  page  lists  the  test-day  and  reference-day  meteorological 
data  at  the  heights  above  ground  level  corresponding  to  the  midpoints  of 
the  layers  of  the  atmosphere  defined  by  the  heights  in  input  data  array 
HM(K) .  At  the  foot  of  the  listings  is  the  station  barometric  pressure; 
the  air  temperature,  relative  humidity,  and  molar  concentration  of  water 
vapor  interpolated  at  a  height  of  10.0  m  from  the  actual  measured  data; 
and  the  corresponding  uniform-atmosphere  values  for  acoustical-reference-dav 
cond it  ions . 


r 

\ 


The  final  cwo  pages  contain  the  last  two  pages  of  summary  output  data. 
The  seventh  page  contains  the  test-day  SPLs  at  the  time  of  occurrence  of 
Ph’LTM,  the  corresponding  reference-day  SPLS  adjusted  by  the  four  methods, 
and  the  corresponding  PNL,  PNLT,  and  EPNL  values.  The  eighth  page  contains 
similar  data  except  at  the  time  of  occurrence  of  ALM  with  corresponding 
reference-day  values  of  AL  and  SEL.  At  the  foot  of  each  page  is  the  height 
of  the  airplane  over  the  microphone  and  the  sound  propagation  pachlength 
and  directivity  angle  at  the  time  of  the  test-day  PNLTM  or  ALM. 
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